j  REPORT  DOCUMENTATiON  PAGE 

A  PH  ^  A  n  tUMM  to  I  Mr  ttr  mmemt 

AD-A250  566  \7z:^r::.z:zn:':ss^s:z^ 

m^\m^  a  Mmm'^  w  ^  CW«>*<tHnmnwn»KlyML< 


OCUMENTATiON  PAGE  vS 

It  nnMM  IS  I  itssr  tsf  istsoft.  Mousmt  ■•"S  ler  'wywt  nomiwwt.  issronws  tiitmy  mu  tu»i«i. 

I  twwpuws  ins  mutmn  nu uissniew  o<  •Mo»m*uos.  i*"S  nysKS^  out  swstn  ttumsn  or  »r*T  o«or  tiaoci  o<  u« 

lor  iiSuBwt  ion aufsso.  le  wo— miQw  rnosiMUtort  torwak  QironoriU  lot  wniroumu  OsoftiwM  orts  MeeHL  HM  roftomo 
a>M.  MU  IS  itic  OtHcoo<ini<nii«itiiM<lMSoM.>»sor»iO«i  »«*»«lso>ro>«tlie>»«-aHB.w»tMr»9»ti.  OC  2S4SJ. 


i  2.  Rf  PORT  OATI  TX  MlPOnT  TY«  AMD  OATIS  CCVtMtD 


April,  1992 


Final 


Soot  Particle  Inception  and  Growth  Processes  in 
Combustion 


t.  AUTMOUtf) 

Robert  J.  Santoro 


7.  PlMfi 


MMMG  ORGANIZATION  NAME(S)  ANO  AOORfSSttS) 

The  Pennsylvania  State  University 
240  Research  Building  E,  Bigler  Road 
University  Park,  PA  16802 


1/15/87-1/15/92 


S.  PUNOMG  NUMREIIS 

PE  -  61102P 
PR  -  2308 
SA  -  BS 

AFOSR  -  87  -  0145 


•.  REVORMMC  ORGANIZATION 
RIRORT  NUMRCR 


AFOSR-TR.  cJ  2  0  3  96 


I.  SRONSORMG/ MONITORING  AGENCY  NAME(S}  ANO  i 

AFOSR/ NA 
Building  410 

Bolling  AFB  DC  20332-6448 


1C 

ELECTE 
MAYi  8  1992 


10.  SRONSORviG/ MONITORING 
A6INO  REKMT  NUMBER 


12*.  OISTminOH/AV 


TATiMINT 


12b.  OlSTIUBUTiON  CODE 


Approved  for  public  release;  dlscribucion  Is 
unllmiced 


12.  abstract  (MuMTMmZOOuroraW 


A  Study  of  soot  particle  inception  and  growth  has  been  completed  which  considered  fuel  molecular  structure, 
fuel  concentration,  temperature  and  operating  pressure  effects.  These  studies  indicated  that  fuel  species  most 
strongly  affected  the  particle  inception  process,  as  opposed  to  the  surface  growth  process,  and  support  an 
interpretation  that  inception  controls  the  maximum  amount  of  soot  formed.  Studies  of  concentration  and 
temperature  variations  indicated  that  temperature  effects  dominate.  The  analysis  yielded  an  apparent 
activation  energy  of  94.5  kcal/mole  for  the  temperature  dependence,  while  the  fuel  concentration  dependence, 
represented  as  [XJ®,  was  given  by  n  =  0.3.  Studies  of  soot  aggregates  found  in  these  flames  yielded  higher 
values  of  volume-mean  diameter,  a  larger  surface  area  per  unit  volume,  and  lower  values  of  the  aggregate 
number  concentration  as  compared  to  spherical  particle  assumptions.  Operating  pressure  studies  indicated 
that  a  power  law  representation  of  the  form  P",  represented  soot  volume  fraction  dependence  on  the  pressure, 
P,  and  confirmed  the  strong  sensitivity  of  soot  formation  to  pressure.  Finally,  high  pressure  diffusion  flame 
studies  revealed  the  onset  of  buoyant  instabilities  induced  by  changes  in  the  pressure.  A  joint  series  of 
experiments  and  computations  provided  strong  evidence  to  support  that  these  studies  isolate  Richardson 
number  as  the  only  variable  parameter.  _ 
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SUMMARY 


A  fundamental  study  of  soot  particle  inception  and  growth  processes  has  been  completed  which 
specifically  considered  the  effects  of  fuel  molecular  structure,  fuel  concentration,  temperature  and 
operating  pressure  on  the  formation  of  soot  particles.  Coannular  laminar  diffusion  flames  burning  in 
air  were  studied  for  pure  fuels  (methane  and  ethene)  as  well  as  fuel  mixtures  involving  ethene 
propene,  butane,  1-butene  and  1,3  butadiene  individually  mixed  with  methane  or  ethene.  The  fuel 
structure  studies  indicated  that  the  fuel  species  most  strongly  affected  the  particle  inception  process  for 
the  fuels  studied,  as  opposed  to  the  surface  growth  process.  Surface  growth  rate  constants,  when 
normalized  by  the  available  surface  area,  were  typically  within  a  factor  of  two  for  the  fuels  studied. 
These  findings  support  an  interpretation  in  which  inception  controls  the  maximum  amount  of  soot 
formed  in  the  flame.  Studies  in  which  fuel  concentration  and  flame  temperature  were  varied  through 
dilution  of  fuel  with  an  inert  species  indicated  that  temperature  effects  dominate  changes  due  to 
variations  in  the  concentration  under  most  conditions.  Local  temperature  and  concentration 
measurements  were  shown  to  be  absolutely  necessary  to  properly  interpret  the  effects  of  inert  dilution. 
The  analysis  yielded  an  apparent  activation  energy  of  94.5  kcal/mole  for  the  temperature  dependence. 
The  fuel  concentration  dependence,  represented  as  [XJ",  where  Xq  is  the  fuel  concentration  at  the 
exit  of  the  fuel  tube,  was  best  represented  by  n  =  0.3.  Studies  of  the  evolution  of  the  soot  particles 
in  laminar  diffusion  flames  indicated  that  clusters  or  aggregates  are  found  in  the  higher  parts  of  the 
flame.  These  clusters,  which  are  composed  of  nearly  spherical  primary  particles,  possess  different 
optical  properties  from  their  spherical  equivalents.  As  compared  to  the  Rayleigh  data  reduction  for 
spherical  particles,  the  aggregate  analysis  yielded  higher  values  of  volume-mean  diameter,  a  larger 
surface  area  per  unit  volume,  and  lower  values  of  the  aggregate  number  concentration.  The  present 
results  indicate  that  a  self-consistent  interpretation  of  the  light  scattering  properties  of  the  soot  aerosol 
is  afforded  by  recognizing  its  aggregate  structure.  Operating  pressure  effects  on  soot  formation  were 
also  studied  and  indicated  that  a  power  law  representation  of  the  form  of  P°,  where  P  is  the  operating 
pressure,  accurately  represented  the  behavior  of  the  soot  volume  fraction  as  the  pressure  was  varied. 
For  studies  involving  the  alkene  species,  ethene  and  propene,  the  value  of  n  was  found  to  be  close  to 
first  order  (l.OS  ±  0.06  and  1.16  ±  0.01,  respectively).  Studies  of  an  alkane  fuel,  ethane,  resulted 
in  a  significantly  larger  value  for  the  power  dependence  closer  to  second  order  in  n  (1.62  -  1.88, 
depending  on  the  flow  rate).  These  results  confirm  the  strong  dependence  of  soot  formation  on  the 
operating  pressure  and  furthermore  point  to  potential  fuel  structure  effects.  Finally,  studies  of  these 
high  pressure  diffusion  flames  have  resulted  in  the  observation  of  the  onset  of  buoyant  instabilities 
induced  by  changes  in  the  pressure.  A  simple  procedure  has  been  described  for  isolating  these 
buoyancy  effects  in  jet  diffusion  flame  experiments.  All  that  is  necessary  is  that  background  pressure 
be  varied  while  maintaining  constant  mass  flows  of  fuel  and  oxidizer  into  the  burner.  A  theoretical 
model  for  these  flames  indicates  that  this  procedure  isolates  Richardson  number  (or  relative  buoyancy 
force)  as  the  only  variable  parameter,  which  is  equivalent  to  varying  the  gravitational  acceleration.  A 
joint  series  of  experiments  and  computations  involving  a  pressurized  low  speed  methane/air  diffusion 
flame  has  provided  strong  evidence  to  support  the  theory. 
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1.0  RESEARCH  OBJECTIVES 

The  primary  objective  of  this  research  program  was  to  provide  additional  fundamental 
understanding  of  the  processes  which  control  soot  formation  under  conditions  applicable  to  future  gas 
turbine  engine  operation.  The  emphasis  of  the  studies  was  to  examine  the  physical  and  chemical 
mechanisms  which  control  soot  particle  inception  and  surface  growth.  The  inception  region  is  of 
particular  interest  because  it  has  been  argued  to  control  the  amount  of  soot  formed  during  the 
combustion  process.  In  the  present  studies,  detailed  measurements  of  both  the  inception  and  surface 
growth  phases  were  undertaken  to  provide  the  necessary  basis  for  evaluating  the  relative  importance 
of  these  fundamental  mechanisms  on  the  formation  of  soot  particles. 

One  particular  challenge  for  the  present  studies  was  to  acquire  data  over  a  wide  range  of  flame 
conditions.  In  this  respect,  fuel  species,  concentration,  temperature  and  operating  pressures  were 
selected  as  parameters  for  study.  The  results  and  conclusions  regarding  the  effects  of  each  of  these 
parameters  will  be  presented  in  the  following  sections.  These  parameters  were  selected  because  they 
are  viewed  as  being  critical  to  controlling  soot  formation  in  future  gas  turbine  engines  which  will  bum 
broader  specification  fuels  at  higher  pressures  than  currently  employed  in  such  engines.  Through  the 
understanding  gained  in  these  studies,  improved  prediction  capabilities  are  resulting  which  will  impact 
the  design  methodology  for  future  gas  turbine  combustors. 

2.0  RESEARCH  APPROACH 

The  basic  nature  of  the  objective  of  the  present  research  program  required  that  fundamental 
studies  be  conducted  under  well  controlled  conditions.  Additionally,  a  capability  to  provide  extensive 
measurements  of  the  soot  particle,  velocity,  temperature  and  gas  concentration  fields  was  also 
desirable.  For  these  reasons,  the  studies  were  conducted  in  a  laminar  diffusion  flame  environment 
burning  in  air  to  which  a  variety  of  diagnostic  techniques  were  applied. 
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Fuel  mixture  studies  employing  a  fuel  addition  approach  have  proven  useful  for  the  study  of  fuel 
molecular  structure  and  pressure  effects.  In  these  studies,  well  characterized  methane  and  ethene 
baseline  diffusion  flames  have  been  used.  To  these  baseline  flames,  controlled  amounts  of  an 
additional  fuel  species  are  added  under  conditions  in  which  the  total  carbon  atom  flow  rate  is  held 
constant.  In  this  manner,  the  effects  of  a  particular  fuel  molecular  structure  can  be  investigated  under 
similar  flame  conditions.  These  studies  were  carried  out  in  a  coannular  diffusion  flame  configuration 
under  both  atmospheric  and  elevated  pressure  conditions.  A  high  pressure  diffusion  flame  facility 
capable  of  operating  at  pressures  up  to  twenty  atmospheres  was  assembled  for  these  high  pressure 
studies,  and  investigations  have  been  conducted  for  a  pressure  range  of  one  to  ten  atmospheres. 

Laser-based  diagnostics  have  been  emphasized  for  the  soot  particle  and  velocity  field 
measurements  in  these  studies.  An  extensive  laser  light  scattering/extinction  apparatus  was  utilized 
for  the  soot  particle  measurements.  This  apparatus  has  the  capability  to  provide  information  on  the 
volume  fraction,  particle  diameter  and  number  concentration  of  soot  particles  present  in  the  flame. 
Using  light  extinction  and  scattering  measurements,  additional  extensive  scattering  measurements 
obtained  as  a  function  of  scattering  angle  and  laser  polarization  provided  a  redundant  measurement 
capability  to  the  extinction/scattering  ratio  technique.  These  capabilities  were  useful  in  studies 
attempting  to  account  for  the  aggregate  nature  of  the  mature  soot  particles.  Velocity  measurements 
using  a  laser  velocimetry  technique  were  obtained  to  determine  the  particle  paths  and  temporal  history 
of  the  soot  formation  process.  From  the  combined  laser  velocimetry  and  light  scattering 
measurements,  information  on  soot  particle  inception,  surface  growth  and  oxidation  were  obtained. 

These  particle  and  velocity  measurements  were  complemented  by  thermocouple  measurements  for 
temperature.  Additionally,  a  quadrapole  mass  spectrometer,  wavelength  selective  optical  absorption 
system  and  laser-induced  fluorescence  apparatus  have  been  employed  to  provide  information  on  gas 
phase  species  present  in  these  flames.  A  related  activity  involved  the  development  of  a  sonic  quartz 
microprobe  for  sampling  in  particle  laden  flows.  This  technique  provides  new  capabilities  needed  to 
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measure  gas  phase  species  in  both  the  precursor  and  surface  growth  regions.  In  the  particle 
diagnostic  area,  scanning  and  transmission  electron  microscopy  have  been  developed  to  allow 
complementary  measurements  of  primary  particle  size  in  soot  aggregates.  Thus,  an  approach, 
involving  carefully  controlled  experiments  to  which  extensive  diagnostic  measurement  techniques  were 
applied,  was  used  to  investigate  basic  mechanisms  important  in  soot  formation.  Through  this  process, 
effects  controlling  the  production  and  destruction  of  soot  particles  in  combustion  systems  could  be 
examined  in  a  quantitative  manner. 

3.0  RESEARCH  ACCOMPLISHMENTS 

3.1  Introduction 

In  the  following  sections,  the  accomplishments  and  results  obtained  during  the  current  research 
program  are  reviewed.  These  include  results  of  the  effects  of  fuel  molecular  structure,  concentration 
and  temperature  variations,  soot  aggregate  processes  and  operating  pressure.  An  additional  section 
dealing  with  buoyancy  effects  in  jet  diffusion  flames,  which  arose  out  of  the  high  pressure  studies, 
will  also  be  reviewed.  For  those  results  for  which  journal  articles  have  been  prepared,  a  summary  of 
the  results  will  be  presented  with  the  complete  article  included  as  an  appendix.  For  the  readers’ 
convenience,  a  section  describing  the  diagnostics  employed  in  these  studies  will  be  presented  first. 

This  section  will  serve  to  provide  the  details  of  the  individual  diagnostic  approaches  for  those 
interested  while  allowing  others  to  proceed  directly  to  the  results  of  interest. 

3.2  Experimental  Apparatus  and  Techniques 
3.2.1.  Burner  and  Flow  Apparatus 

Two  diffusion  flame  facilities  were  used  in  the  present  studies.  The  first  involves  an  atmospheric 
pressure  diffusion  flame  facility  and  the  second  is  a  high  pressure  diffusion  flame  apparatus. 
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3.2.I.I.  Atmospheric  Diffusion  Flame  Facility 

The  atmospheric  diffusion  flame  facility  consists  of  coannular  diffusion  flame  burner,  burner 
chimney,  positioning  system  and  gas  metering  system.  The  burner  has  a  coannular  configuration 
consisting  of  a  1.1  cm  fuel  tube  surrounded  by  a  10  cm  air  annulus.  The  air  passage  is  partially  filled 
with  glass  beads  followed  by  a  series  of  fine  screens  to  provide  flow  conditioning.  A  ceramic 
honeycomb  2.54  cm  in  thickness  is  used  at  the  exit  to  provide  a  uniform  flow  field.  The  fuel  tube 
which  extends  4.8  mm  above  the  ceramic  honeycomb  also  is  partly  filled  with  glass  beads  to 
condition  the  flow.  The  fuel  flow  can  consist  of  up  to  three  gases,  each  metered  with  a  separate 
rotameter.  This  allows  for  mixtures  of  fuels  as  well  a  nitrogen  dilution  of  the  fuel  for  temperature 
and  concentration  variations.  These  rotameters  have  been  calibrated  for  various  gases  using  a  soap 
bubble  meter  technique.  The  air  flow  is  metered  using  a  mass  flowmeter  which  can  monitor  flows  up 
to  5  SCFM  of  air.  To  protect  the  flame  from  room  disturbances,  a  metal  chimney  has  been 
incorporated  into  the  burner  facility.  The  chimney  translates  horizontally  with  the  burner  while 
sliding  vertically  within  the  chimney.  Slots  machined  in  the  chimney  provide  for  optical  access. 

The  burner  is  mounted  on  a  pair  of  motorized  translating  stages  which  provide  for  vertical  and 
horizontal  motion.  A  manual  translation  stage  is  also  included  to  allow  for  adjustment  in  the  second 
horizontal  direction  and  is  used  to  align  the  burner  with  the  laser  scattering  system.  The  motorized 
stages  are  used  to  traverse  the  burner  through  the  laser  beam  to  obtain  measurements  over  the  cross 
section  of  the  flame  at  a  particular  height  in  the  flame.  The  motorized  translation  stages  have  a 
position  resolution  of  0.0127  mm  which  is  sufficient  for  the  present  experiments.  Both  motorized 
stages  are  interfaced  to  and  controlled  by  an  IBM-XT  computer  using  the  general  purpose  interface 
bus  (GPIB-IEEE  488).  The  software  to  control  the  translation  stages  is  incorporated  into  the  data 
acquisition  program  for  the  laser  scattering  measurements. 
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3.2.1. 2.  High  Pressure  Diffusion  Flame  Facility 

In  order  to  provide  for  studies  at  elevated  pressure,  a  high  pressure  diffusion  flaj  "^e  facility  was 
constructed.  This  facility  is  composed  of  a  coannular  burner,  pressure  vessel,  positioning  system  and 
gas  metering  apparatus.  The  design  chosen  provides  for  a  significant  degree  of  similarity  between  the 
atmospheric  and  high  pressure  diffusion  flame  burners.  A  burner  identical  to  that  previously 
described  above  was  constructed.  This  burner  is  mounted  inside  a  high  pressure  vessel  capable  of 
withstanding  the  required  operating  pressures.  The  high  pressure  diffusion  flame  burner  is  shown  in 
figure  1. 

Because  of  the  mass  of  the  pressure  vessel,  the  approach  chosen  for  the  atmospheric  system  of 
vertically  translating  the  entire  burner  and  pressure  vessel  is  not  feasible.  Rather,  the  burner  is 
mounted  on  a  motorized  translation  stage  located  internal  to  the  pressure  vessel.  Electrical 
connections  are  made  through  the  base  of  the  pressure  vessel  allowing  external  control  of  the  vertical 
position  of  the  burner.  To  provide  horizontal  movement,  the  base  of  the  pressure  vessel  is  attached  to 
a  precision  ball  bearing  stage  which  requires  only  a  few  pounds  of  force  to  move  horizontally  with 
weights  as  large  as  1500  pounds  attached.  Thus,  a  relatively  small  motorized  translation  stage  can  be 
used  to  provide  horizontal  movement.  As  with  the  atmospheric  burner,  these  motorized  stages  are 
capable  of  computer  control  using  the  laboratory  IBM-XT  personal  computer.  The  only  disadvantage 
with  its  approach  is  that  sufficient  space  must  be  provided  within  the  pressure  vessel  to  house  the 
vertical  motion  assembly  and  the  flexible  gas  supply  lines  for  the  burner. 

To  accommodate  this  limitation,  the  pressure  vessel  has  been  constructed  as  two  separate 
sections.  The  lower  section  is  constructed  of  8  inch  diameter  carbon  steel  seamless  pipe  which  is  24 
inches  in  length.  The  vertical  translating  mechanism  is  mounted  to  the  base  of  this  lower  section.  In 
addition  to  the  electrical  connection  for  the  motorized  translation  stage,  provisions  for  connection  of 
fuel  and  air  lines  are  also  made  through  this  base.  Internal  connections  to  the  diffusion  flame  burner 
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Figure  1 .  High  pressure  diffusion  flame  apparatus 


can  be  made  with  low  pressure  tygon  tubing  since  no  pressure  difference  exists  across  the  tubing  wall 
inside  the  pressure  vessel. 

The  upper  section  of  the  pressure  vessel,  which  provides  for  optical  access  to  the  flame  region,  is 
30  inches  in  length  and  constructed  from  6  inch  diameter  carbon  steel  pipe.  This  section  includes 
four  windows  which  are  two  inches  in  diameter  and  1/2  inch  thick.  These  windows  allow  sufficient 
optical  access  for  both  the  laser  light  scattering  and  the  laser  velocimetry  measurements.  The 
operating  pressure  in  the  burner  is  adjusted  using  manual  valves  located  in  the  exhaust  line  of  the 
burner.  Connection  to  this  exhaust  system  is  made  through  a  second  flange  to  which  a  2  inch 
diameter  exhaust  line  has  been  connected.  At  this  time,  no  cooling  is  required  for  the  burner  or 
exhaust  section.  Temperature  measurements  indicated  that  for  conditions  typical  for  the  present 
experiments,  the  burner  does  not  heat  up  significantly. 

Control  and  measurement  of  the  fuel  and  air  flow  rates  are  accomplished  using  mass  flow  meters 
and  controllers.  These  meters  are  insensitive  to  the  operating  pressure  and  thus  can  be  calibrated  at 
atmospheric  pressure  while  providing  accurate  metering  at  elevated  pressure.  The  fuel  metering 
system  allows  up  to  three  gases  to  be  mixed  using  the  independent  mass  flow  controllers.  These 
controllers  maintain  a  constant  mass  flow  independent  of  the  pressure  drop  across  the  meter.  This 
assures  a  constant  mass  flow  rate  of  fuel  as  well  as  a  convenient  start-up  capability.  Otherwise,  the 
fuel  flow  rate  would  decrease  as  the  operating  pressure  was  increased,  requiring  continual 
readjustment.  Because  of  the  relatively  large  air  flow  rate  required,  a  mass  controller  approach  is  not 
feasible  for  the  air  flow  metering.  The  air  flow  rate  is  monitored  by  a  mass  flowmeter  which  requires 
adjustment  as  the  operating  pressure  is  varied.  This  present  burner  has  been  designed  to  operate  at 
pressures  as  high  as  twenty  atmospheres.  Stable  flames  at  pressures  up  to  ten  atmospheres  have  been 
studied  to  date. 
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3.2.2.  Non-Intrusive  Optical  Diagnostics 

A  more  complete  understanding  of  soot  formation  and  destruction  processes  requires  the 
measurement  of  the  key  properties  characterizing  the  combusting  flow.  These  property  measurements 
include  soot  particle  size  and  concentration,  temperature,  gas  phase  species  concentration  and 
velocity.  Once  determined,  this  set  of  measurements  provides  the  basis  from  which  to  develop  an 
understanding  of  the  temporal  evolution  of  soot  particle  formation  in  combustion  environments.  To 
address  some  of  these  needs,  non-intrusive  optical  techniques  have  been  employed.  Laser  light 
scattering  has  been  used  to  characterize  the  soot  particle  field  while  laser  velocimetry  has  been 
utilized  for  velocity  measurements.  An  optical  absorption  apparatus  has  been  assembled  to  allow 
wavelength  selective  absorption  measurements  to  be  obtained.  These  techniques  were  complemented 
by  intrusive  probe  sampling  techniques  to  be  described  in  a  later  section. 

3. 2 .2.1.  Laser  Light  Scattering  Apparatus 

The  laser  light  scattering  apparatus,  utilizing  a  4W  argon  ion  laser  as  the  light  source,  provided 
for  particle  extinction  and  scattering  measurements.  Scattering  measurements  were  made  primarily  at 
90®  although  the  system  can  be  used  to  obtain  measurements  at  45®  and  135®  (see  figure  2).  The 
laser  source  was  modulated  using  a  mechanical  chopper  operating  at  1  kHz  to  allow  for  synchronous 
detection  of  the  transmitted  and  scattered  light  signals.  A  polarization  rotator  was  also  incorporated 
in  the  system  to  allow  adjustment  of  the  polarization  of  the  incident  light  beam.  The  laser  beam  was 
focused  in  the  burner  using  a  30  cm  focal  length  lens  which  results  in  a  probe  beam  diameter  of  0.02 
cm.  Typically  the  laser  was  operated  at  the  514.5  nm  laser  line  with  an  output  power  of  about  0.5W. 

The  transmitted  light  signal  was  detected  using  a  silicon  photodiode.  The  laser  light  intensity  was 
reduced  by  a  neutral  density  filter  (N.D.  2.0)  to  a  level  suitable  for  linear  photodiode  response.  The 
scattered  light  was  detected  using  a  photomultiplier  detector  (PMT).  The  PMT  had  a  narrowband 
filter  center  at  514.5  nm  with  a  1  nm  bandwidth  incorporated  in  the  PMT  housing  to  help  reject  light 
other  than  that  generated  from  the  particle  scattering  event.  A  pinhole  with  a  diameter  of  1  mm, 
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Figure  2.  Schematic  of  laser  scattering/extinction  system: 

C  =  chopper,  M  =  mirror,  PR  =  polarization  rotator,  L  =lens, 

S  =  spatial  filter  (circular  aperture),  NDF  =  neutral  density 
filter,  PD  =  photodiode,  P  =  polarization,  F  =  narrowband  filter, 
PMT  =  photomultiplier 


located  in  front  on  this  filter,  limited  the  measurement  volume  length  along  the  beam  to  0. 1  cm.  The 
IS  cm  focal  length  collection  lens  used  to  focus  the  scattered  light  onto  the  PMT  was  arranged  to 
provide  unity  magnification.  Thus,  the  collection  volume  was  a  cylinder  approximately  0.02  cm  in 
diameter  and  .1  cm  in  length  resulting  in  a  volume  of  3.1  x  10  cm^.  The  collection  lens  was 
preceded  by  a  polarization  filter  to  allow  polarization  discrimination  for  the  scattered  light  detected. 
The  collection  solid  angle  was  determined  by  a  1.27  cm  aperture  which  is  located  between  the 
polarization  filter  and  the  collection  lens.  This  aperture  limited  the  collection  angle  to  approximately 
2“. 

The  output  from  each  of  the  detectors  was  input  into  separate  lock-in  amplifiers  which  were 
interfaced  to  the  IBM-XT  computer  over  the  GPIB  bus.  These  units  are  capable  of  full  computer 
control  which  allows  for  autoranging  of  the  lock-in  amplifiers.  Computer  software  to  control  the 
lock-in  amplifiers  has  been  integrated  with  the  stepper  motor  control  programs  to  provide  a  complete 
data  acquisition  routine.  These  procedures  provide  a  high  degree  of  user  independent  data 
acquisition. 

Data  taken  from  the  lock-in  amplifiers  was  stored  on  the  IBM-XT  internal  disk.  The  light 
scattering  measurements  required  extensive  data  reaction  to  yield  particle  size  information.  Such 
calculations  were  performed  on  computers  more  capable  than  the  laboratory  computer  used  for  the 
data  acquisition.  Thus,  the  IBM-XT  was  linked  over  the  Mechanical  Engineering  Department’s 
network  to  the  University’s  VAX  computer  facilities.  This  allowed  for  rapid  data  analysis  and 
graphical  output.  The  necessary  data  reduction  and  plotting  software  have  been  developed  for  this 
facility. 

The  scattered  light  detection  system  was  calibrated  to  account  for  effects  of  the  incident  laser 
power,  sample  volume,  light  collection  efficiency,  photomultiplier  sensitivity  and  electronic  gain  of 
the  system.  The  calibration  was  accomplished  by  passing  ethene,  a  gas  with  a  known  Rayleigh 
scattering  cross  section,  through  the  fuel  passage  of  the  burner  and  measuring  the  resulting  scattered 
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light.  This  procedure  allowed  for  an  absolute  determination  of  the  differential  scattering  cross  section 
per  unit  volume,  Q(0),  which  is  the  power  scattered  in  the  direction  6  per  unit  incident  flux. 

Several  approaches  can  be  taken  to  determine  the  particle  size  from  the  above  defined  cross 
sections.  Three  are  of  interest  with  respect  to  the  instrumental  set  up  discussed  previously.  These  are 
the  ratio  of  scattering  to  extinction,  (90°)/Kgjj,  the  ratios  of  the  scattering  signal  at  two  angles 
(dissynunetry  ratio)  Qw(45‘’)/Qw(90°)  and  Qw(45‘’)/Qyy(135®),  and  the  polarization  ratio, 
Q|jj(90°)/Qyy(90®);  the  scattering  measurement  angle  is  specified  for  each  quantity.  Each  of  these 
ratios  has  a  strong  dependence  on  particle  size  and  represents  a  redundant  set  of  data  in  terms  of  the 
particle  properties  for  particles  in  the  Mie  size  region  (D  >  X).  For  particles  in  the  Rayleigh  size 
region  (D  <  <  X),  scattering  is  isotropic  so  the  dissymmetry  ratio  has  a  value  of  unity  and  the 
defined  polarization  ration  approaches  zero.  Thus,  in  this  small  particle  limit,  only  the 
scattering/extinction  ratio  yields  size  information.  However,  for  soot  formation  processes,  the  rapid 
coagulation  and  surface  growth  processes  lead  quickly  to  particle  sizes  characteristic  of  the  Mie  size 
region.  Because  of  its  applicability  throughout  the  Rayleigh  and  Mie  size  regimes, 
scattering/extinction  measurements  have  proved  most  useful  in  prior  studies  in  this  laboratory. 
Recently,  efforts  involving  dissymmetry  and  polarization  ratio  have  been  undertaken  to  consider  the 
effects  of  particle  agglomeration  on  light  scattering  measurement  [1]. 

The  laser  light  scattering  apparatus  described  ^^ove  provides  a  versatile  instrument  for 
characterizing  the  volume  fraction,  particle  size  and  number  concentration  of  soot  particles  formed  in 
flames.  Thus,  processes  starting  from  particle  inception  through  surface  growth  and  eventually 
oxidation  can  be  measured.  Application  of  these  techniques  in  the  present  studies  has  allowed  an 
extensive  examination  of  soot  formation  in  diffusion  flames  [2-8]. 
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3.2.2^  Laser  Velocimetry  Apparatus 

Laser  velocimetry  (LV)  has  been  used  to  measure  the  velocity  field  to  investigate  the  transport 
processes  in  the  laminar  diffusion  flame.  The  laser  velocimeter  is  a  conventional  dual  beam  system 
utilizing  an  argon  ion  laser  operated  at  the  514. S  nm  laser  line  to  provide  the  incident  laser  beam. 
This  laser  beam  is  then  split  into  two  parallel  beams  of  equal  intensity  which  are  separated  by  a  SO 
mm  spacing.  These  beams  are  then  focused  to  form  the  LV  probe  volume  using  a  lens.  The 
resulting  ellipsoidal  probe  volume  will  vary  in  size  as  a  function  of  the  focal  length  of  the  focusing 
lens.  Typically,  a  lens  with  a  250  mm  focal  length  is  used  resulting  in  a  calculated  probe  diameter  of 
0.16  mm  and  a  length  of  1.7  mm.  Scattered  light  is  then  collected  with  a  second  lens  and  focused 
onto  a  photomultiplier  tube. 

Velocity  measurements  are  obtained  from  the  modulation  frequency  of  the  intensity  of  light 
scattered  from  particles  traversing  the  probe  volume  formed  by  the  two  laser  beams.  This  scattered 
light  can  be  viewed  as  a  signature  of  the  interference  pattern  resulting  from  the  constructive  and 
destructive  interference  of  the  laser  beams.  The  period  between  the  intensity  maxima  in  the  scattered 
light  signal  is  thus  proportional  to  the  velocity  of  the  particles  in  the  probe  volume.  Since  the  signals 
occur  randomly,  digital  counting  techniques  (burst  processing)  are  used  to  obtain  the  frequency  and 
thus  the  velocity.  The  digital  output  of  the  processor  is  stored  on  an  IBM-XT  computer. 

In  order  to  obtain  measurements  in  situations  in  which  the  direction  of  the  flow  reverses,  a  Bragg 
cell  is  placed  between  the  beam  splitter  and  the  focusing  lens.  This  cell  imparts  a  frequency  offset  to 
one  of  the  two  beams  used  in  forming  the  LV  probe  volume.  This  frequency  shift  allows 
determination  of  the  direction  as  well  as  the  magnitude  of  the  velocity  component  by  providing  a 
frequency  offset  for  zero  velocity.  The  introduction  of  this  frequency  offset  is  also  useful  for 
measurements  in  low  velocity  flows. 

Measurements  of  the  axial  and  radial  components  of  the  velocity  field  have  been  used  previously 
to  construct  the  streak  lines  which  soot  particles  follow  in  the  laminar  diffusion  flame  studies  [3]. 
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This  information  has  been  critical  in  gaining  information  on  the  temporal  evolution  of  the  soot  particle 
field. 

3^23  Optical  Absorption  System 

An  optical  absorption  apparatus  was  assembled  to  provide  for  wavelength  selective  measurements 
of  absorbing  species  in  the  flame  studies.  This  system  consisted  of  a  tungsten  filament  lamp  driven 
with  a  highly  stable  D.C.  power  supply.  The  light  emitted  from  the  lamp  was  focused  on  the  axis  of 
the  burner  and  then  passed  into  a  spectrometer  for  detection.  Both  photomultiplier  and  photodiode 
detectors  were  used  for  signal  detection  with  the  choice  of  detector  subject  to  the  wavelength  range  of 
interest.  The  spatial  resolution  of  measurements  was  adjusted  by  the  use  of  suitable  ^ertures  placed 
before  and  after  the  flame. 

In  order  to  provide  for  rejection  of  background  radiation  and  noise,  the  incident  light  was 
modulated  with  a  mechanical  chopper.  Signals  were  then  measured  using  a  lock-in  detection 
technique.  For  cases  where  extremely  small  signals  are  encountered,  a  difrerential  reference  channel 
technique  was  assembled  as  well.  In  this  approach,  a  portion  of  the  incident  light  was  directed 
around  the  flame  region.  This  beam  was  modulated  using  a  separated  mechanical  chopper  operating 
at  a  different  frequency  from  that  used  for  the  beam  passing  through  the  flame.  This  beam  was  then 
recombined  with  the  main  beam  and  input  into  the  spectrometer.  Two  lock-in  amplifiers,  each 
operating  at  the  separate  frequencies  of  the  mechanical  choppers,  were  used  to  measure  the  absorption 
signal  as  well  as  the  incident  light  signal.  This  approach  eliminates  signal  variations  resulting  from 
changes  in  the  lamp  source  and  provided  sensitivity  of  three  parts  in  a  thousand.  Wavelength 
dqiendent  measurements  of  the  absorption  process  were  then  obtained  by  scanning  the  spectrometer 
over  a  desired  wavelength  region  and  recording  the  resulting  signal.  Data  acquisition  and  control  of 
the  spectrometer  were  provided  by  an  IBM  XT  laboratory  computer. 

Analysis  of  the  measurements  was  intended  to  provide  information  on  both  gas  phase  and  soot 
particle  concentration.  The  measurement  of  gas  phase  species  can  be  based  on  the  observation  that. 
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as  molecular  weight  increases,  absorption  shifts  to  long  wavelengths.  Thus,  by  a  systematic 
evaluation  of  the  wavelength  dependence  of  the  absorption  signals  as  a  function  of  location  in  the 
flame,  the  evolution  of  large  molecular  weight  species  can  be  followed. 

3^3  Intrusive  Probe  Sampling  for  Species  Concentration  Measurements 
3.2.3. 1  Ouadranole  Mass  Spectrometer  and  Conventional  Quartz  Microprobe 

In  order  to  provide  more  quantitative  information  on  the  species  concentration,  a  mass 
spectrometer  system  was  assembled  for  these  studies.  This  facility  was  obtained  under  support 
through  a  DoD  University  Research  Instrumentation  Award  by  AFOSR. 

Gas  samples  from  the  flames  under  study  were  sampled  using  a  conventional  quartz  microprobe 
technique  similar  to  the  approach  taken  by  Smyth  et  al.  [9].  The  mass  spectrometer  selected  for  these 
studies  utilizes  a  quadrapole  mass  filter  with  a  mass  range  of  1  to  SOO  amu.  This  unit  (an  Extrel 
Model  C-SO)  has  excellent  mass  resolution  (0. 1  amu)  and  high  sensitivity.  The  probe  and  mass 
spectrom^er  assembly  were  differentially  pumped.  A  pressure  of  between  1  and  10  torrs  using  a 
vacuum  roughing  pump  was  achieved  directly  behind  the  quartz  orifice  which  is  nominally  1(X)  urn  in 
diameter.  Sampling  from  this  low  pressure  region  was  provided  by  a  second  probe  of  identical 
construction.  This  second  probe  was  used  to  introduce  gases  into  the  mass  spectrometer  which  was 
maintained  at  a  pressure  near  10*^  torr  using  a  turbomolecular  pumping  system. 

Radial  profiles  of  the  species  concentration  as  a  function  of  axial  position  were  obtained  by 
transversing  the  flame  in  the  radial  direction.  Similar  to  the  optical  diagnostic  ^proaches,  data 
acquisition  and  the  control  of  the  mass  spectrometer  have  been  implemented  using  a  laboratory 
personnel  computer  system.  This  capability  allowed  for  rapid  and  convenient  data  acquisition,  while 
also  affording  opportunities  for  other  probing  approaches,  such  as  rapid  insertion  sampling.  For 
situations  where  more  species  specific  information  is  required,  gas  chromatography  techniques  were 
also  available  within  the  laboratory.  These  techniques,  which  were  of  limited  use  in  terms  of  the 
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flame  sampling,  do  provide  capabilities  useful  in  unambiguously  identifying  the  specific  nature  of  an 
observed  mass  peak,  as  well  as  for  comparing  calibration  of  particular  species. 

This  system  was  used  extensively  to  obtain  concentration  information  in  diluted  methane  flames. 
These  results  will  be  described  below  in  detail.  Because  the  quartz  microprobe  could  not  be 
employed  in  particle  laden  regions  of  the  flame  due  to  probe  clogging,  a  separate  sonic  probe  for 
sampling  in  these  regions  was  developed  and  is  described  next. 

3. 2. 3. 2  Electromechanical  Sonic  Probe 

As  mentioned  above,  the  need  to  determine  gas  species  concentrations  in  regions  containing  soot 
particles  presents  serious  diagnostic  challenges.  In  general,  optical  diagnostic  techniques  can  measure 
only  a  limited  number  of  species  and  in  some  cases  can  not  be  reliably  applied  in  the  presence  of  soot 
particles.  Sampling  probe  techniques  usually  are  limited  by  orifice  clogging  problems  which  can  only 
be  overcome  by  increasing  the  orifice  dimensions  at  the  sacrifice  of  spatial  resolution  of  the 
measurements.  To  overcome  the  aforementioned  orifice  clogging  problem,  a  novel  sonic  sampling 
probe  has  been  developed  and  applied  to  a  series  of  laminar  diffusion  flames  containing  various 
amounts  of  soot.  A  patent  disclosure  for  this  probe  was  submitted;  and  subsequently,  patent  rights 
have  been  transferred  to  the  U.S.  government.  The  details  of  the  design  and  operation  of  this  probe 
were  presented  at  the  Eastern  Section  Meeting  of  The  Combustion  Institute  [10].  A  summary  of  the 
operation  of  the  sonic  probe  along  with  the  results  of  a  set  of  comparisons  obtained  to  validate  its 
operation  is  given  below. 

The  approach  employed  to  prevent  clogging  of  the  probe  orifice  involves  mechanically  oscillating 
a  wire  through  the  orifice  region  using  a  spring  loaded  solenoid  plunger  whose  driving  circuit  is 
electronically  interrupted  periodically  (see  figure  3).  The  constant  motion  of  the  wire  relative  to  the 
quartz  tube  is  responsible  for  keeping  the  orifice  open.  For  the  present  probe,  the  effective  orifice  is 
the  annulus  formed  by  an  oscillating  380  nm  diameter  Nichrome  wire  of  uniform  cross  section  and 
the  400  fim  diameter  orifice  in  the  quartz  tube.  Therefore,  the  effective  orifice  is  equivalent  to  a  125 
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fim  circular  hole.  The  wire  is  soldered  into  a  1.59  mm  (1/16  inch)  tube  which  is  attached  to  a  spring 
loaded  soft  iron  solenoid  plunger.  The  solenoid  is  activated  with  a  24  volt  power  supply  providing 
0.3  amperes  current.  The  solenoid  circuit  is  interrupted  with  the  aid  of  a  mini  relay  which  is 
activated  by  a  regular  series  of  5  volt  pulses  from  a  function  generator.  When  the  solenoid  circuit  is 
complete,  the  solenoid  retracts  the  soft  iron  plunger  and  consequently  the  wire  is  partially  withdrawn. 
A  small  portion  of  the  wire  still  extends  beyond  the  orifice  because  a  constant  effective  orifice  size  is 
desired.  The  inward  motion  of  the  plunger  compresses  the  spring  behind  the  plunger.  When  the  mini 
relay  interrupts  the  solenoid  circuit,  the  compressed  spring  expands  and  returns  the  plunger  and  the 
wire  to  its  original  position.  Thus,  the  wire  can  be  made  to  oscillate  through  the  orifice  region  at  the 
desired  frequency.  For  the  present  system,  an  oscillation  frequency  of  30  Hz  has  been  selected. 

The  operation  of  the  probe  was  tested  by  sampling  from  a  flame  burning  a  mixture  of  methane 
(5.6  cm^/s)  and  1-butene  (1.05  cm^/s)  in  air  (2.75  scfin)  where  the  local  soot  volume  fraction  was 
4.0  X  10*^.  Previous  studies  of  this  flame  using  a  probe  with  a  1  mm  diameter  orifice  resulted  in 
clogging  within  15  seconds  under  isokinetic  sampling  conditions.  During  initial  testing,  the  pressure 
downstream  of  the  probe  was  monitored  as  an  indicator  of  clogging  problems.  The  pressure  remained 
constant  at  5  torr  for  a  period  of  5  minutes  indicating  the  probe  did  not  clog  during  this  time  period. 
By  manually  blocking  the  probe  orifice,  it  was  established  that  should  the  probe  clog,  the  pressure 
inside  the  probe  would  fall  below  I  torr.  To  further  validate  operation  of  the  probe,  species 
measurements  with  gas  chromatography  were  carried  out  in  the  flame  described  above.  In  order  to 
use  this  low  pressure  sampling  probe  in  conjunction  with  gas  chromatography,  a  cylinder  piston 
device  (figure  4)  was  designed  and  fabricated  to  compress  the  sample  from  5  torr  to  15  psi.  The 
measurements  from  this  probe  were  compared  with  results  obtained  from  a  water  cooled  stainless  steel 
probe  having  a  2  mm  orifice  (which  clogs  within  5  minutes  under  isokinetic  sampling  conditions). 
Good  agreement  («  10%)  is  observed  between  the  two  probes  for  measurements  obtained  along  the 
centerline  from  the  upper  regions  of  this  flame.  Figure  5  presents  this  comparison  in  a  plot  of  CO 
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Axial  Position,  z  (cm) 


Figure  5.  Comparison  of  CO  measurements  from  the  Electromechanical  Sonic  Probe  with 
similar  measurements  from  a  water  cooled  stainless  steel  probe. 


Radial  Position,  mm 


Figure  6.  Radial  profiles  of  CO  mole  fraction  from  probes  with  three  different  orifice  sizes 
demonstrate  gradient  broadening.  The  symbols  identified  on  the  plot  refer  to  the  orifice 
diameters. 
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mole  fraction  versus  axial  position  in  the  flame.  Other  species  also  compare  as  well  (data  not 
shown).  This  probe  was  also  c'^mpared  with  a  conventional  sonic  probe  (orifice  85  /rm)  in  a 
relatively  soot  free  methane  flame  and  again  the  measurements  compared  well  for  the  upper  regions 
of  the  flame.  However,  some  differences  have  been  observed  in  the  lower  parts  of  the  flame  where 
steep  concentration  gradients  are  present. 

In  order  to  more  closely  examine  these  differences,  a  pure  methane  flame  (CH4:  5.7  cc/s.  Air: 
795  cc/s),  that  was  previously  studied  by  Mitchell  [11],  was  selected  .  Radial  particle  measurements 
of  the  CO  mole  fraction  at  a  height  of  12  mm  were  made  with  three  different  probes  -  two 
conventional  sonic  probes  and  the  electromechanical  sonic  probe  (EMS).  The  probe  orifice  sizes 
were  estimated  by  measuring  the  flow  rate  of  room  air  through  the  probes  under  choked  flow 
conditions.  The  EMS  probe  has  the  largest  orifice  with  an  equivalent  circular  diameter  of  0.184  mm. 
The  differences  between  various  probe  measurements  that  were  observed  in  the  lower  part  of  the 
flame  are  related  to  the  positioning  uncertainty  of  the  probe  and  differences  in  the  spatial  resolutions 
of  each  probe.  In  the  lower  region  of  diffusion  flames,  the  species  concentration  gradients  are  ste^ 
on  the  air  side  of  the  flame  front.  Therefore,  in  this  region  uncertainties  in  probe  positioning  and 
differences  in  spatial  resolutions  have  to  be  kept  in  mind  when  making  comparisons.  The  effect  of 
the  finite  probe  orifice  diameter  results  in  "broadening"  or  decreasing  the  steepness  of  the  species 
gradient  and  will  be  referred  to  as  "gradient  broadening". 

Figure  6  compares  the  radial  measurements  of  CO  mole  fraction  with  the  three  different  probes. 
Since  complete  radial  profiles  across  the  entire  flame  were  not  obtained,  some  uncertainty  in  the 
location  of  the  centerline  of  the  flame  exists.  In  ordei  to  examine  the  probe  orifice  effects  on  gradient 
broadening,  the  data  in  figure  6  are  shifted  slightly  to  provide  an  overlap  of  the  region  of  steepest 
gradient  for  each  data  set.  For  the  region  between  the  centerline  and  the  peak  CO  mole  fraction 
location,  there  is  no  significant  difference  between  the  measurements  from  the  three  probes.  On  the 
air  side  of  the  peak  CO  mole  fraction,  the  probe  with  the  smaller  orifice  gave  CO  values  displaying 
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steeper  radial  gradients.  This  is  understandable  since  a  smaller  orifice  results  in  better  spatial 
resolution  and  less  gradient  broadening. 

For  the  present  EMS  probe,  a  nichrome  wire  is  used  to  keep  the  orifice  of  the  probe  open. 
Attempts  at  using  a  quartz  fiber  have  not  been  successful,  since  this  fiber  becomes  brittle  once  it  is 
exposed  to  the  flame  and  breaks  easily.  Since  the  measurements  from  the  probe  with  the  nichrome 
wire  compare  well  with  measurements  using  the  other  probes,  this  probe  has  been  used  to  obtain 
radial  profiles  of  gas  phase  species  in  the  sooty  regions  of  the  flames  currently  under  study.  Some 
concerns  remain  about  potential  catalytic  effects  due  to  the  nichrome  wire  and  must  be  considered  as 
the  EMS  probe  is  used  in  other  environments. 

Table  1  indicates  the  fuel  and  air  flow  rates  for  the  profiles  studied  with  the  various  probes, 
while  figures  7  and  8  show  typical  radial  profile  measurements  of  N2,  CO2,  O2  and  CO  in  these 
flames.  The  symbols  represent  actual  measurements  and  the  lines  through  the  symbols  are  extended 
across  the  flame  centerline  by  symmetry.  Fuel  species  measurements  and  moisture  estimates  are  not 
plotted. 


Table  1:  Fuel  and  air  flow  rates  (cm^/s)  for  the  overventilated  flame  studies 


FLAME 

CH4 

C4H,o 

C4H8 

AIR 

Pure  Methane 

9.8 

- 

- 

1298 

Methane/Butane 

5.6 

1.05 

“ 

1298 

Methane/ 1 -Butene 

5.6 

- 

1.05 

1298 

Mitchell 

5.7 

795 

25 


0J10 


Figure  7.  Radial  profiles  of  Nj,  COj,  Oj  and  CO  in  the  pure  methane  flame  at  a  height  of  7.874 
cm. 


Figure  8.  Radial  profiles  of  N2,  COj,  O2  and  CO  in  the  pure  methane  flame  at  a  height  of  9.906 
cm. 
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3.2.4  Diagnostics  for  Temperature  Measurements 

Besides  the  particle,  velocity  and  species  concentration  measurements,  temperature  measurements 
were  required  to  complete  the  data  base  needed  to  analyze  the  soot  formation  and  destruction 
processes.  For  these  measurements,  an  intrusive  measurement  approach  using  fine  wire 
thermocouples  was  utilized. 

Temperature  profiles  were  obtained  using  Pt/Pt  -  10%  Rh  fine  wire  thermocouples.  A  rapid 
insertion  technique  similar  to  that  described  by  Kent  and  Wagner  [12]  was  employed  to  reduce  the 
effects  of  soot  particle  disposition.  A  motorized  translation  stage,  which  is  controlled  by  a  computer, 
allowed  rapid  positioning  of  the  thermocouple  in  the  flame.  After  inserting  the  thermocouple  and 
obtaining  the  reading,  the  thermocouple  was  retracted  to  a  position  on  the  lean  (air)  side  of  the  flame 
where  soot  particles  deposited  on  the  thermocouple  can  be  oxidized.  By  repeating  the  same  sequence 
for  each  measurement  position,  a  complete  radial  profile  of  the  temperature  field  was  obtained. 

Recent  measurements  have  shown  that  this  rapid  insertion  technique  significantly  reduces  the  effects 
of  soot  particle  deposition  and  improves  agreement  with  other  techniques  [3,12].  The  difficulties  of 
using  thermocouples,  particularly  in  sooting  flames,  are  well  documented.  Nonetheless,  these 
thermocouple  measurements  provided  a  suitable  basis  from  which  to  investigate  the  evolution  of  the 
temperature  field. 

3.2.5  Particle  Size  Analysis 

The  laser  light  scattering  and  extinction  apparatus  provides  the  particle  measurement  capability 
for  the  present  studies.  Analysis  of  the  scattering  measurements  is  based  on  a  MIE  theory  approach 
for  spherical  particles.  For  completeness,  a  brief  review  of  the  theoretical  basis  for  this  analysis  is 
presented. 

In  general,  the  interaction  of  a  particle  and  light  wave  is  dependent  on: 

1 .  index  of  refraction  of  the  particle,  m 

2.  particle  size,  D 
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3.  particle  number  density,  N 

4.  particle  size  distribution,  P(D) 

5.  particle  shape 

6.  wavelength  of  the  scattered  light 

7.  scattering  geometry 

Obviously,  measurement  of  any  one  of  the  above  requires  knowledge  of  all  of  the  other  variables.  In 
a  typical  combustion  environment  containing  soot  particles,  only  the  scattering  geometry  and  the 
wavelength  of  the  light  can  be  specified  beforehand.  The  remaining  unknowns  are  determined  by  an 
appropriate  combination  of  measurements  to  reduce  the  number  of  free  parameters  along  with  some 
reasonable  assumptions.  It  is  these  assumptions  which  can  result  in  inaccuracies  and  therefore  must 
be  considered  carefully. 

Combinations  of  scattering  measurements  utilizing  the  angular  and  polarization  dependencies  of 
the  scattering  process,  or  scattering  and  absorption  measurements,  have  been  used  to  determine 
particle  number  density  and  size  [2-4].  In  such  an  approach,  the  index  of  refraction  is  taken  to  be 
known,  the  particles  are  assumed  to  be  spherical  and  a  size  distribution  is  assumed  (e.g.,  log-normal). 

The  determination  of  particle  size  and  number  concentration  from  light  scattering  data  depends 
primarily  on  our  capability  to  relate  the  measurement  of  scattered  and  absorbed  light  intensity  to  the 
particle  properties.  Appropriate  theories  have  been  developed  for  spherical  particles  and  have  been 
widely  applied  for  aerosol  measurements.  The  particle  properties  of  interest  are  the  differential 
scattering  cross  section,  Cjj  (0),  for  production  of  scattered  light  at  a  specified  direction  and  total 
cross  section  for  a  specified  particle  size,  The  subscripts  i  and  j  for  the  scattering  cross  section 
assume  letters  v  or  h  according  to  whether  the  state  of  polarization  of  the  scattered  (i)  and  incident  (j) 
radiation  is  perpendicular  or  parallel,  respectively,  to  the  plane  of  observation.  For  spherical  particles 
of  isotropic  material,  only  the  case  of  i  =  j  is  of  interest  because  cross  polarization  effects  are  absent. 
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The  relationship  between  these  cross  sections  and  the  experimentally  measured  quantities  are 
expressed  as 

Is  =  loNCii  G  (1) 

and 


lo 


L 

exp-  J 
O 


N  Ce^^t  (^.m)  ds 


(2) 


where  Ig,  1  and  1^  are  the  scattered,  transmitted  and  incident  intensities  respectively,  N  is  number 
concentration  of  particles,  D  is  diameter,  m  is  the  refractive  index,  and  G  is  a  constant  which 
involves  factors  related  to  the  sample  volume  and  detection  instrumentation.  It  should  be  noted  that 
the  scattering  cross  section  measurements  given  by  equation  (1)  are  point  measurements,  while  the 
extinction  cross  section  is  an  integrated  quantity  over  the  path  length  L.  Thus,  to  obtain  local  particle 
extinction  values,  a  data  inversion  technique  must  be  utilized  [19], 

In  an  experiment,  the  quantities  which  are  determined,  once  appropriate  calibration  factors  have 
been  introduced,  are 

Qii  =  N  Cii  (3) 

and 

Kext  =  N  (4) 

where  Qy  is  the  volumetric  scattering  cross  section  and  is  the  extinction  coefficient. 

Several  approaches  can  be  taken  to  determine  the  particle  size  from  the  above  defined  cross 
sections.  Three  are  of  interest  with  respect  to  the  instrumental  set  up  discussed  previously.  These  are 
the  ratio  of  scattering  to  extinction,  Qw(90‘’))/Kgjj|,  the  ratios  of  the  scattering  signal  at  two  angles 
(dissymmetry  ratio)  Qw(45‘’)/Qw(90°)  and  Qw(45‘’)/Qw(^35'’),  and  the  polarization  ratio 
Qjjj(90“)/Qyy(90®);  the  scattering  measurement  angle  is  specified  for  each  quantity.  Each  of  these 


29. 


ratios  has  a  strong  dependence  on  particle  size  and  represents  a  redundant  set  of  data  in  terms  of  the 
particle  properties  for  particles  in  the  Mie  size  region  (D  >  X).  For  particles  in  the  Rayleigh  size 
region  (D  <  <  X),  scattering  is  isotropic  so  the  dissymmetry  ratio  has  a  value  of  unity  and  the 
defined  polarization  ratio  approaches  zero.  Thus,  in  this  small  particle  limit,  only  the 
scattering/extinction  ratio  yields  size  information.  However,  for  soot  formation  processes,  the  rapid 
coagulation  and  surface  growth  processes  lead  quickly  to  particles  in  the  Mie  size  region. 

With  the  particle  size  established,  any  one  of  the  scattering  cross  sections  can  be  used  to  find  N 
or  fy,  the  soot  volume  fraction.  The  soot  volume  fraction,  particle  size  and  number  density  are 
related  by 


NxD^ 

6 


(5) 


For  a  system  in  which  simultaneous  nucleation,  particle  growth  and  coagulation  are  present,  a 
particle  size  distribution  can  be  assumed  to  exist.  In  this  case,  the  previously  defined  particle 
scattering  cross  sections  must  be  averaged  over  the  size  distribution  function,  P(D),  to  yield  mean 
values. 
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A  widely  used  expression  for  P(D)  is  the  logarithmic  normal  distribution  that  is  given  by 

exp[-(ln  D/D/l2ah 
P(D)  =  - 2 - L 

^xoTgD 

where  Dg  and  are  the  geometric  mean  diameter  and  geometric  mean  standard  deviation, 
respectively. 
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33  Fuel  Molecular  Structure  Effects 


33.1  Introduction 

Fuel  molecular  structure  has  long  been  recognized  to  have  important  effects  on  the  production  of 
soot  particles  in  combustion  systems  [13,14].  In  particular  for  gas  turbine  engines,  the  advent  of 
broad  specification  and  alternative  fuels  has  heightened  the  interest  in  the  importance  of  fuel  structure 
effects  [15-18].  General  correlations  involving  the  carbon-to-hydrogen  (C/H)  ratio  have  been  shown 
to  grossly  represent  the  increasing  tendency  to  form  soot  as  the  fuel  hydrogen  content  decreases  [18]. 
Sooting  height  measurements  have  also  been  shown  to  strongly  depend  on  the  fuel  molecular  structure 
[14].  However,  the  direct  investigations  of  the  specific  processes  affected  by  changes  in  the  fuel 
structure  have  been  much  less  numerous,  particularly  for  mixing  dominated  conditions.  The  present 
study  has  focused  on  this  aspect  of  the  soot  formation  process  for  diffusion  flame  conditions. 

A  series  of  laminar  diffusion  flame  studies  was  conducted  to  examine  the  effects  of  fuel  molecular 
structure  on  the  formation  of  soot  particles.  In  order  to  maintain  similar  flame  conditions,  a  fuel 
mixture  sqjproach  was  employed.  This  fuel  mixture  approach,  which  involves  introducing  a  known 
amount  of  hydrocarbon  fuel  into  a  well  characterized  baseline  flame,  allows  the  soot  production  in  the 
flame  to  be  significantly  varied  while  maintaining  similar  flame  size  and  shape.  In  order  to  provide  a 
basis  for  comparison,  the  total  carbon  flow  rate  in  each  flame  is  kept  constant.  This  means  that  an 
appropriate  fraction  of  the  baseline  fuel  is  replaced  by  the  added  fuel  species.  With  the  soot 
contribution  from  the  baseline  fuel  known,  the  effect  of  the  added  fuel  on  soot  formation  can  be 
examined.  Methane  (CH4)  was  employed  as  the  baseline  fuel  because  it  displays  relatively  low  soot 
particle  formation  and,  thus,  provides  greater  sensitivity  for  measuring  changes  introduced  by  varying 
the  fuel  molecular  structure.  A  few  experiments  were  also  conducted  with  ethene  (C2H4)  as  the 
baseline  fuel  to  compare  the  effects  between  two  different  baseline  fuels. 

The  experiments  have  emphasized  the  methane  baseline  flames.  Table  2  summarizes  the  flame 
conditions  which  have  been  studied.  The  experiments  involving  the  fuel  addition  to  methane  and 
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ethene  flames  were  selected  to  examine  three  effects.  Specifically,  the  experiments  were  intended: 

(1)  To  systematically  examine  the  conversion  of  fuel  carbon  to  soot  as  function  of  fuel  flow  rate. 
The  experiments  in  which  ethene  or  butene  are  added  to  the  methane  flame  addressed  this 
point. 

(2)  To  investigate  the  relative  effects  of  changing  the  fuel  structure  in  a  series  of  fuels  which 
involved  an  alkane,  alkene,  and  alkyne  species.  The  experiments  in  which  butane,  1-butene 
and  1,3-butadiene  were  added  to  the  flame  addressed  this  point. 

(3)  To  compare  the  effects  of  the  baseline  fuel  on  the  conversion  of  fuel  carbon  to  soot  for  the 
fuel  addition  approach.  Experiments  in  which  ethene  and  1 -butene  have  been  added  to  both 
baseline  flames  have  been  utilized  for  this  comparison. 


Table  2.  Flow  conditions  for  the  fuel  addition  studies 


Baseline  Fuel 

Flow  Rate 
(cm^/s) 

Fuel  Added 

Flow  Rate 
(cm^/s) 

Air  Flow  Rate 
SCFM 

CH4 

8.75 

C2H4 

0.525 

2.75 

CH4 

7.7 

C2H4 

1.05 

2.75 

CH4 

5.6 

C2H4 

2.10 

2.75 

CH4 

5.6 

C4H8 

1.05 

2.75 

CH4 

5.6 

C4H8 

0.872 

2.75 

CH4 

7.7 

C4H8 

0.525 

2.75 

CH4 

5.6 

C4H10 

1.05 

2.75 

CH^ 

5.6 

1.05 

2.75 

The  analysis  which  is  presented  below  focuses  on  the  carbon  conversion  percentage  for  eadi  of 
the  flames  based  on  the  amount  of  fuel  added  to  the  baseline  flame.  The  contribution  of  the  baseline 
flame  is  subtracted  from  the  soot  present  in  the  flame  based  on  the  studies  of  the  atmospheric  ethene 
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and  methane  baseline  flames.  The  carbon  conversion  percentage  is  calculated  for  the  axial  location 
displaying  the  maximum  value  of  Fy,  the  integrated  soot  volume  fraction. 


R 

Fy  =  2x  J  fy  rdr 
0 


(9) 


where  R  is  the  radius  of  soot  particle  field.  The  carbon  conversion  percentage  can  then  be  expressed 
as 


%  conversion  = 


ihj.-mB 


100 


*hs  ~^sB 
dladd 


X  100 


(10) 


where  m^  is  the  soot  mass  flow  rate  at  the  location  of  maximum  Fy  for  the  flame  containing  the  fuel 
addition,  ihgB  >s  the  corresponding  soot  mass  flow  rate  for  the  pure  baseline  flame,  ihj.  is  the  mass 


flow  rate  of  carbon  entering  the  burner,  ihg  is  the  mass  flow  rate  of  carbon  contained  as  baseline  fuel 
and  ihmid  carbon  corresponding  to  the  fuel  species  added  to  the  flame.  The 

determination  of  the  soot  mass  flow  rate  at  a  particular  height  in  the  flame  can  be  calculated  from: 


R 

«hs(Zni)  =  2xp  J  v(r,  fy(r,  Zq,)  rdr 
o 


(11) 


where  z^,  is  the  axial  location  where  Fy  is  a  maximum,  p  is  the  density  of  soot  and  v  is  the  velocity. 
Thus,  the  velocity  profile  must  be  known  to  precisely  calculate  the  value  of  ihg.  For  the  present 
analysis,  the  velocity  has  been  assumed  to  be  independent  of  r  with  a  value  for  axial  location  Zj^  taken 
from  Ref.  3.  This  is  a  reasonable  assumption  for  the  locations  for  which  Fy  is  observed  to  be  a 
maximum  [3].  With  this  assumption  m^  can  be  expressed  as 

*s  (^m)  =  P  Fv(Zm) 

8 


where  p  is  taken  to  be  1.8  g/cm^.  Tables  3-5  tabulate  the  results  of  this  analysis  for  several  of  the 
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flames  studied.  Included  along  with  the  percent  conversion  result  is  the  calculated  adiabatic  flame 
temperature  for  each  flame. 

Table  3.  Flow  conditions  for  the  studies  of  the  fuel  flow  rate  effects 


Baseline  Fuel 

Flow  Rate 
(cm^/s) 

Fuel  Added 

Flow  Rate 
(cm^/s) 

%  Conversion* 

Tgd 

(K) 

CH4 

8.75 

C2H4 

0.525 

21.5 

2239 

CH4 

7.70 

C2H4 

1.05 

19.0 

2252 

CH4 

5.60 

C2H4 

2.10 

16.2 

2280 

CH4 

7.70 

C4H8 

0.525 

36.6 

2242 

CH4 

5.60 

C4H8 

0.872 

33.3 

2256 

CH4 

5.60 

C4H8 

1.05 

35.1 

2260 

*  p  -  1.8  gm/cm^ 


As  stated  earlier,  in  the  fuel  addition  approach,  a  fraction  of  the  baseline  is  replaced  with  the 
fuel  species  of  interest  subject  to  the  constraint  that  the  total  flow  rate  of  carbon  into  the  flame  is 
maintained  constant.  It  is  of  interest  to  ascertain  if  the  percentage  of  fuel  carbon  converted  to  soot  is 
dependent  on  the  amount  of  the  baseline  fuel  replaced.  The  results  tabulated  in  Table  3,  for  a 
methane  baseline  flame  in  which  ethene  or  butene  was  introduced,  provided  information  on  this  point. 
For  the  ethene  fuel  addition  case,  the  ethene  flow  rate  was  varied  by  a  factor  4,  whereas  for  the  1- 
butene  case  the  flow  rate  was  changed  by  a  factor  of  2.  In  the  case  of  the  1 -butene  studies,  the 
percent  conversion  remained  relatively  constant.  For  the  ethene  case,  a  small  systematic  decrease  in 
the  conversion  percentage  of  fuel  carbon  to  soot  is  observed  with  increasing  ethene  flow  rate. 
Consideration  of  the  variation  in  the  calculated  adiabatic  flame  temperature  (T^j)  does  not  explain  the 
observed  results  since  a  decrease  in  the  conversion  percentage  is  opposite  to  the  current  view  of  the 
effect  of  temperature  on  soot  production.  It  is  conceivable  that  the  observed  variation  in  the 


34. 


conversion  percentage  is  within  the  experimental  error.  At  present,  the  date  supports  the  conclusion 
that  the  conversion  percentage  is  only  weakly  dependent  on  the  fuel  addition  flow  rate.  The 
differences  observed  between  species  with  different  fuel  molecular  structure  are  significantly  greater 
than  the  observed  variations  with  respect  to  flow  rate.  This  provides  an  experimental  justification  to 
extend  the  results  of  these  studies  to  more  general  fuel  mixture  conditions. 

Table  4.  Flow  conditions  for  the  butane,  butene,  and  butadiene  addition  studies 


Baseline  Fuel 

Flow  Rate 
(cm^/s) 

Fuel  Added 
(cm^/s) 

Flow  Rate 

%  Conversion* 

T«i 

(K) 

CH4 

5.6 

C4H10 

1.05 

12.2 

2243 

CH4 

5.6 

C4H8 

1.05 

35.1 

2260 

CH4 

5.6 

C4H6 

1.05 

44.4 

2285 

•p  =  1-8  gm/cm^ 

To  specifically  examine  the  effect  of  fuel  molecular  structure  on  the  soot  formation  process,  a 
series  of  flames  were  examined  involving  C4  species.  The  flames  studied  involved  additions  of 
butane,  1-butene  and  1,3-butadiene  to  a  methane  baseline  flame.  In  this  series,  the  fuel  structure  is 
varied  in  terms  of  the  arrangement  of  carbon  bonds  (single  bonds  and  double  bonds)  resulting,  of 
course,  in  a  variation  of  the  carbon  to  hydrogen  ratio.  Table  4  tabulates  the  results  of  these 
experiments  in  terms  of  conversion  percentage  of  fuel  carbon  to  soot.  Figure  9  shows  the  integrated 
soot  volume  fraction,  Fy,  as  a  function  of  the  axial  coordinate,  z. 

Clearly,  the  results  in  Table  4  and  figure  9  indicate  a  strong  fuel  structure  effect  with  the 
conversion  percentage  more  than  tripling.  In  addition,  there  is  a  systematic  decrease  in  the  axial 
position  where  soot  is  first  observed  as  the  more  sooty  fuels  are  considered.  This  indicates  that  the 
reactions  leading  to  the  first  soot  particles  occur  more  rapidly  in  these  flames. 
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Figure  9.  Soot  volume  fraction  integrated  over  the  flame  diameter 


The  discussion  above  has  focused  on  the  effects  of  fuel  molecular  structure  and  flow  rate.  An 


additional  important  consideration  is  the  importance  of  the  baseline  fuel  which  is  used.  To  examine 
this  point,  two  fuel  addition  studies  involving  ethene  and  butene  were  compared  in  which  the  baseline 
fuels  were  varied  using  methane  and  ethene.  The  results  of  these  studies  in  terms  of  conversion 
percentage  are  tabulated  in  Table  S. 

Table  S.  Flow  conditions  for  flame  studies  of  baseline  fuel  synergism 


Baseline  Fuel 

Flow  Rate 
(cm^/s) 

Fuel  Added 

Flow  Rate 
(cm^/s) 

%  Conversion* 

Tad 

(K) 

CH4 

7.7 

C2H4 

1.05 

19.0 

2252 

CH4 

7.7 

C4H8 

0.525 

36.6 

2242 

C2H4 

3.85 

C2H4 

1.05 

16.2 

2369 

C2H4 

3.85 

C4H8 

0.525 

46.8 

2359 

•p  =  1.8  gmIcTV? 

These  results,  as  with  the  flow  rate  results  previously  discussed  (see  Table  3),  do  show  some 
sensitivity  to  a  variation  in  the  flame  conditions.  However,  again,  the  differences  are  significantly 
smaller  than  the  variation  resulting  from  the  fuel  structure  variation.  For  the  experiments  shown  in 
Table  S,  the  soot  conversion  percentage  approximately  doubles  for  the  change  in  fuel  species. 
However,  differences  in  the  results  for  the  various  baseline  flames  for  a  particular  fuel  (ethene  or 
butene)  is  typically  20-25%.  It  is  also  worth  mentioning  that  the  calculated  flame  temperatures  vary 
by  about  lOOK  for  the  fuels  studied.  Thus,  some  of  the  variation  may  be  a  result  of  the  temperature 
difference.  The  comparison  between  the  ethene  and  butene  results  may  indicate  different  temperature 
sensitivities  since  the  observed  effects  of  the  baseline  flame  are  reversed. 
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The  above  discussion  has  emphasized  the  conversion  percentage  of  fuel  carbon  to  soot  particles. 
This  quantity  is  very  useful  in  illustrating  the  impact  of  the  effect  of  fuel  molecular  structure  on  soot 
formation.  However,  it  represents  a  global  measurement  of  the  soot  formation  process.  More 
fundamentally  significant  results  are  realized  when  the  detailed  temperature  time  history  characterizing 
these  chemically  reacting  systems  are  analyzed. 

ZJiJl  Soot  Particle  Precursor  and  Inception  Studies 

In  order  to  examine  the  details  of  the  soot  formation  process,  it  is  useful  to  follow  the  temporal 
evolution  of  the  soot  particle  field.  Thus,  measurements  resolved  along  specific  particle  paths  are  of 
more  use  in  developing  quantitative  results.  From  measurements  of  the  axial  and  radial  velocity 
components,  the  necessary  data  to  provide  this  information  can  be  obtained  [3]. 

Laser  light  scattering  measurements  for  soot  particle  characterization  and  fluorescence 
measurements  to  characterize  the  soot  precursor  field  prior  to  particle  inception  were  obtained  in 
butane,  1-butene  and  1,3  butadiene  mixture  studies  described  above.  The  fluores«.mce  measurements 
were  obtained  by  tuning  the  laser  to  the  488  nm  laser  line  and  observing  the  fluorescence  a  514.5  nm. 
Previous  workers  have  shown  that  this  fluorescence  likely  results  from  large  hydrocarbon  species, 
probably  aromatic  in  nature,  which  are  present  in  the  flame  [20-23].  Comparisons  of  the  observed 
fluorescence  intensity  as  the  fuel  is  varied  provide  a  semi-quantitative  indication  of  the  effects  on  the 
gas  phase  soot  precursors.  Unfortunately,  such  measurements  are  not  species  specific  and,  thus,  can 
not  add  direct  insight  into  the  chemical  mechanism  leading  to  soot  particle  inception. 

The  results  of  these  fluorescence  and  soot  volume  fraction  measurements  are  shown  in  figure  10. 
Figure  10  displays  the  results  along  the  particle  path  exhibiting  the  maximum  soot  formation.  The 
fluorescence  measurements  are  observed  to  precede  the  formation  of  soot  particles  and  then  to 
increase  rapidly  in  the  early  formation  period.  The  magnitude  of  the  fluorescence  signals  is  observed 
to  be  strongly  dependent  on  the  fuel  molecular  structure  and  vary  in  a  manner  very  similar  to  the  soot 
volume  fraction  as  the  added  fuel  is  varied.  Comparisons  along  individual  particle  paths  indicate  that 
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Figure  10.  The  soot  volume  fraction  and  fluorescence  measurements  along  the  particle  paths  containing 
the  maximum  soot  concentration  in  the  annular  region  of  the  flame. 


between  20%  and  70%  of  the  soot  growth  occurs  prior  to  the  fluorescence  achieving  a  maximum 
value  depending  on  the  particle  path  location  in  the  flame.  This  point  raises  the  question  as  to  the 
participation  of  large  hydrocarbon  species  in  the  surface  growth  process  since  such  species  are 
responsible  for  the  observed  fluorescence.  This  point  will  be  addressed  again  later  in  this  report. 

The  results  described  above  lead  to  the  following  observations: 

1 .  The  species  responsible  for  the  observed  fluorescence  increase  in  proportion  to  the  sooting 
tendency  of  the  fuels  studied.  In  this  respect,  such  species  satisfy  one  criteria  for  representing 
the  soot  precursors. 

2.  These  fluorescing  species  precede  the  particle  inception  region  throughout  the  flame  zone. 

This  is,  of  course,  a  minimal  condition  for  a  precursor  species. 

3.  Following  particle  inception,  these  species  increase  rapidly  during  a  significant  fraction  of  the 
particle  surface  growth  process.  This  leads  to  the  possibility  that  such  species  contribute  to 
the  early  surface  growth  process  as  well  as  the  inception  process. 

The  results  described  above  are  quite  consistent  with  results  obtained  recently  in  premixed  flames 
with  respect  to  the  likely  importance  of  PAH  species  on  soot  particle  inception  [32,33].  Although 
some  work  has  been  obtained  under  diffusion  flame  conditions,  present  results  lack  the  temporal  and 
spatial  detail  necessary  to  compare  to  models  of  diffusion  flames  [2,20,23,24].  As  the  models 
successfully  developed  for  premixed  flames  are  extended  to  mixing  dominated  combustion  conditions, 
there  will  be  a  need  for  a  more  exhaustive  measurement  base  to  provide  appropriate  comparisons  with 
model  prediction.  Recent  experiences  with  modeling  of  soot  particle  surface  growth  and  oxidation 
has  shown  the  value  of  detailed  laminar  diffusion  flame  measurements  for  model  validation  [1,25- 
29,30].  A  similar  effort  in  the  area  of  soot  precursor  and  particle  inception  phenomena  is  certainly 
warranted.  Results  similar  to  those  described  above  in  figure  10  can  be  used  presently  for 
comparisons  with  models  for  total  PAH  evolution.  However,  such  comparisons  would  benefit  greatly 
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from  additional  measurements  of  stable  as  well  as  radical  species  in  the  flame.  This  should  be  a 
priority  area  for  research  involving  soot  formation  in  diffusion  flames. 

3.3  J  Surface  Growth 

Although  particle  inception  is  one  of  the  key  processes  in  the  formation  of  soot,  it  actually 
accounts  only  for  a  small  fraction  of  the  mass  addition.  The  major  mass  addition  process  involves  the 
surface  growth  stage  which  follows  inception.  A  great  deal  of  attention  has  been  given  to  this  process 
in  both  premixed  and  diffusion  flames  [5,31,32],  Major  contributions  to  the  present  understanding  of 
surface  growth  have  been  made  by  Wagner  [48]  and  Harris  [31].  The  mechanisms  developed  in  these 
studies  differ  significantly  in  concept  and  work  is  continuing  to  further  understand  the  surface  growth 
process  [37-43]. 

Most  approaches  tc  soot  surface  growth  phenomena  argue  that  acetylene,  or  a  species  closely 
associated  with  acetylene,  provides  the  major  surface  reactant.  Simple  mass  balance  considerations 
provide  the  basis  for  this  observation,  since  only  acetylene  appears  to  be  present  in  sufficient 
concentration  to  account  for  the  amount  of  mass  added  to  the  soot  particles  [31].  However,  other 
workers  have  pointed  out  that  under  suitable  conditions  large  PAH  species  may  directly  add  to  the 
growing  soot  particles  [32-34].  Since  measurements  made  in  diffusion  flame  are  limited  to  a  small 
number  of  fuels,  it  is  difficult  to  presently  assess  the  importance  of  these  PAH  species  to  surface 
growth  processes  in  soot  formation. 

Before  proceeding  further  it  is  worthwhile  to  review  the  current  understanding  of  soot  particle 
surface  growth  process  in  laminar  diffusion  flames.  In  particular,  the  effect  of  fuel  structure  on  the 
amount  of  soot  formed  in  the  flame  is  useful  in  this  discussion. 

The  present  work  at  Penn  State  investigating  laminar  diffusion  flames  has  shown  that  the  fuel 
structure  effects  can  be  taken  into  ac  unt  if  variations  in  the  available  surface  area  are  considered 
[4,5,35].  This  analysis  implies  that  the  initial  surface  area  formed  during  the  particle  inception 
process  controls  subsequent  growth.  Returning  to  the  results  shown  in  figure  10,  the  analysis  can  be 
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extended  to  consideration  of  the  surface  growth  and  oxidation  regions  of  the  flame.  The  discussion 
above  focused  on  fuel  structure  effects  on  particle  inception  where  clearly  a  strong  sensitivity  to  the 
fuel  is  observed.  It  is  equally  interesting  to  determine  the  effect  of  fuel  molecular  structure  on  the 
surface  growth  process  for  soot  particles  formed  in  diffusion  flames.  An  additional  point  of  interest 
lies  in  the  role  that  large  hydrocarbon  species  may  have  in  the  early  surface  growth  process.  This 
aspect  arises  from  the  results  described  above,  in  which  the  fluorescence  signals  attributed  to  large 
PAH  species  maximized  at  some  significant  fraction  of  the  mass  addition  to  the  soot  particles,  thus, 
allowing  the  possibility  that  these  species  contribute  significantly  to  early  soot  growth. 

The  analysis  has  been  carried  out  along  the  same  particle  path  for  which  measurements  are  shown 
in  figure  10.  Similar  results  have  also  been  obtained  for  the  centerline  region  of  these  flames  and 
discussion  of  those  results  is  also  included  below.  The  present  approach  to  the  analysis  is  to  represent 
the  soot  mass  growth  rate  as 

^  =  k  S[C*  Hy]  (13) 

where  m  is  the  mass  of  soot  (g/cm^),  k  is  the  surface  growth  rate  constant  (g/cm^-s-atm),  S  is  the 
specific  surface  area  (cm^/cm^)  and  [Cjj  Hy]  is  the  species  concentration  reacting  with  the  surface  to 
add  mass  (atm).  In  the  experiments,  the  product  of  k  and  [C^  Hy]  can  be  determined  since  S  and 
dm/dt  are  determined  from  the  combined  light  scattering  and  velocity  measurements.  This  product, 
k[C,H  y],  is  termed  the  specific  surface  growth  rate  constant.  The  specific  surface  area  is  calculated 
from  the  particle  diameter  and  number  density  measurements,  while  the  soot  mass  growth  rate  can  be 
found  from  the  measured  soot  volume  fraction  assuming  the  density  of  soot  particles  is  1.8  g/cm^. 
Figures  1 1  and  12  show  the  specific  surface  growth  rate  constant  for  the  C4  fuel  studies  in  the 
aimulus  and  along  the  centerline,  respectively.  Shortly  after  the  particles  are  formed,  the  specific 
surface  growth  is  close  to  2x10'^  g/cm^-s.  After  a  period  for  which  the  rate  constant  is 
approximately  constant,  the  surface  growth  rate  drops  rapidly.  It  should  be  pointed  out  that  the 
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Figure  11.  The  specific  surface  growth  rate  as  a  function  of  time  for  the  butane,  1 -butene  and  1,3 

butadiene  addition  studies  to  the  methane  baseline  flame.  Results  are  for  the  particle  path 
containing  the  maximum  value  of  the  soot  volume  fraction. 
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Figure  12.  The  specific  surface  growth  rate  constant  as  a  function  of  time  for  the  butane,  1-butene  and 
1,3  butadiene  addition  studies  to  the  methane  baseline  flame.  Results  are  for  the  particle  path 
along  the  centerline  of  the  flame. 


present  analysis  has  been  limited  to  the  region  of  the  flame  before  the  soot  volume  fraction  reaches  a 
maximum.  Thus,  regions  where  soot  particle  oxidation  dominates  have  been  excluded  from 
consideration. 

Figure  13  shows  the  results  for  the  ethene  baseline  flame  studies  where  ethene,  propene,  butene 
and  toluene  were  added  to  the  flame.  Clearly  the  same  behavior  is  observed  in  these  flames.  Also 
shown  in  figure  13  are  the  results  of  a  rich  ethene/air  premixed  flame  studied  by  Harris  and  co- 
workers  [31].  It  is  interesting  to  note  that  although  the  temporal  behavior  of  the  specific  surface 
growth  rate  differs,  there  is  order  of  magnitude  agreement  between  the  value  for  this  rate  constant  in 
both  flame  systems. 

From  these  extensive  studies,  it  is  clear  that  the  surface  growth  process  is  similar  in  all  the 
flames  studied  here.  Furthermore,  since  accounting  for  the  differences  in  available  surface  area 
substantially  explains  differences  in  the  growth  processes  for  flames  of  widely  differing  soot 
concentrations,  the  value  of  the  initial  surface  area  appears  to  be  critical  to  subsequent  surface 
growth.  Additionally,  these  studies  show  that  similar  values  of  specific  surface  growth  rate  are  found 
throughout  the  flame  for  the  fuels  studied. 

Before  making  further  observations  concerning  the  present  status  of  understanding  in  the  area  of 
soot  particle  surface  growth,  some  discussion  of  the  effect  of  the  fluorescing  species  in  this  process  is 
appropriate.  Clearly,  if  the  initial  surface  area  is  important  in  all  subsequent  growth,  species  involved 
in  the  first  stages  of  surface  mass  addition  should  be  carefully  examined.  The  results  described 
earlier  regarding  soot  particle  inception  showed  that  the  fluorescence  measurements  attained  a 
maximum  after  some  significant  fraction  of  the  mass  due  to  surface  growth  had  occurred.  To 
quantify  the  contribution  of  large  hydrocarbon  species  attributed  as  the  source  of  the  observed 
fluorescence,  a  series  of  highly  spatially  resolved  radial  profiles  of  the  soot  particle  field  were 
obtained  for  closely  spaced  intervals  along  the  axis  of  the  flame.  These  profiles  correspond  to  the 
earliest  data  points  shown  in  figure  11.  With  the  exception  of  the  butadiene  flame,  no  systematic 
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increase  in  the  specific  surface  growth  rate  constant  is  observed.  Thus,  the  present  data  do  not 
support  a  model  in  which  the  species  observed  in  the  fluorescence  measurements  contribute 
extensively  to  the  surface  growth  in  the  region  immediately  following  particle  inception.  However, 
the  fluorescence  results  represent  a  global  measurement  of  the  large  hydrocarbon  species  present  in 
the  flame.  Recent  work  in  premixed  flames  by  Harris  [36]  and  McKinnon  [34]  supports  a  more 
careful  consideration  of  this  matter. 

Finally,  the  results  shown  in  figures  11  through  13  indicate  a  decrease  in  the  specific  surface 
growth  rate  constant  as  reaction  time  proceeds.  This  observation  has  been  noted  for  both  premixed 
and  diffusion  flames.  To  date,  however,  there  is  no  satisfactory  explanation  for  this  effect.  In 
general,  three  mechanisms  can  be  considered  as  likely  candidates.  First,  the  hydrocarbon  growth 
species  could  be  depleted  resulting  in  decreased  surface  growth.  Measurements  in  fuel-rich  premixed 
flames  do  not  support  this  mechanism,  however,  a  similar  data  base  for  diffusion  flames  is  lacking. 

A  second  mechanism  involves  a  decrease  in  the  particle  surface  reactivity.  In  this  process,  as  the 
particle  ages  in  the  high  temperature  flame  environment,  a  loss  of  radical  reactive  sites  (active  sites) 
on  surface  results,  leading  to  a  decrease  in  the  surface  reactivity.  Finally,  in  the  case  of  diffusion 
flames,  there  is  a  distinct  possibility  that  oxidation  begins  to  compete  effectively  with  soot  growth, 
causing  the  net  growth  rate  to  approach  zero  or  even  become  negative.  Previous  studies  have  shown 
that  the  emission  of  soot  particles  from  laminar  diffusion  flames  is  strongly  affected  by  temperature, 
which  soot  particles  influence  through  radiative  transfer  [12,24].  The  specific  local  conditions  which 
govern  the  transition  from  growth  to  oxidative  destruction  is,  however,  presently  not  known.  An 
understanding  of  the  transition  region  between  the  cessation  of  soot  mass  growth  and  the  onset  of 
oxidation  would  have  a  significant  impact  on  predictive  modeling  efforts. 

Recently  there  has  been  increased  interest  in  the  specific  phenomena  controlling  the  surface 
reactivity  in  soot  surface  growth  processes  [32,33,37-40].  Each  of  these  studies  has  focussed  on  the 
concept  of  active  surface  sites  as  reaction  centers  for  the  addition  of  carbon  through  reactions  with  gas 
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SPECIFIC  SURFACE  GROWTH 


Figure  13.  The  specific  surface  growth  rate  along  the  particle  path  exhibiting  the  maximum  soot  volume 
fraction.  The  dashed  line  ( — )  is  the  best  faired  curve  the  ethene  flame  results. 
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phase  species  such  as  acetylene.  The  objectives  of  these  studies  differ,  but  they  do  raise  a  set  of 
common  features.  First,  soot  particles  possess  highly  reactive  surfaces,  exceeding  by  two  orders  of 
magnitude  the  reactivity  of  pyrolytically  formed  carbon  particles  (pyrocarbon)  [39].  The  regeneration 
of  active  sites  appears  to  be  related  to  atomic  hydrogen  reactions  with  the  surface  to  form  new 
reactive  sites  [32,33],  which  have  a  relativity  constant  concentration  in  many  situations  [38,40].  This 
approach  has  been  used  to  reconcile  the  disparate  observations  in  the  literature  regarding  the 
sensitivity  of  the  surface  growth  process  to  the  available  soot  particle  surface  area.  The  results  of 
Harris  and  Weiner  [31,41]  indicate  a  direct  proportionality  to  available  surface  area  as  indicated  by 
eq.  13.  Results  from  the  work  of  Haynes  et  al.  [42]  and  Weischnowsky  et  al.  [43]  indicate  that  the 
surface  growth  process  is  independent  of  the  available  surface  area.  In  the  recent  work,  Harris  [38], 
Haynes  [40]  and  Frenklach  [32]  have  replace  or  augmented  the  surface  area  arguments  with  the 
addition  of  a  surface  active  site  concentration.  Although  the  details  of  the  representations  are 
different,  the  basic  approach  involves  similar  concepts  with  regard  to  the  importance  of  these  surface 
sites.  For  example,  Haynes  [40]  represents  the  increase  in  soot  volume  fraction  with  time  as: 


dfy  g 

—  =a-  N,A«v„N 


dt 


s 


(14) 


where  fy  is  the  soot  volume  fraction,  a  is  the  growth  reaction  probability  on  collision,  Ng  is  the 
number  density  of  gas  phase  growth  species,  Vg  is  the  molecular  volume  of  the  growth  species  when 
incorporated  in  the  soot  particle,  c  is  the  mean  free  speed  of  the  gas  phase  growth  species  (8kT/x' 
mg)*^^,  Ng  is  the  total  concentration  of  active  sites  and  Ag  is  the  active  area  of  a  single  site.  In  this 
representation  the  difference  in  the  experimental  results  is  determined  by  the  behavior  of  the  active 
sites  in  the  various  flames.  In  the  results  of  Harris  and  Weiner  [31,41]  the  active  site  concentration 
remains  nearly  constant  and,  thus,  scales  with  particle  surface  area.  In  contrast,  for  the  higher 
temperature  flame  studies  by  Haynes  et  al.  [42]  and  Weischnowsky  et  al.  [43],  the  surface  site 
concentration  is  determined  to  decrease  with  time.  The  basis  for  the  decay  in  the  active  surface  site  is 


48. 


«  > 


not  known,  although  the  higher  temperature  flames  can  be  argued  to  provide  an  environment  in  which 
more  rapid  annealing  of  the  soot  particles  may  be  present.  Frenklach  [32]  has  argued  in  a  somewhat 
similar  manner  that  the  state  of  the  initial  surfaces  of  the  soot  particles  may  differ  in  the  two  flames 
with  the  higher  temperature  flames  having  fewer  active  sites. 

The  above  work  illustrates  the  detail  to  which  present  soot  models  and  experiments  can  consider 
the  fundamental  mechanisms  controlling  individual  processes.  Although  this  modeling  capability  is 
impressive,  there  is  little  data  particularly  under  diffusion  flame  conditions  to  compare  with  model 
prediction.  Clearly,  the  results  of  the  present  project  point  to  the  need  for  further  measurements 
obtained  under  conditions  which  allow  detailed  comparisons  with  mechanistic  models  for  surface 
growth. 

Despite  the  questions  concerning  the  specifics  of  the  surface  growth  mechanism,  recent  efforts  to 
model  soot  particle  growth  in  diffusion  flames  have  met  with  some  success  in  reproducing  the  general 
features  of  the  soot  formation  process  [26].  However,  the  model  representation  for  the  surface 
growth  process  shows  a  strong  sensitivity  to  the  specific  surface  growth  rate  constant.  In  this  case, 
small  differences  in  the  growth  rate  constant  could  be  argued  to  be  responsible  for  variations  in  the 
sooting  propensity  of  various  fuels.  In  the  discussion  of  the  results  shown  in  figures  1 1  through  13, 
emphasis  was  given  to  the  importance  of  the  available  surface  area  and  the  critical  nature  of  the 
inception  process  which  controls  initial  soot  surface  area.  Thus,  we  have  argued  along  with  others 
[44]  that  inception  is  controlling  the  amount  of  soot  ultimately  formed  in  the  flame.  From  the 
modeling  studies  of  Kennedy  [26],  this  conclusion  has  been  challenged  based  on  the  observation  that 
in  heavily  sooting  flames  such  as  those  involving  ethene  [2,3]  surface  growth  controls  the  amount  of 
soot  formed.  Since  these  modeling  studies  used  the  ethene  data  shown  in  figure  13  to  validate  the 
model,  the  observations  by  Kennedy  provide  a  good  justification  for  carefully  reconsidering  the 
analysis  of  these  flames  and  assessing  the  need  for  further  study  of  surface  growth  in  laminar  flames. 
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In  fact,  it  is  clear  that  more  comparisons  with  the  model  for  a  wider  range  of  fuels  would  be 
extremely  valuable.  For  such  comparisons,  not  only  must  the  soot  particle  field  be  determined,  but 
the  major  gas  phase  species  also  need  to  be  determined.  This  latter  requirement  is  due  to  the  fact  that 
the  Kennedy  model  uses  a  measured  mixture  fraction  relationship  for  each  fuel  studied  [26]. 

3.4  Concentration  and  Temperature  Effects 

A  detailed  study  has  been  made  of  the  relative  effects  of  concentration  and  temperature  on  soot 
formation  in  diluted  laminar  ethene/air  diffusion  flames  employing:  (1)  argon  and  nitrogen  diluted 
flames  at  equal  calculated  adiabatic  flame  temperatures  but  different  initial  concentrations  to  isolate 
the  effect  of  concentration  and  2)  argon  and  nitrogen  diluted  flames  with  equal  initial  fuel 
concentrations  but  different  temperatures  to  isolate  the  effect  of  temperature.  Total  integrated  soot 
volume  fraction  measurements  show  that  in  argon  and  nitrogen  diluted  flames  with  equal  calculated 
adiabatic  flame  temperatures,  the  more  diluted  argon  flames  consistently  display  lower  soot 
concentrations.  However,  in  flames  of  equal  dilution,  argon  flames  consistently  display  higher  soot 
concentrations  than  the  slightly  cooler  nitrogen  diluted  flames.  Local  temperature  measurements  show 
that  for  nitrogen  and  argon  diluted  flames  with  equal  calculated  adiabatic  flame  temperature,  argon 
diluted  flames  display  lower  temperatures  in  the  region  where  soot  is  first  formed.  Mass 
spectrometric  measurements  of  gas  concentrations  in  diluted  and  undiluted  methane/air  diflusion 
flames  were  obtained  and  compared  with  a  numerical  flame  model.  These  results  show  that  the  initial 
difference  in  fuel  concentration  in  diluted  and  undiluted  flames  diminishes  rapidly  with  height. 
Furthermore,  laser  light  scattering  measurements  show  that  as  inert  diluent  is  added  to  the  flame,  soot 
inception  is  delayed  and  consequently  less  time  is  available  for  soot  growth.  Based  on  this  extensive 
set  of  data,  a  quantitative  assessment  was  made  of  the  relative  effects  of  temperature  and 
concentration  indicating  that  for  coflow  diffusion  flames,  the  role  of  temperature  is  more  important 
than  the  decrease  in  reaction  rate  due  to  the  t  eduction  in  fuel  concentration  when  an  inert  diluent  is 
added  to  the  fuel  flow. 


50. 


A  complete  treatment  of  the  results  of  these  studies  has  been  submitted  to  Combustion  Science 
and  Technology  for  publication  and  is  attached  as  Appendix  1. 

3.5  Analysis  of  Light  Scattering  From  Soot  Particles  Using  Optical  Cross  Sections  for 
Aggregates 

Soot  formed  in  flames  usually  consists  of  aggregates  (clusters  or  agglomerates)  of  a  variable 
number  of  nearly  spherical,  monodisperse  primary  particles  (monomers  or  spherules).  In  this  work, 
the  optical  properties  of  polydisperse  aggregates  are  used  to  analyze  light  scattering  data  from  a 
coannular  ethene  diffusion  flame.  In  previously  reported  studies,  data  have  been  obtained  on  the  local 
extinction  and  volumetric  scattering  cross  sections  from  laser  scattering  experiments,  on  the  flame 
velocity  field  from  laser  velocimetry,  and  on  the  primary  particle  sizes  determined  by  electron 
microscopy.  The  present  analysis  yields  the  average  number  of  primary  particles  per  aggregate,  the 
mean-square  radius  of  gyration,  the  soot  volume  fraction  and  the  aggregation  rate.  It  is  found  that 
sustained  collisional  growth  of  the  aggregates  occurs  while  their  primary  particles  grow  through 
heterogeneous  reactions  low  in  the  flame,  and  contract  through  surface  oxidation  in  the  upper  half  of 
the  flame.  A  recent  value  of  the  refractive  index  gives  internally  consistent  moment  ratios  of  the 
aggregate  size  distribution  function.  This  method  of  analysis  provides  a  more  detailed  and  complete 
description  of  the  formation,  growth  and  oxidation  of  soot  aggregates  in  a  diffusion  flame. 

A  complete  treatment  of  the  problem  of  aggregate  structure  on  interpretation  of  soot  particle 
properties  is  contained  in  Appendix  2.  This  appendix  contains  a  reprint  of  a  recent  article  published 
in  the  Twenty-Third  Symposium  (International)  on  Combustion. 

3.6  Operating  Pressure  Effects  on  Soot  Formation 

A  series  of  laminar  diffusion  flames  have  been  studied  over  a  range  of  pressures  between  1  and 
10  atm  (.1  MPa  to  1.0  MPa).  The  coannular  diffusion  flame  burner  and  laser  diagnostics  for  soot 
particle  measurements  were  similar  to  those  used  for  atmospheric  flame  studies  previously  reported. 
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The  coannular  burner  has  been  mounted  in  a  high  pressure  cell  which  was  extensively  described 
earlier.  The  studies  undertaken  involved  laminar  dif^sion  flames  burning  ethene  (C2H4),  ethane 
(C2Hg)  and  a  mixture  of  ethane  and  propene  (C3H6)  in  air.  The  measured  soot  volume  fraction  is 
obtained  from  laser  extinction  measurements  and  is  represented  by  [45]: 

-j*  fvdx  =  C(X,m)ln(I/Io)  (15) 

where  is  the  local  soot  volume  fraction,  x  lies  along  a  laser  path  through  the  flame  on  a  diameter, 
C(X,m)  is  a  constant  which  can  be  determined  from  Rayleigh  light  scattering  theory  to  be 
1.05x10*5  for  m=  1.57-0.56  i  at  X=514.5nm  [45],  and  finally  I/Iq  is  the  ratio  of  the  transmitted  to 
incident  laser  power.  This  quantity  is  proportional  to  the  total  soot  present  at  a  particular  axial 
location  in  the  flame.  Figure  14  shows  this  quantity  as  a  function  of  the  non-dimensional  axial 
location  for  the  ethene/air  flames  at  a  series  of  operating  pressures.  The  non-dimensional  axial 
position  is  expressed  as 

=  (ZD/Q)  In  (1  +  1/s)  (16) 

where  Z  is  the  axial  position,  D  is  the  diffusion  coefficient,  Q  is  the  volumetric  fuel  flow  rate  and  S 
is  the  volume  of  air  to  the  volume  of  fuel  required  for  complete  combustion. 

If  the  increase  in  the  maximum  soot  volume  fraction  at  a  particular  operating  pressure  is  related 
to  that  pressure  through  a  power  law  dependence,  that  is 

(f  fv'llWoP" 

then  a  fit  to  the  data  can  be  used  to  yield  a  value  for  n.  Figure  15  shows  the  maximum  soot  volume 
fractions  as  a  function  of  the  operating  pressure  along  with  the  best  fit  value  for  n  determined  from  a 
linear  least  square  procedure.  For  the  ethene  flame,  the  fuel  flow  rate  was  3.9  cm  /s  and  n  was 
found  to  be  1.05±  .06  in  good  agreement  with  the  results  of  Flower  and  Bowman  [45]  who  found 
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Figure  14.  C(X,in)  In  (I/Iq)  as  a  function  of  the  non-dimensional  position,  ri,  for  the  ethene/air  laminar 
diffusion  flames  for  operating  pressures  between  1  and  10  atm.  The  non-dimem  ional  axial 
position  is  given  by  jj  =  (zD/Q)  In  (1  +  1/s),  where  z  is  the  axial  position,  D  is  the  diffusion 
coefficient  (0.156  cm^s**),  Q  is  the  fuel  flow  rate  (3.9  cm^/s)  and  s  is  the  volume  of  air  to 
the  volume  of  fuel  required  for  complete  combustion.  (S=  14.28  for  C2H4). 


C(X.m)  ln(l/lo)mQx 


1  .OOOE' 

1  .OOOE 

5.000E 


Figure  15.  C(X,in)  In  (I/Io)/inax  “  ^  function  of  operating  pressure  for  a  series  of  ethene,  ethane  and 
ethane/propene  laminar  diffusion  flames  burning  in  air.  Solid  lines  are  fits  to  P". 


n=  1. 2+0.1.  It  should  be  noted  that  results  shown  in  figures  14  and  15  have  taken  into  account  the 
slight  burner  diameter  difference  for  the  two  studies.  In  the  present  study,  the  diameter  was  1 . 1  cm 
where  as  Flower  and  Bowman  employed  a  burner  with  a  1.27  cm  diameter.  Also  shown  on  figure  15 
are  the  results  for  an  ethene/air  diffusion  flame  where  a  copper  (Cu)  rather  than  a  stainless  steel  (SS) 
fuel  tube  was  utilized.  The  results  are  very  similar  for  both  fuel  tube  materials.  Although  the  higher 
heat  loss  characteristics  of  the  copper  fuel  tube  were  observed  to  slightly  reduce  the  amount  of  soot 
formed  as  compared  to  the  stainless  steel  case,  the  pressure  dependence  is  largely  unaffected  by  the 
fuel  tube  material. 

Figure  15  also  shows  the  results  obtained  for  flames  burning  ethane  or  an  ethane/propene 
mixture.  For  these  flames,  the  pressure  could  only  be  varied  over  a  range  of  1  to  3  atm  (.1  to  .3 
MPa)  before  a  buoyancy  driven  flame  instability  was  observed.  These  buoyancy  effects,  at  elevated 
pressure,  have  recently  been  the  subject  of  a  collaborative  investigation  with  workers  at  the  National 
Institute  of  Standards  and  Technology  (NIST)  [46]  and  are  discussed  elsewhere  in  this  report.  An 
analysis  similar  to  that  applied  to  the  ethene  flame  studies  yields  pressure  power  dependences  of 
n=1.88±0.13  and  n=1.62±0.24  for  ethane  fuel  flow  rates  of  2.80  and  3.85  cm^/s,  respectively. 
These  higher  values  of  the  pressure  power  dependence  appear  to  be  related  to  fuel  structure  effects, 
the  direct  nature  of  which  remains  to  be  understood.  These  observations  of  a  fuel  structure 
dependence  are  consistent  with  previous  studies  [47]. 

As  a  further  investigation  of  fuel  molecular  structure  effects,  a  mixture  of  ethane  (2.80  cm^/s) 
and  propene  (0.7  cm^/s)  was  studied  and  these  results  are  also  shown  in  figure  15.  In  the 
ethane/propene  study,  the  total  carbon  atom  fuel  flow  rate  was  maintained  constant  with  respect  to  the 
3.85  cm^/s  ethane  flame.  For  these  conditions,  the  value  of  n  is  observed  to  be  1.16+0.10,  a  value 
similar  to  that  observed  for  the  ethene  studies.  Thus,  the  present  studies  with  an  alkane  species 
appear  to  yield  a  higher  pressure  dependence  than  observed  for  either  of  the  alkene  species.  This 
occurrence  of  a  variation  in  the  pressure  dependence  may  be  important  in  providing  an  understanding 
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of  the  differences  observed  in  premixed  and  diffusion  flame  studies  at  high  pressure.  For  premixed 
flames,  pressure  dependence  studies  have  yielded  values  of  n=2  [48]  as  compared  to  the  values  of  n 
between  1  and  1.2  reported  for  ethene.  The  higher  values  of  n  for  the  ethane  flames  are  closer  to  the 
values  observed  in  premixed  flames.  The  cause  of  the  higher  pressure  dependence  obtained  for  the 
ethane  diffusion  flames  remains  to  be  established.  Fuel  structure  effects  may  be  reflected  in  the  lower 
soot  volume  fractions  rather  than  explicitly  different  chemical  mechanisms.  Premixed  flame  studies 
have  examined  a  variety  of  fuels,  including  ethene  and  benzene,  and  obtained  a  dependence  on 
pressure.  The  present  experiments  provide  some  important  additional  information  for  consideration  in 
resolving  the  actual  pressure  dependence  for  soot  formation  processes  and  indicate  further  study  of 
high  pressure  diffusion  flames  is  required. 

3.7  Isolation  of  Buovancv  Effect.s  in  .let  Diffusion  Flame  Experiments 

The  occurrence  of  instability  phenomena  in  flames  has  been  a  topic  of  interest  in  combustion  for 
a  variety  of  reasons.  Instability  processes  have  been  observed  to  enhance  mixing  and,  thus,  accelerate 
chemical  reaction  progress.  Less  favorable  results  have  been  observed  to  accompany  combustion 
instabilities  as  well.  For  gas  turbine  and  rocket  engines,  combustion  instabilities  can  present  serious 
operational  problems  and  may,  in  fact,  lead  to  catastrophic  engine  failure.  In  such  cases,  expensive 
engine  redesign  or  implementation  of  damping  devices  has  been  necessary  to  resolve  the  problem. 
Clearly,  a  desirable  situation  would  be  one  in  which  the  potential  benefits  to  enhanced  mixing  could 
be  realized  without  the  penalties  of  hardware  performance  loss.  This  requires  a  fundamental 
understanding  of  phenomena  which  controls  the  initiation  and  growth  of  the  instability.  One  of  the 
obstacles  to  achieving  such  an  understanding  lies  in  the  inability  to  study  instability  phenomena  under 
well  controlled  conditions. 

In  addition  to  the  combustion  instability  interests,  turbulence  combustion  phenomena  may  also  be 
influenced  by  instability  phenomena.  The  occurrence  of  well  ordered  structures  in  turbulent  flames 
has  been  a  subject  of  great  interest  over  the  preceding  decade.  Although  often  observed,  the 
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importance  of  such  features  with  respect  to  the  transition  to  turbulent  combustion  conditions  is  still  an 
area  of  active  investigation.  An  understanding  of  the  basic  mechanisms  controlling  the  onset  of  these 
structures  is  still  lacking. 

Recently,  an  interesting  flame  instability  phenomena  was  observed  during  the  study  of  soot 
formation  in  high  pressure  diffusion  flames  which  possesses  many  of  the  desired  features  for  the  study 
of  this  phenomena.  The  flame  instability  was  observed  to  evolve  from  a  stable  methane/air  diffusion 
flame  as  the  pressure  was  raised  from  one  atmosphere  (14.7  psia)  to  approximately  6  psig  (20.7  psia). 
At  this  pressure,  the  flame  exhibited  a  regular  pulsation  or  oscillation.  In  order  to  examine  this 
phenomena  further,  a  series  of  high  speed  video  images  of  this  methane  flame  were  recorded.  These 
images  were  obtained  at  a  framing  rate  of  SOO  frames  per  second  and  simply  recorded  the  flame 
luminosity  from  the  soot  present  in  the  flame.  Analysis  of  these  video  recordings  revealed  the  onset 
and  evolution  of  a  very  regular  flame  structure  with  a  repetition  rate  of  approximately  15  Hz.  Several 
attempts  were  made  to  eliminate  these  pulsating  flames  by  modifying  the  pressure  chamber  with 
baffles  or  screens.  These  changes  were  intended  to  alter  the  acoustic  modes  of  the  chamber  or  to 
impede  flow  recirculation.  None  of  these  approaches  significantly  affected  the  phenomena.  Thus,  the 
observed  oscillations  are  not  believed  to  be  an  artifact  of  the  apparatus. 

Previous  researchers  have  observed  similar  self-excited  instabilities  in  laminar  fuel  jets  [49,50]. 

In  these  studies,  the  fuel  flow  rate  was  varied  to  produce  the  onset  of  the  instability.  Analysis 
explaining  the  oscillatory  behavior  emphasized  perturbations  introduced  into  the  velocity  field  which 
were  amplified  through  either  linear  or  non-linear  mechanisms  [49,50].  More  recent  work  has  placed 
a  stronger  emphasis  on  the  effect  of  buoyancy  in  both  laminar  and  turbulent  flow  regimes  [51,52-55]. 
The  general  trends  indicate  that  as  Reynolds  number  increases  (Re  >7000),  the  effect  of  these 
buoyancy  related  instabilities  is  less  important  [51,55].  The  present  evidence  is  not  conclusive  on  the 
potential  importance  of  such  low  frequency  instabilities  with  regard  to  turbulent  flows.  There  are 
studies  which  indicate  that  the  above  observed  instabilities  can  couple  to  the  inner  core  flow  of  the  jet 
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and  may  be  important  in  the  dynamic  processes  which  occur  during  transition  to  fully  turbulent 
conditions  [51,54].  Recent  visualization  studies  of  atmospheric  pressure  difhision  flames  have  shown 
that  these  flames  are  characterized  by  complex  vortical  structures  [54]. 

The  general  flame  structure  of  the  present  high  pressure  oscillating  flames  is  similar  to  that 
reported  in  the  studies  referenced  in  the  above  discussion.  The  observed  frequency  is  quite  close  to 
the  17  Hz  frequency  predicted  by  the  analysis  of  Buckmaster  and  Peters  [52].  Thus,  it  seems  clear 
that  these  instability  features  are  buoyancy  driven.  Recent  modelling  results  for  unsteady  buoyancy- 
driven  jet  diffusion  flames  conducted  at  NIST  have  been  quite  successful  in  reproducing  the  general 
characteristics  of  the  observed  flames  [46].  In  the  formulation  of  this  model,  three  dimensionless 
parameters  enter  into  the  formulation:  the  Reynolds  number,  the  Peclet  number  and  the  Richardson 
number.  In  the  present  experiments,  because  the  mass  flow  rates  of  friel  and  oxidizer  are  maintained 
at  a  constant  value,  the  product  of  the  density  and  inlet  gas  velocity  remains  constant  as  the  pressure 
is  changed.  Since  the  viscosity  is  independent  of  the  pressure,  the  Reynolds  number  also  remains 
constant  as  the  pressure  is  varied.  In  the  limit  of  unity  Lewis  number,  the  Peclet  number  is  equal  to 
the  fuel  jet  radius  multiplied  by  the  ratio  of  inlet  gas  velocity  and  the  binary  diffusion  coefflcient  of 
the  fuel.  Since  both  the  inlet  gas  velocity  and  the  diffusion  coefficient  are  inversely  proportional  to 
the  pressure,  the  Peclet  number  also  remains  constant.  Finally,  the  Richardson  number,  which 
expresses  the  ratio  of  buoyancy  forces  to  inertial  forces,  can  be  expressed  as  gLAJ^  where  g  is  the 
gravitational  constant,  L  is  the  fuel  jet  radius  and  U  is  taken  as  the  air  inlet  velocity.  Clearly,  this 
number  will  vary  with  the  square  of  the  pressure  through  the  velocity  term  in  this  ratio.  Thus,  the 
effective  gravitational  acceleration  will  vary  as  the  pressure  squared.  Consequently,  in  the  model,  a 
pressure  of  two  atmospheres  can  be  simply  simulated  by  choosing  a  gravitational  constant  of  4  g. 

The  significance  of  the  present  experimental  configuration  is  that  relative  effects  of  buoyancy  can 
be  examined  in  isolation  from  the  other  parameters.  This  situation  is  unique  to  this  burner 
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configuration  and  offers  an  excellent  opportunity  for  investigating  the  fundamental  aspects  of  the 
instability  onset  from  both  experimental  and  theoretical  modelling  aspects. 

Model  comparisons  with  the  imaging  results  described  above  show  good  agreement  with  respect 
to  the  instability  frequency  and  tip  cutting  location  [46].  This  agreement  includes  the  predication  of 
the  observation  that  the  tip  cutting  location  decreases  as  the  pressure  increases.  One  area  where 
a  ment  has  not  been  obtained  to  date,  involves  the  pressure  at  which  the  onset  of  the  instability  is 
observed.  The  model  predicts  a  pressure  of  about  one  half  that  observed  experimentally.  Appendix  3 
contains  a  detailed  discussion  of  the  model  and  experiment  comparisons. 

Additional  experiments  using  ethane  and  ethene  have  also  displayed  similar  flame  instability 
behavior  as  the  pressure  is  increased.  For  each  flame,  however,  the  onset  of  the  instability  occurs  at 
distinctly  different  pressures,  even  though  the  gas  flow  rates  are  identical.  Thus,  a  measurable  fuel 
dependence  is  observed  with  respect  to  the  onset  conditions.  Previous  studies  at  atmospheric 
conditions  have  not  investigated  fuel  related  effects,  because  the  instability  frequency  generally  is  not 
strongly  affected  [52].  For  the  high  pressure  initiated  instabilities,  however,  measurable  fuel 
constituent  effects  are  observed  and,  thus,  provide  a  potential  for  further  insights  into  the  instability 
phenomena. 
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4.0  CONCLUSIONS 


A  fundamental  study  of  soot  particle  inception  and  growth  processes  has  been  completed  which 
specifically  considered  the  effects  of  fuel  molecular  structure,  fuel  concentration,  temperature  and 
operating  pressure  on  the  formation  of  soot  particles.  Coannular  laminar  diffusion  flames  burning  in 
air  were  studied  for  pure  fuels  (methane  and  ethene)  as  well  as  fuel  mixtures  involving  ethene 
propene,  butane,  1-butene  and  1,3  butadiene  individually  mixed  with  methane  or  ethene.  The  fuel 
structure  studies  indicated  that  the  fuel  species  most  strongly  affected  the  particle  inception  process  for 
the  fuels  studied,  as  opposed  to  the  surface  growth  process.  Laser  fluorescence  measurements 
indicate  that  the  soot  particle  precursors  concentrations  increased  and  were  observed  sooner  in  the 
flame  as  the  sooting  propensity  of  the  fuel  increased.  Surface  growth  rate  constants,  when  normalized 
by  the  available  surface  area,  were  typically  within  a  factor  of  two  for  the  fuels  studied.  These 
findings  support  an  interpretation  in  which  inception  controls  the  maximum  amount  of  soot  formed  in 
the  flame. 

Studies  in  which  fuel  concentration  and  flame  temperature  were  varied  through  dilution  of  fuel 
with  an  inert  species  indicated  that  temperature  effects  dominate  changes  due  to  variations  in  the 
concentration  under  most  conditions.  Local  temperature  and  concentration  measurements  were  shown 
to  be  absolutely  necessary  to  properly  interpret  the  effects  of  inert  dilution.  In  particular,  small 
temperature  differences  in  the  initial  soot  formation  region,  as  well  as  the  mitigation  of  the  initial 
concentration  at  the  fuel  tube  exit  by  inter-diffiision  of  nitrogen  from  the  air  stream  and  fuel,  were 
shown  to  make  simple  interpretation  of  the  effects  of  dilution  problematic.  Based  on  studies  over  a 
wide  range  of  flow  conditions,  it  was  possible  to  estimated  the  separate  effects  of  dilution  and 
concentration.  This  estimate  showed  that  for  the  conditions  studied,  the  temperature  effects  dominate 
over  the  concentration  effects.  The  analysis  yielded  an  apparent  activation  energy  of  94.5  kcal/mole 
for  the  temperature  dependence.  The  fuel  concentration  dependence,  represented  as  where  Xq 
is  the  fuel  concentration  at  the  exit  of  the  fuel  tube,  was  best  represented  by  n  =  0.3. 


60. 


Studies  of  the  evolution  of  the  soot  particles  in  laminar  diffusion  flames  indicated  that  clusters  or 
aggregates  are  found  in  the  higher  parts  of  the  flame.  These  clusters,  which  are  composed  of  nearly 
spherical  primary  particles,  possess  different  optical  properties  from  their  spherical  equivalents. 
Through  a  combination  of  transmission  electron  microscope  and  laser  light  scattering  measurements, 
an  analysis  of  the  optical  properties  of  these  aggregates  was  developed. 

A  data  reduction  method  that  is  based  on  the  use  of  optical  cross  sections  for  aggregates  was 
applied  to  measurements  of  local  extinction  and  differential  scattering  at  three  angles  by  soot 
aggregates  in  a  laminar  ethene  diffusion  flame.  As  compared  to  the  Rayleigh  data  reduction  for 
spherical  particles,  this  method  yields  higher  values  of  volume-mean  diameter,  a  larger  surface  area 
per  unit  volume,  and  lower  values  of  the  aggregate  number  concentration.  The  analysis  also  provides 
the  mean-square  radius  of  gyration  and  the  average  number  of  primary  particles  per  aggregate.  Both 
these  quantities  increase  monotonically  with  time  as  is  consistent  with  the  occurrence  of  cluster-cluster 
aggregation.  Primary  particle  diameter  is  also  obtained  from  the  optical  observations  and  it  has  been 
independently  measured  by  thermophoretic  sampling  followed  by  TEM  analysis.  By  requiring  the 
rms  difference  between  the  two  observations  of  primary  particle  size  to  be  a  minimum,  information 
on  the  allowable  values  of  refractive  index  was  obtained.  Higher  values  of  the  real  portion  of  the 
refractive  index  give  more  plausible  moment  ratios  of  the  size  distribution  of  the  aggregates.  The 
present  results  indicate  that  a  self-consistent  interpretation  of  the  light  scattering  properties  of  the  soot 
aerosol  is  afforded  by  recognizing  its  aggregate  structure.  In  particular,  significant  differences  in  the 
surface  area  and  number  concentration  are  observed  for  soot  aerosols  consisting  of  aggregates  of 
primary  particles.  These  differences  have  important  ramifications  in  deriving  chemical  and  physical 
rate  data  from  the  laser  scattering  extinction  experiment.  Furthermore,  when  combined  with  TEM 
analysis  of  the  primary  particles  which  constitute  the  aggregates,  this  analysis  provides  the  basis  for 
evaluating  the  importance  of  other  particle  properties,  such  as  refractive  index.  Thus,  this  method  of 
data  analysis  provides  a  more  detailed  and  comprehensive  description  of  the  development  and 
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character  of  the  soot  aggregate  population  in  laminar  diffusion  flames. 

Operating  pressure  effects  on  soot  formation  were  also  studied.  These  studies  indicated  that  a 
power  law  representation  of  the  form  of  P“,  where  P  is  the  operating  pressure,  accurately  represented 
the  behavior  of  the  soot  volume  fraction  as  the  pressure  was  varied.  For  studies  involving  the  alkene 
species,  ethene  and  propene,  the  value  of  n  was  found  to  be  close  to  first  order  (l.OS  ±  0.06  and 
1.16  ±  0.01,  respectively).  Studies  of  an  alkane  fuel,  ethane,  resulted  in  a  significantly  larger  value 
for  the  power  dependence  closer  to  second  order  in  n  (1.62  -  1.88,  depending  on  the  flow  rate). 

These  results  confirm  the  strong  dependence  of  soot  formation  on  the  operating  pressure  and 
furthermore  point  to  potential  fuel  structure  effects. 

Finally,  studies  of  these  high  pressure  diffusion  flames  have  resulted  in  the  observation  of  the 
onset  of  buoyant  instabilities  induced  by  changes  in  the  pressure.  A  simple  procedure  has  been 
described  for  isolating  these  buoyancy  effects  in  jet  diffusion  flame  experiments.  All  that  is  necessary 
is  that  background  pressure  be  varied  while  maintaining  constant  mass  flows  of  fuel  and  oxidizer  into 
the  burner.  A  theoretical  model  for  these  flames  indicates  that  this  procedure  isolates  Richardson 
number  (or  relative  buoyancy  force)  as  the  only  variable  parameter,  which  is  equivalent  to  varying  the 
gravitational  acceleration.  A  joint  series  of  experiments  and  computations  involving  a  pressurized  low 
speed  methane/air  diffusion  flame  has  provided  strong  evidence  to  support  the  theory. 

It  is  clear  that  much  work  needs  to  be  done  in  order  to  clarify  the  important  role  that  buoyancy 
plays  in  jet  diffusion  flame  dynamics.  The  results  of  the  present  studies  only  open  the  door  to  future 
studies,  which  can  now  be  performed  with  reasonable  effort  and  cost-effectiveness.  Future  research  is 
also  needed  in  order  to  determine  the  limits  of  usefulness  of  the  buoyancy-isolation  procedure 
described  here.  Clearly,  use  of  this  procedure  to  attain  very  small  effective  gravitational  accelerations 
appears  unreasonable,  as  chemical  kinetics  could  certainly  not  be  ignored  at  the  low  pressures 
required.  Additionally,  enhanced  soot  production  at  low  g  would  not  be  simulated  by  this  pressure 
variation  technique.  However,  relaxation  of  the  requirement  that  Lewis  number  be  unity  appears 
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feasible,  since  Prandtl  and  Schmidt  numbers  can  both  be  considered  as  independent  of  pressure. 
Finally,  since  the  basic  model  described  here  is  not  dependent  on  a  particular  geometry,  other  types 
of  diffusion  flames  in  different  burner  configurations  could  also  be  considered  as  candidates  for 
buoyancy-isolation  experiments.  Thus,  the  possibilities  for  future  research  efforts  in  the  area  of 
buoyancy  effects  in  diffusion  flames  appear  rather  substantial. 
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Abstract 


A  detailed  study  has  been  made  of  the  relative  effects  of  concentration  and  temperature  on 
soot  formation  in  diluted  laminar  ethene/air  diffusion  flames  employing:  (1)  argon  and  nitrogen 
diluted  flames  at  equal  calculated  adiabatic  flame  temperatures  but  different  initial  concentrations  to 
isolate  the  effect  of  concentration  and  (2)  argon  and  nitrogen  diluted  flames  with  equal  initial  fuel 
concentrations  but  different  temperatures  to  isolate  the  effect  of  temperature.  Total  integrated  soot 
volume  fraction  measurements  show  that  in  argon  and  nitrogen  diluted  flames  with  equal  calculated 
adiabatic  flame  temperatures,  the  more  diluted  argon  flames  consistently  display  lower  soot 
concentrations.  However,  in  flames  of  equal  dilution,  argon  flames  consistently  display  higher  soot 
concentrations  than  the  slightly  cooler  nitrogen  diluted  flames.  Local  temperature  measurements 
show  that  for  nitrogen  and  argon  diluted  flames  with  equal  calculated  adiabatic  flame  temperatures, 
argon  diluted  flames  display  lower  temperatures  in  the  region  where  soot  is  first  formed.  Mass 
spectrometric  measurements  of  gas  species  concentrations  in  diluted  and  undiluted  methane/air 
diffusion  flames  were  obtained  and  compared  with  a  numerical  flame  model.  These  results  show  that 
the  initial  difference  in  fuel  concentration  in  diluted  and  undiluted  flames  diminishes  rapidly  with 
height.  Furthermore,  laser  light  scattering  measurements  show  that  as  inert  diluent  is  added  to  the 
flame,  soot  inception  is  delayed  and  consequently  less  time  is  available  for  soot  growth.  Based  on  this 
extensive  set  of  data,  a  quantitative  assessment  was  made  of  the  relative  effects  of  temperature  and 
concentration  indicating  that  for  coflow  diffusion  flames  the  role  of  temperature  is  more  important 
than  the  decrease  in  reaction  rate  due  to  the  reduction  in  fuel  concentration  when  an  inert  diluent  is 
added  to  the  fuel  flow. 

Introduction 

The  formation  of  soot  particles  in  combustion  environments  involves  a  complex  series  of 
chemical  and  physical  processes  which  result  in  the  conversion  of  carbon  contained  in  the  fuel  to 
carbonaceous  soot  particles.  A  variety  of  combustion  systems  have  been  previously  examined  to 
investigate  the  specific  details  of  the  soot  formation  process.  These  studies  range  from  simple 
laboratory  premixed  and  diffusion  flames  to  more  practical  scale  combustion  involving,  for  example, 
gas  turbines,  diesel  engines  and  furnaces.  Extensive  reviews  of  soot  formation  phenomena  exist 
summarizing  many  of  these  studies  [Glassman,  1988;  Haynes  and  Wagner,  1981;  Wagner,  1978, 
1981).  These  previous  studies  investigated  the  basic  importance  of  parameters  such  as  pressure, 
temperature,  and  fuel  molecular  structure.  An  additional  parameter  of  potential  importance  to  which 
recent  attention  has  been  given  is  fuel  concentration. 
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One  means  to  investigate  the  effect  of  fuel  concentration  is  through  the  introduction  of  an 
inert  diluent  into  the  fuel  stream.  Several  studies  note  that  the  addition  of  an  inert  diluent  to  the  fuel 
reduces  soot  concentration  in  laminar  diffusion  flames  [e.g..  Deardon  and  Long.  1968;  McLintock, 
1968).  This  reduction  is  in  part  a  result  of  lower  flame  temperatures  achieved  in  diluted  flames  as 
compared  to  the  undiluted  case.  These  lower  temperatures  reduce  the  pyrolysis  reaction  rates  leading 
to  decreased  soot  formation  rates  and  consequently  lower  soot  concentrations  [Classman  and 
Yaccarino.1981].  Reducing  the  fuel  concentration  may  also  have  an  effect  on  soot  concentration 
through  a  proportional  reduction  in  fuel  species  concentration  in  the  soot  formation  region,  i.e., 
reaction  rates  are  decreased  by  lower  species  concentrations.  The  objective  of  the  present  paper  is  to 
examine  the  effects  of  changes  in  temperature  and  fuel  concentration  on  soot  formation  when  an 
inert  diluent  gas  is  mixed  into  the  fuel  stream  of  a  coarmular  laminar  diffusion  flame.  This  study  also 
considers  other  processes  that  are  influenced  by  inert  dilution. 

Although  reduction  in  soot  through  inert  diluent  addition  has  been  observed  by  numerous 
researchers,  only  a  few  studies  exist  which  investigate  the  relative  roles  of  fuel  concentration  and 
flame  temperature.  Previous  studies  of  laminar  diffusion  flames  have  employed  dilution  as  a  means 
to  study  temperature  effects  on  soot  formation  [Classman  and  Yaccarino.  1981;  Santoro  and 
Semerjian.  1984].  In  these  cases  the  effects  of  decreasing  fuel  concentration  were  argued  to  be  small 
relative  to  the  temperature  reduction  contribution.  However.  Kent  observed  that  soot  concentrations 
in  diluted  flames  are  reduced  by  an  amount  that  can  not  be  explained  by  the  reduction  in  flame 
temperature  alone  [Kent  and  Wagner,  1984b].  Kent  and  Wagner  [1984b]  concluded  that  fuel 
concentration  appears  to  be  a  significant  factor  influencing  soot  production  rates.  Markstein  [1984, 
1986]  in  a  series  of  studies  examining  the  radiative  characteristics  of  laminar  diffusion  flames  also 
commented  on  the  relative  contributions  of  temperature  and  concentration  effects  when  diluents  are 
added  to  the  flame  [Markstein,  1986].  Axelbaum  and  et  al.  addressed  this  issue  quantitatively  in  a 
series  of  experiments  in  counterflow  diffusion  flames  [Axelbaum  et  al..  1988a,  1988b].  The 
important  conclusion  from  these  studies  is  that  soot  formation  rates  depend  linearly  on  the  initial  fuel 
concentration.  Axelbaum  and  Law  recently  extended  their  dilution  studies  to  include  coflow 
diffusion  flames  [Axelbaum  and  Law,  1990].  They  report  a  dilution  effect  for  coflow  flames  that  is 
consistent  with  the  counterflow  flame  results;  that  is,  the  maximum  soot  volume  fraction  increases 
linearly  with  the  fuel  concentration  at  the  burner  exit.  These  results  imply  that  when  diluents  arc 
added  to  the  fuel  the  role  of  temperature  in  soot  formation,  commonly  believed  to  be  the  governing 
parameter  of  soot  formation  rates,  is  in  some  situations  less  important  than  fuel  concentration.  Recent 
studies  in  similar  coflow  diffusion  flames  have  been  reported  in  which  the  temperature  field  is  varied 
through  preheating  of  the  fuel  and  air  streams  [Giilder  and  Snelling.  1990;  Giilder,  1992).  The 
results  of  these  studies  have  been  interpreted  to  support  the  dominant  role  of  dilution  effects  on  soot 
formation  in  these  laminar  diffusion  flame  experiments.  In  the  present  investigation,  the  'olc  of 
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dilution  on  soot  reduction  is  investigated  by  making  detailed  measurements  of  the  critical  quantities: 
soot  particle  concentration,  gas  species  concentrations  and  temperature. 

Experimental 

Investigating  the  role  of  temperature  and  concentration  on  soot  formation  in  laminar 
diffusion  flames  requires  careful  selection  of  experimental  conditions  and  measurement  approaches. 
Following  a  method  similar  to  the  approach  used  by  Axelbaum  and  Law  [1990],  dilution  is  studied 
through  a  series  of  flames  at  identical  calculated  adiabatic  flame  temperatures,  but  which  have 
different  initial  fuel  concentrations.  The  differences  in  heat  capacities  of  argon  and  nitrogen  in 
theory  permit  comparison  of  flames  with  identical  temperatures,  but  different  dilutions.  Therefore, 
the  effect  of  concentration  is  systematically  isolated  from  temperature.  Similarly,  flame  conditions 
can  be  selected  with  identical  dilutions.  e.g..{50%  N2  and  50%  Ar),  but  different  temperatures,  thus 
isolating  the  effect  of  temperature  on  soot  formation.  The  present  study  examines  both  flame  diluent 
conditions.  Desired  flame  temperature  conditions  are  calculated  with  the  NASA  Chemical 
Equilibrium  Code  [Gordon  and  McBride,  1976)  based  on  the  assumption  that  the  temperatures  in  the 
formation  region  of  laminar  diffusion  flames  scale  with  calculated  adiabatic  temperatures  [Boedeker 
and  Dobbs,  1986a]. 

A  coannular  burner  operating  at  atmospheric  pressure  conditions  was  used  to  study  laminar 
diffusion  flames  burning  ethene  or  methane.  A  laser  scattering/extinction  system  was  used  to  obtain 
data  on  the  soot  particle  field.  The  burner  and  light  scattering  apparatus  have  been  previously 
described  in  detail  [Santoro  et  al.,  1983]  and  are  only  briefly  reviewed  here.  The  coannular  burner 
consists  of  an  inner  brass  fuel  tube  (1.1  cm  id)  surrounded  by  an  outer  tube  (10.0  cm  id)  for  air  flow. 
A  40.5  cm  long  brass  cylinder  was  used  as  a  chimney  to  shield  the  flame  from  laboratory  air 
currents.  Slots  0.476  cm  high  and  2.54  cm  wide  were  machined  in  the  chimney  providing  access  for 
the  incident,  transmitted  and  scattered  light.  These  slots  also  allowed  access  for  the  thermocouple 
probe.  Two  of  the  slots  were  modified  to  provide  access  for  a  mass  spectrometer  sampling  probe 
which  required  a  1.252  cm  circular  opening  in  the  chimney.  The  fuel  or  fuel  inert  diluent  mixture  is 
burned  in  a  highly  over-ventilated  air  flow,  minimizing  the  effect  of  air  flow  rate  [Roper  et  al.,  1977a, 
1977b].  Fuel  flow  rates  are  monitored  using  calibrated  rotameters,  while  a  mass  flow  meter  measured 
the  air  flow  rates.  Fuel  gases  had  a  stated  purity  of  99.0%  and  99.5%  for  methane  and  ethene 
respectively.  Air  was  supplied  from  an  in-l.'-  se  compressor  which  was  filtered  to  remove  particles 
and  moisture.  Nitrogen  and  argon  were  added  to  the  fuel  flow  using  separately  calibrated  rotameters 
and  had  .stated  purities  of  99.99%  and  99.995%  respectively.  The  air  flow  rate  for  the  flames  studied 
was  1062  cm^/s.  Table  1  summarizes  the  ethene  flame  flow  conditions,  which  span  a  range  of  flow 
rates,  dilutions,  and  calculated  adiabatic  flame  temperatures. 
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The  laser  light  scattering/extinclion  apparatus  used  a  4W  argon  ion  laser  operating  at  the 
514.5  nxn  laser  line.  A  mechanical  chopper  operating  at  1  kHz  modulated  the  laser  source  allowing 
synchronous  detection  of  the  transmitted  light  signals.  The  laser  beam  was  focused  into  the  burner 
using  a  40  cm  focal  length  lens  and  was  typically  operated  with  an  output  power  of  0.5W.  The 
transmitted  light  was  detected  using  a  silicon  photodiode.  The  transmitted  laser  light  intensity  was 
reduced  by  a  neutral  density  filter  (N.D.=2)  to  a  level  suitable  for  linear  photodiode  response.  The 
output  of  this  detector  was  input  into  a  two  phase  lock-in  amplifier  which  was  interfaced  to  an  IBM- 
XT  computer. 

Laser  light  scattering  measurements  at  a  scattering  angle  of  90°  with  respect  to  the  incident 
beam  could  also  be  obtained  with  this  system.  However,  in  the  present  experiments  such 
measurements  were  only  used  for  comparisons  with  previous  smdies  and  to  determine  the  location 
where  soot  particles  first  form. 

The  entire  atmospheric  burner  was  mounted  on  a  three-dimensional  translating  stage  system. 
Computer  controlled  stepper  motors  were  used  to  adjust  the  vertical  and  one  of  the  horizontal 
coordinates.  This  allowed  radial  profiles  of  the  laser  extinction  to  be  achieved  at  various  axial 
positions  in  the  flame. 

Rapid  insertion  thermocouple  measurements  provided  measurements  of  the  temperature 
distributions  in  the  flames  studied.  The  rapid  insertion  thermocouple  technique  [Kent  and  Wagner. 
1984a]  avoids  continuous  coating  with  soot  by  burning  off  deposited  soot  on  the  thermocouple 
surface  by  positioning  the  thermocouple  in  the  oxidation  zone  of  the  flame  prior  to  making 
measurements  at  a  new  location.  Pt/Pt-10%  Rh  fine  wire  thermocouples  were  made  from  127  pm 
diameter  wire.  The  resulting  bead  diameter  was  160  pm.  A  limited  number  of  temperature  profiles 
were  also  taken  with  a  thermocouple  whose  Junction  was  made  with  63.5  pm  wire  (95  pm  bead 
diameter)  to  verify  some  of  the  trends.  The  uncoated  thermocouples  were  mounted  on  a  stepper 
motor  stage  for  rapid  positioning.  This  technique  provides  accurate  measurement  of  relative  changes 
in  the  temperature  as  diluent  spiecies  and  concentrations  were  varied. 

In  order  to  examine  the  effect  of  diluent  addition  on  the  local  gas  species  concentrations  in 
the  flame,  a  series  of  mass  spectrometric  measurements  were  obtained.  These  measurements  are 
similar  in  approach  to  previous  studies  [Smyth  et  al..  1985).  Coflowing  methane/air  flames  were 
studied  in  order  to  avoid  the  difficulties  introduced  by  the  presence  of  large  soot  concentrations 
which  clog  the  quartz  microprobe  orifice  used  in  the  sampling  procedure.  The  quartz  microprobe 
used  in  these  studies  had  a  90  pm  orifice.  This  probe  was  connected  to  an  Extrel  model  C50 
quadropole  mass  spectrometer.  Species  mole  fraction  information  from  the  measurements  was 
obtained  using  a  procedure  in  which  the  spectrometer  sensitivity  to  a  particular  species  is  measured 
using  a  calibration  mixture  at  room  temperature.  The  relative  sensitivities  are  assumed  to  be 
independent  of  temperature  allowing  calculation  of  the  actual  mole  fraction  by  assuming  the  mole 
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fraction  of  measured  species  sums  to  unity  at  each  measurement  location  (Bittner.  1981],  An 
undiluted  methane  flame  and  a  50%  diluted  methane  (with  nitrogen)  flame  were  studied.  The 
methane  fuel  flow  rate  was  5.70  cm^/s  for  both  flames  with  an  equal  amount  of  nitrogen  added  for 
the  diluted  case.  The  air  flow  rate  for  these  studies  was  1062  cm^/s. 

Analysis 


In  the  present  study,  the  quantity  of  interest  is  the  change  in  soot  volume  fraction,  resulting 
from  dilution  of  the  fuel.  In  previous  studies  of  laminar  diffusion  flames,  tomographic 
reconstruction  approaches  have  proven  useful  in  providing  spatially  resolved  measurements  of 
from  line-of-sight  extinction  measurements  [Santoro  et  al..  1983].  In  the  present  case,  however, 
errors  associated  with  the  tomographic  reconstruction  approach  detract  from  its  usefulness  since  small 
relative  variations  in  need  to  be  measured. 

An  alternate  approach  is  to  relate  the  extinction  measurement  to  an  appropriate  spatial 
integral  of  the  soot  volume  fraction.  For  soot  particles  in  the  Rayleigh  size  limit  (d«A  )  the 
extinction  (///<,)  is  defined  as  the  ratio  of  the  transmitted  laser  intensity  (I)  to  incident  laser  intensity 
(lo) .  and  is  related  to  as: 


-c(A,m)ln 


\hj 


Eq-l’ 


where  x  is  the  direction  along  which  the  laser  beam  propagates,  c(A,m)  is  a  constant  determined  from 
Rayleigh  theory  which  depends  on  the  wavelength.  A,  of  the  laser  and  the  refractive  index  ,  m,  of  soot 
particles  and  x  is  the  path  length  through  the  flame.  As  in  previous  work,  the  value  of  complex  index 
of  refraction  used  is  m=1.57-0.56i  [Dalzell  and  Sarofim,  1969].  The  constant,  c(A,m)  calculated 
from  this  value  is  1.05  x  lO"^  cm.  A  measure  of  the  total  integrated  soot  volume  fraction,  F^,  can  be 
obtained  by  integrating  equation  (1)  along  the  direction  perpendicular  to  x.  This  integrated  volume 
fraction  is  a  measure  of  the  total  amount  of  soot  at  a  particular  axial  location  and  is  given  by 


F",  =  -  j  c(A,/7j)ln 


\^oJ 


\dy 


Eq-2. 


where  y  is  the  direction  perpendicular  to  x.  The  total  integrated  soot  volume  fraction,  has  been 
selected  because  it  also  incorporates  any  contributions  due  to  changes  in  flame  diameter  as  the 
fuel/diluent  flow  rate  is  varied. 
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Results 


The  ensuing  sections  present  and  discuss  the  experimental  results  from  the  dilution  studies. 
These  studies  address  four  aspects  of  the  influence  of  inert  additives  on  the  reduction  in  soot 
formation;  1)  the  relation  between  fuel  species  concentration  and  soot  concentration;  2)  the  relation 
between  temperature  and  soot  formation;  3)  the  evolution  of  species  concentrations  in  both  diluted 
and  undiluted  flames;  and  4)  variation  in  soot  particle  residence  time  associated  with  dilution.  Using 
the  information  from  these  experiments,  the  relative  influence  of  fuel  concentration  and  temperature 
is  then  considered  for  coflowing  laminar  diffusion  flames. 

Constant  Adiabatic  Flame  Temperature  Studies 
Temperature  Results 

The  approach  used  to  determine  the  relationship  between  the  degree  of  dilution  (i.e.,  fuel 
concentration)  on  soot  formation  involved  comparing  argon  and  nitrogen  diluted  flames  having  the 
same  calculated  adiabatic  flame  temperature.  Since  argon  and  nitrogen  have  different  molar  heat 
capacities,  the  initial  fuel  concentration  of  these  flames  must  be  different  if  the  temperature  is  to  be 
maintained  constant.  To  ensure  that  the  observed  differences  in  soot  formation  rates  in  these  flames 
of  equal  calculated  adiabatic  flame  temperature  are  a  result  of  differences  in  fuel  concentration  and 
not  temperature,  the  temperature  field  was  determined  from  rapid  insertion  thermocouple 
measurements.  In  addition,  these  measurements  allow  a  determination  of  whether  calculated  adiabatic 
flame  temperatures  may  be  used  to  represent  the  changes  in  temperature  expected  with  dilution.  This 
may  not  be  self  evident  since  calculated  adiabatic  flame  temperatures  are  based  on  the  fuel/diluent 
mixture  being  supplied  to  the  burner  and  not  the  local  fuel  concentration  of  the  flames. 

Thermocouple  measurements  of  the  temperature  field  were  obtained  in  the  4.90  cm^/s  ethene 
flames  (see  Table  1)  whose  calculated  adiabatic  flame  temperatures  were  2369K  (flame  12),  2346K 
(flames  5  and  6)  and  2333K  (flames  7  and  8).  The  two  lower  temperature  conditions  were 
established  individually  for  nitrogen  (flames  5  and  7)  and  argon  (flames  6  and  8)  diluent.  Radial 
temperature  profiles  at  several  axial  locations  were  obtained  for  flames  12  and  7  while  the  other 
flames  were  examined  only  at  an  axial  location  5  mm  above  the  fuel  tube  exit.  Uncorrected 
measured  temperatures  at  four  selected  axial  locations  are  shown  in  Figure  1  for  an  undiluted  and 
diluted  flame  (flames  12  and  7).  The  radial  temperature  profiles  of  flames  12  (undiluted)  and  8 
(diluted),  shown  in  Figure  1  as  a  function  of  axial  position,  reveal  a  trend  expected  in  the  temperature 
field.  For  the  displayed  heights  less  than  50  mm  from  the  burner  exit,  temperatures  are  reduced  in 
the  diluted  flame.  Above  50  mm  the  ordering  of  the  flames  reverses  and  the  diluted  flame  has  a 
higher  temperature.  This  reversal  in  temperature  is  due  to  more  soot  being  formed  in  the  initially 
hotter  undiluted  flame,  increasing  radiative  losses,  thus  cooling  the  flame  to  a  greater  extent  at  higher 
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locations  (Kent  and  Wagner,  1984a.  1984b;  Santoro  and  Semerjian.  1984;  Boedeker  and  Dobbs. 
1986a,  1986bl. 

Careful  inspection  of  the  temperature  profiles  (Figure  1)  reveals  that  low  in  the  flame,  dilution 
reduces  temperatures  at  the  centerline  to  a  greater  extent  than  in  the  annular  region  [Boedeker  and 
Dobbs,  1986a.  1986bl.  This  difference  in  the  effect  of  the  diluent  as  a  function  of  position  is  related 
to  the  diffusional  nature  of  these  flames.  For  both  the  nitrogen  and  argon  diluent  cases  as  well  as  the 
undiluted  flame  cases,  nitrogen  from  the  air  stream  is  diffusing  towards  the  centerline  while  the  "fuel" 
(diluted  or  undiluted)  diffuses  towards  the  flame  front.  This  interdiffusional  aspect  of  the  flame 
tends  to  mitigate  the  initial  fuel  concentration  differences  resulting  from  the  addition  of  the  diluent 
species.  This  point  will  be  discussed  in  detail  in  a  subsequent  section. 

A  comparison  of  the  effect  of  diluent  addition  on  the  temperature  profiles  is  further 
illustrated  in  Figure  2.  For  these  measurements  an  axial  location  5  mm  above  the  burner  exit  has 
been  selected.  This  is  near  the  location  where  soot  is  first  observed  to  form  in  these  flames  based  on 
light  scattering  profile  measurements.  Since  soot  particles  were  not  observed  or  exist  in  very  low 
concentration  at  this  location  [Santoro  and  Semerjian,  1984],  radiation  effect:  due  to  the  presence  of 
soot  particles  are  absent  and  do  not  complicate  the  interpretation  of  aicasured  temperature 
differences  for  various  flame  conditions.  Furthermore,  by  selecting  a  location  close  to  ue  fuel  tube 
exit,  potential  concentration  effects  introduced  by  diluting  the  fuel  flow  are  more  likely  to  be 
observed. 

From  the  temperature  profiles  it  is  observed  that  in  the  region  near  the  maximum 
temperature,  low  in  the  flame,  the  measured  temperature  differences  between  the  undiluted  and 
diluted  flames  are  similar  to  the  differences  in  calculated  adiabatic  flame  temperature  (see  Figure  2, 
insert).  Similar  results  using  CARS  temperature  measurements  were  obtained  by  Boedeker  and 
Dobbs  [1986a,  1986b].  Furthermore,  temperatures  observed  in  flames  in  which  diluent  conditions 
are  chosen  to  result  in  the  same  calculated  adiabatic  flame  temperature  differ  by  less  than  lOK  in  the 
region  near  the  temperature  maximum.  However,  as  the  centerline  is  approached  from  the 
temperature  maximum,  differences  between  the  Nt  and  Ar  diluted  flames  for  identical  calculated 
adiabatic  flame  conditions  increase  with  the  Nt  diluted  flames  having  slightly  higher  temperatures. 
Near  the  radial  location  (r=5  mm)  where  soot  is  formed  [Santoro  and  Semerjian,  1984;  Santoro  et  al., 
1983],  this  difference  can  be  as  large  as  33-56K  (see  Table  2).  Positioning  uncertainty  with  respect 
to  the  thermocouple  measurement  location  is  not  believed  to  be  the  source  for  these  observed 
differences,  since  both  the  centerline  minimum  and  the  symmetric  temperature  maximums  in  the 
annular  region  show  good  positional  agreement  for  each  of  the  flames  (see  Figure  2).  These 
observed  differences  are  likely  due  to  the  differences  in  the  thermal  transport  properties  in  the 
individual  flames.  Since  the  thermocouple  measurements  are  also  influenced  by  local  thermal 
transport  phenomena,  differences  in  the  observed  temperatures  could  be  a  result  of  this  systematic 
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variation.  To  investigate  this  possibility,  an  estimate  of  the  effect  of  the  difference  in  the  transport 
properties  between  argon  and  nitrogen  on  the  thermocouple  measurements  was  made  [Bradley  and 
Mathews,  1968].  This  estimate  indicated  that  at  most  one  third  of  the  measured  temperature 
difference  was  due  to  variation  in  transport  properties.  Furthermore,  the  thermocouple  measurements 
were  repeated  for  the  5  mm  axial  location  using  a  smaller  bead  (95  pm)  thermocouple.  These 
measurement  showed  very  similar  temperature  differences  as  compared  with  the  larger  bead 
thermocouple  (160  pm)  measurements.  Since  variations  in  the  bead  diameter  result  in  significant 
differences  in  the  heat  transfer  and  radiative  processes  involving  the  thermocouple  bead,  the 
insensitivity  of  the  measured  temperature  aifferences  to  variations  in  the  bead  size  suggests  that  the 
measured  differences  are  indeed  representative  of  actual  temperature  differences  in  the  flame. 

The  temperature  measurement  results  comparing  the  argon  and  nitrogen  diluted  cases  at  the  5 
mm  axial  location  for  several  flames  are  summarized  in  Table  2.  Listed  in  this  table  are  the 
calculated  adiabatic  flame  temperature  (Tad),  the  maximum  temperature  observed  (Tmax. 
uncorrected  for  radiative  effects),  and  the  difference  in  the  temperature  (ATr=5  „,„,)  observed  at  the  5 
mm  radial  position  between  the  Ni  diluted  and  Ar  diluted  flames  which  is  near  where  soot  is  first 
observed. 

The  following  observations  can  be  made  based  on  the  results  shown  in  Table  2.  Flames  with 
similar  calculated  adiabatic  flame  temperatures  achieve  similar  measured  maximum  temperatures. 
For  the  nitrogen  and  argon  flames  with  similar  calculated  temperatures,  the  observed  temperature 
difference  between  flames  increases  with  dilution  in  the  interior  region  of  the  flame  (i.e.  r=  5  mm  and 
AT).  Thus,  although  dilution  with  differing  diluents  can  be  used  to  achieve  similar  maximum 
temperature  conditions,  measurable  differences  in  the  temperature  field  will  result  in  regions  interior 
to  the  maximum  temperature  region.  These  observations  underscore  the  need  to  consider  the  local 
temperature  and  concentration  as  weU  as  the  overall  expected  behavior  of  the  flame  in  assessing  the 
effects  of  dilution. 

Soot  Formation  Results 

To  investigate  the  variation  in  soot  formation  rate  with  the  initial  fuel  concentration,  the  soot 
yield  in  flames  of  similar  temperatures  but  different  dilutions  were  undertaken  using  three  pairs  of 
flames,  each  pair  at  a  different  temperature.  The  variation  in  the  total  integrated  soot  volume  fraction, 
F^,  at  various  axial  locations  is  shown  in  Figure  3  for  each  of  the  ethene  flames  with  a  flow  rate  of 
4.90  cm^/s.  It  is  interesting  to  note  that,  although  addition  of  diluent  increases  the  burner  exit 
velocities,  the  height  at  which  the  peak  volume  fraction  is  reached  is  nearly  constant  in  these  flames 
with  possibly  the  exception  of  flame  4  which  has  the  highest  dilution.  The  height  at  which  the  soot 
maximum  is  reached  is  determined  by  the  location  where  the  oxidizer  replaces  the  fuel  at  the 
centerline  [Roper.  1977a,  1977bj.  Thus,  the  rate  of  diffusion  of  oxidizer  into  the  flame  front  is  only 
weakly  dependent  on  the  amount  of  fuel  dilution. 


8 


82, 


The  maximum  F^.  is  used  as  a  readily  determined  measure  of  the  sooting  propensity  of  these 
flames.  The  reduction  in  soot  formed  is  calculated  as  the  ratio  of  the  maximum  F^,  in  the  diluted  and 
undiluted  flames  (see  Table  1).  in  this  comparison,  the  reduction  in  soot  may  be  due  to  temperature 
and/or  concentration  effects.  As  expected,  the  maximum  soot  concentration  ratios  show  that  reducing 
fuel  concentration  and  temjjerature  reduces  the  amount  of  soot  formed.  The  flames  which  show  the 
largest  relative  reduction  in  soot  (flames  3  and  4).  the  most  diluted  flames,  show  a  dependence  on  the 
soot  formation  rate  which  is  directly  proponional  to  the  initial  fuel  concentration.  This  implies  soot 
formation  rates  vary  with  initial  fuel  concentration  in  less  than  a  first  order  dependence,  since  a 
temperature  reduction  effect  on  soot  formation  has  not  been  taken  into  account. 

Of  course  the  motivation  for  establishing  flames  with  the  same  temperature  field  but  different 
fuel  concentration  is  to  presumably  isolate  concentration  effects  from  temperature  effects.  The 
following  discussion  therefore  emphasizes  comparisons  between  diluted  flames  having  identical 
calculated  adiabatic  flame  temperatures.  Under  these  conditions  the  peak  soot  concentrations  are 
observed  to  decrease  to  a  greater  extent  for  the  more  heavily  diluted  argon  flames.  This  result  is 
similar  to  previous  observations  and  is  argued  to  be  a  measure  of  the  dilution  effect  [Axelbaum  and 
Law,  1990].  In  order  to  quantify  the  effects  of  fuel  concentration  on  soot  formation  rates,  the 
maximum  F^  was  taken  to  be  related  to  the  initial  mole  fraction  by: 

Eq-3. 

Based  on  the  ratio  of  the  maximum  values  and  the  known  initial  concentrations,  a  value  for  b  can 
be  obtained  for  flames  with  equal  calculated  adiabatic  flame  temperatures.  The  results  of  these 
calculations  are  shown  in  Table  3  for  several  flames.  As  can  be  seen  no  unique  value  of  b  results  and 
a  first  order  dependence  does  not  appear  to  be  generally  valid.  It  should  be  mentioned  that  relatively 
small  changes  in  F^  are  involved  in  some  of  the  results  and  thus  the  accuracy  in  b  may  not  be  high. 

The  above  results  indicate  that  the  concentration  effect  of  soot  formation  may  not  be  first 
order  in  the  initial  fuel  mole  fraction.  Furthermore,  for  flames  with  the  different  diluents  but  equal 
calculated  adiabatic  flame  temperature,  measurements  indicate  that  differences  in  temperature  exist  in 
the  region  in  which  soot  inception  occurs.  The  potential  significance  of  these  temperature 
differences  will  be  addressed  in  the  discussion  section. 

Constant  Concentration  Studies 

As  a  complement  to  the  studies  at  constant  temperature,  a  series  of  flames  was  also  studied  in 
which  the  initial  concentration  of  the  fuel  was  held  constant  and  the  temperature  was  varied.  In  this 
approach  the  variation  of  soot  concentration  due  to  temperature  differences  in  pairs  of  ethcne  flames 
equally  diluted  with  nitrogen  or  argon  were  measured.  For  these  constant  concentration  studies,  the 
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argon  flames  have  higher  calculated  adiabatic  flame  temperatures  (see  Table  1).  Since  argon  and 
nitrogen  diffuse  in  a  similar  manner,  these  flames  should  have  similar  concentration  fields.  Three 
pairs  of  flames  (Flames  1  and  la.  Flames  3  and  3a,  and  Flames  6  and  7)  were  compared  with  fuel 
mole  fractions  of  0.5,  0.5  and  0.64  respectively. 

Temperature  profile  measurements  at  an  axial  location  5  mm  above  the  fuel  tube  exit  show 
the  maximum  temperatures  in  the  flames  differed  by  16-25K  in  agreement  with  adiabatic  flame 
temperature  calculations.  In  addition,  centerline  temperatures  were  found  to  be  10-22K  larger  for  the 
argon  diluted  flames.  Any  meaningful  assessment  of  the  temperature  differences  between  these  two 
locations  was  not  possible  since  the  argon  profiles  were  observed  to  extend  to  a  slightly  larger  radial 
coordinate. 

Measurements  of  the  soot  concentrations  of  these  flames  with  the  same  dilution  show  that  the 
nitrogen  flames  with  lower  calculated  adiabatic  flame  temperature  produce  less  soot  (Figure  4).  Thus, 
reducing  temperature  clearly  reduces  the  production  of  soot.  In  conjunction  with  the  constant 
temperature  studies,  these  results  can  be  used  to  assess  the  relative  contributions  of  temperature  and 
concentration  on  the  soot  formation  process  in  these  diluted  flames. 

Concentration  Measurement  Results 

In  any  chemically  reacting  system  the  reaction  rates  must  clearly  be  sensitive  to  the 
concentration  of  the  reactants.  In  studies  where  the  initial  fuel  concentration  is  varied  through 
dilution,  the  occurrence  of  concentration  effects  on  the  soot  formation  process  appears  to  be  self 
evident.  However,  such  concentration  effects  should  be  related  to  the  local  concentration  in  the 
region  where  soot  formation  and  growth  occur.  Thus,  it  is  important  to  understand  the  marmer  in 
which  variations  in  the  initial  fuel  concentration  affect  species  concentration  at  other  locations  in  the 
flame.  For  this  reason  a  series  of  mass  spectrometric  profiles  were  obtained  in  an  undiluted  methane 
flame  and  a  methane  flame  in  which  the  fuel  was  diluted  to  50%  with  nitrogen.  The  objective  of 
these  studies  was  to  quantify  specifically  the  variation  in  local  species  concentration  as  the  initial  fuel 
concentration  was  changed.  Methane/air  flames  were  studied  in  order  to  avoid  probe  plugging 
problems  associated  with  regions  of  high  soot  concentration  typical  of  ethene  diffusion  flames. 

Radial  mole  fraction  profiles  of  major  species  (CH4,  N2,  O2,  CO2,  CO,  H2O.  and  H2)  were 
obtained  at  axial  positions  of  5  and  12  mm.  The  results  for  CH4,  N2,  O2,  CO2  at  12  mm  were 
compared  with  gas  chromatography  analysis  and  showed  agreement  to  better  than  10%.  Figures  5 
and  6  show  the  radial  concentration  profile  results  for  both  flame  conditions  at  an  axial  location  of 
12  mm.  As  is  clear  from  these  results,  the  initial  difference  in  the  fuel  concentration  rapidly  narrows 
as  diffusion  of  fuel  and  nitrogen,  as  well  as  combustion  products,  proceeds.  Comparisons  of  the 
measured  concentrations  for  both  flames  have  been  made  at  the  centerline  and  at  a  radial  position  of 
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5  mm.  The  5  mm  radial  position  is  again  taken  to  be  characteristic  of  the  soot  precursor  region.  For 
the  5  mm  axial  location,  the  ratio  of  the  measured  concentrations  in  the  diluted  and  undiluted 
methane  flames  are  0.71  and  0.69  for  the  5  mm  radial  location  and  centerline  (r  =0.0  mm) 
respectively.  At  the  12  mm  axial  location,  the  values  for  the  same  two  radial  locations  are  0.96  and 
U.92  respectively.  The  values  depart  markedly  from  0.5,  the  ratio  which  represents  the  initial 
concentration  ratio  at  the  fuel  tube  exit  for  the  diluted  and  undiluted  flame  cases.  The  rapid 
interdiffusion  of  nitrogen  from  the  air  stream  and  fuel  has  previously  been  noted  by  Boedeker  and 
Dobbs  (1986b].  It  should  also  be  mentioned  that  the  temperature  as  well  as  the  fuel  concentration  is 
varied  in  these  two  methane/air  diffusion  flame  cases.  However  the  temperature  differences  are  not 
large  enough  to  significantly  affect  the  diffusion  rates.  A  clear  observation  from  these  results  is  that 
variations  in  the  local  fuel  concentration  introduced  through  dilution  dlTer  significantly  from  those 
at  the  fuel  tube  exit. 

Residence  Time  Results 

Residence  time  effects  have  been  shown  to  be  important  in  diffusion  flame  studies  of  soot 
formation  [Santoro  et  al.,  1987;  Honnery  and  Kent,  1990).  In  fact,  the  mass  flow  rate  of  soot  in 
laminar  diffusion  flames  has  been  observed  to  have  a  squared  or  cubed  power  dependence  on 
residence  time  (Honnery  and  Kent.  1990].  Thus,  changes  in  velocity  which  result  from  addition  of 
diluent  species  could  affect  the  residence  time  and  the  ultimate  soot  volume  fraction  achieved  as  a 
result.  However  the  present  measurements  in  these  buoyancy  dominated  diffusion  flames  indicate 
that  the  major  observed  effect  of  dilution  with  respect  to  residence  time  is  a  shift  of  the  location  where 
soot  is  formed  to  higher  axial  positions.  Since  the  location  of  the  maximum  soot  volume  fraction  is 
largely  unaffected  by  dilution  (see  Figure  3).  the  total  soot  residence  time  is  reduced.  A  reduction  in 
residence  time  under  diluted  flame  conditions  may  account  for  a  significant  reduction  in  soot  volume 
fraction  and  should  be  accounted  for  separately  in  studies  where  the  effects  of  temperature  and 
concentration  are  studied.  To  determine  the  location  where  soot  particles  first  form,  the  scattered 
light  intensity  at  a  scattering  angle  of  90°  was  measured  in  the  lower  region  of  flames  1,  2.  9.  10,  11. 
and  13  of  this  study.  Additional  light  scattering  measurements  were  taken  for  3.85  cm^/s  C2H4  flame 
diluted  with  8.82  cm^/s  of  N2  (experiment  no.  7  of  Santoro  and  Semerjian  (1984]).  These  results, 
which  are  shown  in  Figure  7  indicate  that  once  soot  particles  are  present  the  maximum  scattered  light 
intensity  at  a  given  axial  location  increases  rapidly  with  height  above  the  burner.  Furthermore,  this 
figure  clearly  shows  that  as  the  flames  are  diluted  with  inert  gas,  inception  is  delayed,  and  less  time  for 
soot  growth  is  possible.  Note  that  the  shift  in  the  location  of  the  initial  soot  formation  region  could 
be  a  result  of  temperature,  velocity  and/or  concentration  effects. 
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Discussion 


From  the  above  results,  it  is  clear  that  studies  which  utilize  dilution  of  the  fuel  flow  to  var>' 
temperature  and  concentration  parameters,  are  difficult  to  unambiguously  interpret  in  terms  of  global 
effects.  Mixing  and  transport  processes  have  been  shown  in  this  work  to  locally  alter  the  global 
dilution  effects  of  temperature  and  concentration.  In  order  to  understand  the  relative  effects  of 
temperature  and  concentration  variations  on  soot  formation  processes,  it  is  important  to  account  for 
these  local  effects. 

There  are  two  major  observations  based  on  the  present  work  regarding  local  temperature  and 
concentration  effects  when  diluents  are  added  to  the  fuel  flow  in  laminar  diffusion  flames.  The  first 
involves  the  case  when  different  diluents  (e.g.,  N2  and  Ar)  are  used  to  achieve  constant  temperature 
conditions,  but  varying  concentration  conditions  in  the  flame.  The  present  work  indicates  that 
although  nearly  equal  maximum  temperatures  are  achieved,  interior  to  the  flame  front  in  the  region 
where  soot  particles  are  initially  formed,  temperature  differences  of  33-56  K  are  observed.  The 
second  observation  is  that  mixing  effects  due  to  interdiffusion  of  species,  such  as  nitrogen  from  the 
air  flow  and  combustion  products  from  the  reaction  zone,  mitigate  the  effects  of  dilution  rapidly  in 
these  flames.  Thus,  the  ratio  of  the  local  concentration  between  flames  of  varying  fuel  dilution  can 
differ  significantly  from  initial  concentration  ratio  at  the  fuel  tube  exit.  These  observations  imply  a 
more  complicated  interaction  between  temperature  and  concentration  effects  when  fuel  dilution  is 
used  to  study  soot  formation  than  was  anticipated  in  some  previous  work  [e.g.  Santoro  and  Semerjian, 
1984], 

In  the  discussion  below,  supporting  evidence  for  the  importance  of  the  observations  described 
above  will  be  given.  This  will  be  followed  by  an  analysis  which  is  intended  to  estimate  the  relative 
importance  of  the  contributions  of  temperature  and  concentration  effects  on  soot  formation  when 
fuel  dilution  is  used  as  a  means  to  alter  the  soot  formation  characteristics  of  the  flame.  The  later  topic 
is  of  interest  because  of  the  present  diversity  of  results  in  the  literature  regarding  this  matter 
[Classman  and  Yaccarino,  1981;  Santoro  and  Semeijian,  1984;  Axelbaum  and  Law,  1990;  Giilder 
and  Snelling,  1990;  Giilder,  1992], 

It  has  been  previously  mentioned  that  the  present  studies  indicate  that  small  temperature 
differences  are  observed  near  the  soot  formation  region  for  diluted  flame  conditions  of  equal 
calculated  adiabatic  flame  temperatures.  It  is  appropriate,  therefore,  to  inquire  as  to  the  potential 
significance  of  such  small  temperature  variations  (AT=33-56  K)  on  the  soot  fonnation  process. 

Results  from  previous  work  by  Kent  and  Wagner  [1984b]  for  laminar  diffusion  flames  and 
BOhm  et  al.  [1988]  for  premixed  flames  indicate  that  relatively  small  changes  in  temperature  have  a 
significant  affect  on  the  amount  of  soot  produced.  Kent  and  Wagner  [1984b]  investigated  the  effect 
of  fuel  tube  lip  conditions  on  the  production  of  soot  by  replacing  their  water  cooled  fuel  tip  with  an 
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uncooled  glass  tip  while  holding  all  other  conditions  constant.  They  found  for  fuel  flow  rates 
comparable  to  the  present  study  that  the  maximum  soot  volume  fraction  increased  by  approximately 
20%  for  the  glass  fuel  tip.  Accompanying  temperature  measurements  made  5  mm  above  the  lip  of 
the  fuel  tube  showed  very  small  differences  for  the  two  fuel  tube  configurations  [Kent  and  Wagner, 
1984b].  These  results  indicate  that  the  soot  particle  formation  process  is  quite  sensitive  to 
temperature  changes  in  the  region  of  the  flame  near  the  fuel  tube  exit. 

Studies  by  Bdhm  et  al.  [1988]  in  premixed  flames  show  for  ethene/air  flames  that  changes  in 
the  measured  temperature  of  50  K  can  result  in  approximately  a  factor  of  two  change  in  the  soot 
volume  fraction,  /yoo,  for  temperatures  between  1550  K  and  1625  K.  Further  increasing  the 
temperature  results  in /yoo  reaching  a  maximum  at  about  1675  K  and  then  decreasing  for  higher 
temperatures.  These  same  trends  were  observed  for  a  range  of  C/0  ratios  between  0.6  and  0.8  which 
correspond  to  an  equivalence  ratio  between  1.8  and  2.4  for  ethene  burning  in  air.  This  range  of 
temperatures  is  similar  to  the  range  for  which  soot  particles  are  observed  to  form  in  laminar  diffusion 
flames.  An  interesting  result  of  these  experiments  is  that  the  observed  temperature  sensitivity  of /yoo 
varies,  being  larger  for  temperatures  below  1600  K  than  for  temperatures  between  1600  K  and  1675 
K.  Above  1675  K,/yoo  decreases  as  temperature  is  further  increased.  Thus,  the  Bdhm  et  al.  [1988] 
results  clearly  show  that  the  temperature  sensitivity  observed  for  the  soot  formation  process  will 
depend  on  the  temperature  region  studied.  In  the  present  studies  we  argue  that  the  temperature  lies  in 
the  range  between  1350  K  and  1625  K  (see  Figure  2.  r=5  mm)  where  significant  temperature 
sensitivity  is  observed  [Bdhm  et  al.,  1988].  If  temperatures  were  raised  above  1650  K  by  some 
means,  one  would  expect  a  less  significant  temperature  sensitivity  or  possibly  a  decrease  in  the 
amount  of  soot  formed  if  temperatures  above  1675  are  achieved.  It  should  also  be  noted  that  the 
above  observations  are  for  ethene  flames.  Other  fuels  are  observed  to  have  different  temperature 
sensitivities  [Bdhm  et  al.,  1988;  Frenklach  et  al.,  1983]  and  care  must  be  exercised  in  generalizing 
results  obtained  over  a  specific  temperature  range  to  ail  hydrocarbon  fuels.  The  critical  point  for  the 
present  study  is  that  significant  temperature  sensitivity  has  been  observed  for  the  soot  volume  fraction 
for  temperature  changes  of  50  K  in  ethene  flames.  Thus,  the  observed  temperature  variations 
between  the  nitrogen  and  argon  diluted  flames  in  the  interior  regions  where  soot  is  formed  is 
potentially  important  in  terms  of  changes  in  the  observed  soot  volume  fraction. 

The  second  observation  involves  the  effect  of  diffusion  on  the  evolution  of  the  fuel 
concentration  field  when  diluent  is  added  to  the  fuel  flow.  The  previously  presented  results  show  that 
diffusion  processes  diminish  the  effects  of  diluting  the  fuel  at  locations  downstream  of  the  fuel  tube 
exit.  This  is  particularly  true  for  the  region  where  soot  particles  are  first  observed  to  form.  In  order 
to  further  investigate  the  generality  of  these  observations  a  numerical  modeling  study  was  conducted 
and  compared  to  the  mass  spectrometric  measurements  of  species  concentrations  in  the  methane/air 
diffusion  flame  described  earlier. 
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A  generalized  flame  sheet  model  for  laminar  diffusion  flames  developed  by  Mitchell  and 
coworkers  [Mitchell  et  al..  1980)  was  used  to  investigate  the  variation  of  fuel  concentration  as  a 
function  of  position  in  the  flame.  The  calculations  shown  in  Figure  8  were  carried  out  for  undiluted 
and  50%  nitrogen  diluted  ethene/air  flames  at  fuel  flow  rates  of  2.75  and  6.58  cm^/s.  The  lower  fuel 
flow  rate  conditions  correspond  to  flames  11  and  1.  This  figure  plots  the  calculated  fuel  mole 
fraction  as  a  function  of  height  along  the  flame  center  line.  In  each  case,  undiluted  or  diluted,  the 
fuel  mole  fractions  reach  similar  values  at  axial  locations  between  10  to  30  mm  above  the  burner  exit 
depending  on  the  fuel  flow  rate.  These  results  are  further  examined  in  Figure  9  where  radial  profiles 
are  presented  for  Flames  1 1  and  1  at  three  axial  locations  (0.3  mm,  2.7  mm  and  14.5  mm). 
Additionally,  Table  4  shows  results  tabulated  for  other  fuel  flow  rates  at  specific  axial  locations  for 
selected  radial  positions  (4.63  mm  and  the  center  line).  Comparisons  shown  in  Table  4  and  in  Figure 
9  indicate  that  in  the  region  where  soot  is  first  formed,  the  local  fuel  concentration  ratio  between  the 
diluted  and  undiluted  flame  is  typically  between  0.72  and  0.88.  Additionally,  the  effect  of  dilution  is 
slightly  more  pronounced  for  the  lower  volumetric  fuel  flow  rate  case.  For  example,  at  z  =  2.7  mm 
and  R  =  4.63  mm,  the  concentration  ratio  is  0.76  for  the  2.75  cmVs  case  as  compared  to  0.85  for  the 
6.58  cm^/s  flame.  Thus,  these  modeling  results  show  that  in  the  region  where  soot  is  first  formed,  the 
ratio  of  the  fuel  concentration  in  the  diluted  and  undiluted  flames  achieves  similar  values  to  those 
based  on  mass  spectrometric  measurements  in  the  methane  flames  previously  presented.  The 
experimental  data  shows  a  similar  diminishing  effect  of  the  initial  fuel  dilution  as  predicted  using  the 
flame  code.  This  effect  is  a  result  of  the  rapid  interdiffusion  of  nitrogen  from  the  air  with  the  fuel, 
mitigating  the  initial  dilution.  Consequently,  all  laminar  diffusion  flames  burning  in  air  undergo 
"dilution"  by  nitrogen,  rapidly  reducing  the  local  fuel  concentration  [see  also  Boedeker  and  Dobbs, 
1986a.  1986b].  From  these  results,  one  can  conclude  that  the  local  concentration  of  fuel  and 
combustion  products  will  quickly  achieve  similar  values  in  both  diluted  and  undiluted  flames. 
Furthermore,  concentration  effects,  based  on  initial  values  at  the  fuel  tube  exit,  will  not  reflect  the  true 
concentration  effects  in  the  soot  formation  region. 

An  additional  aspect  of  the  soot  formation  process  observed  to  be  affected  by  dilution  is  the 
change  in  time  available  for  soot  growth  or  the  soot  residence  time.  Figure  7  which  shows  the 
maximum  scattered  light  intensity  at  selected  axial  positions,  clearly  indicates  that  soot  inception  is 
delayed  with  increased  dilution  of  the  fuel  flow.  To  calculate  the  soot  residence  time  a  buoyancy 
driven  acceleration  of  3100  cm/s^  is  used  [Flower,  1989;  Santoro  et  al.,1987].  The  soot  residence 
time  is  then  calculated  assuming  a  flow  with  constant  acceleration  and  an  initial  velocity  determined 
by  the  flow  conditions  at  the  fuel  tube  exit.  For  the  diluted  ethene  flames  in  this  study,  a  decrease  in 
soot  residence  time  was  estimated  to  account  .or  between  2%  and  15%  of  the  reduction  in  Fy  from 
the  undiluted  flames,  assuming  the  soot  volume  fraction  increases  linearly  with  time,  with  the  largest 
reductions  corresponding  to  the  6.58  cm^/s  ethene  flames.  In  more  dilute  flames  such  as  the 
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nitrogen  diluted  ethene  flame  of  Santoro  and  Semerjian  [1984]  (see  Figure  7),  the  decrease  in 
residence  time  is  observed  to  be  more  significant,  accounting  for  27%  of  the  reduction  in  soot  again 
assuming  a  linear  relation  between  Fy  and  time. 

Although  the  reduction  in  soot  particle  residence  time  with  increased  dilu'ion  may  be  a 
temperature  or  a  concentration  effect,  two  observations  exist  which  suggest  that  the  decrease  in 
residence  time  is  due  to  a  reduction  in  temperature.  First,  the  height  at  which  soot  inception  occurs, 
as  defined  by  a  rapid  increase  in  scattered  light  intensity,  is  nearly  the  same  for  flames  with  equal 
calculated  adiabatic  flame  temperature  indicating  that  the  soot  particle  inception  location  is  more 
strongly  influenced  by  temperature  (see  Figure  7).  Additional  temperature  measurements  are  needed 
to  fully  confirm  this  effect.  Second,  soot  is  observed  to  form  at  slightly  larger  (Ai^O.25  mm)  radial 
coordinates  for  the  more  dilute  flames.  Larger  radial  coordinates  correspond  to  higher  temperatures 
and  lower  fuel  concentrations.  Thus,  an  increase  in  temperature  is  a  physically  more  satisfying 
explanation  for  the  radial  shift  in  the  soot  particle  inception  location. 

A  more  quantitative  assessment  of  the  effects  of  temperamre  and  fuel  concentration  can  be 
made  based  upon  the  dilution  experiments  previously  discussed.  To  separate  the  temperature  and 
concentration  effects  on  soot  formation,  a  functional  form  for  soot  growth  is  assumed.  The 
functional  dependence  on  temperature  and  concentration  is  taken  to  be  similar  to  the  form  used  by 
Axelbaum  et  al.  [1988] 

Eq-4. 

where  a  is  the  apparent  order  of  the  soot  growth  rate  with  respect  to  the  initial  fuel  mole  fraction  and 
Ea  is  the  apparent  activation  energy.  In  the  present  analysis,  the  dependence  on  initial  fuel 
concentration,  Xpo<  is  not  taken  to  be  first  order,  rather  a  is  estimated  from  a  fit  to  the  experimental 
data. 

Several  cautionary  remarks  need  to  be  made  about  the  use  of  such  an  expression.  Soot 
formation  in  diluted  diffusion  flames  is  affected  by  a  series  of  complex  interactions  due  to 
diffusional  mixing,  temperature  field  variations,  and  soot  formation  chemistry.  Furthermore,  the 
local  soot  volume  fraction  within  the  flames  is  a  result  of  soot  inception  and  subsequent  surface 
growth  which  are  described  by  two  different  temperature  sensitivities,  i.e.  apparent  activation  energies 
[Vandsburger  et  al..  1984]  In  addition,  temperature  measurements  in  this  study  have  shown  that 
even  under  equal  calculated  adiabatic  flame  temperature  conditions,  differences  in  temperature  in 
argon  and  nitrogen  diluted  flames  occur  in  the  region  where  soot  is  first  formed.  Qearly  such  a 
simple  expression  does  not  incorporate  enough  of  the  physical  details  of  these  processes  to  make 
general  use  of  this  relationship  to  describe  soot  growth  in  diffusion  flames  justified.  Instead,  the 
deduced  values  for  a  and  Ea  represent  purely  empirical  values  based  on  a  correlation  to  the  ethene 
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tlames  in  this  study.  These  results  are  then  used  to  assess  the  relative  importance  of  temperature  and 
concentration  effects  in  the  present  study,  which  can  be  extended  to  studies  under  similar  conditions. 

To  evaluate  equation  4  for  values  of  a  and  Eq  .  the  time  rate  of  change  of  the  total  soot 
volume  fraction  is  needed.  The  calculation  of  the  time  coordinate  is  carried  out  in  the  same  manner 
described  for  the  calculation  of  soot  particle  residence  time.  The  slope  of  a  linear  least  squares  fit  of 
the  Fv  vs.  time  data  in  the  soot  growth  region  was  taken  to  represent  dFyIdt.  Temperature  values  are 
the  calculated  adiabatic  flame  values  and  Xpo  is  the  fuel  mole  fraction  at  the  fuel  tube  exit.  Using  a 
nonlinear  least  squares  curve  fitting  procedure  applied  to  the  results  of  the  twelve  diluted  ethene 
flame  conditions  listed  in  Table  1,  values  of  a  and  Ea  were  determined  to  be  0.30±0.18  and 
94.5±25.5  kcal/mole,  respectively.  It  should  be  noted  that  the  use  of  calculated  adiabatic 
temperatures  results  in  apparently  large  activation  energies.  For  example,  if  an  apparent  activation 
energy  was  determined  from  the  premixed  flame  results  of  Bdhm  et  al.  [1988]  based  on  the  variation 
of/voo  over  the  measured  temperature  region  of  1550  K  and  1600  K,  one  calculates  a  value  of  68 
kcal/mole.  If  instead  the  calculated  adiabatic  temperature  of  2369  K  is  used  with  a  AT  of  50  K  as 
appropriate  to  an  undiluted  and  diluted  ethene/air  diffusion  flame,  a  value  of  150  kcal/mole  results. 
Similarly,  Vandsburger  et  al.  [1984]  reported  the  apparent  activation  energy  for  the  overall  soot 
growth  rate  for  an  ethene  counterflow  diffusion  flame  as  41  kcal/mole.  However,  because  actual 
measured  temperatures  in  the  range  1270-1480  K  were  used,  an  assessment  of  the  apparent  activation 
energy  using  the  higher  calculated  adiabatic  flame  temperatures  would  result  in  a  significantly  larger 
value  of  the  apparent  activation  energy.  Thus,  interpretation  of  the  physical  meaning  of  this 
activation  energy  needs  to  be  done  with  care.  Essentially  for  the  present  purposes  Eq  represents  a 
measure  of  the  temperature  sensitivity  of  the  soot  formation  process  as  observed  in  these  studies  and 
has  no  further  fundamental  chemical  kinetic  significance. 

The  relative  effects  of  concentration  and  temperature  can  be  deduced  from  this  expression  by 
making  comparisons  between  an  undiluted  and  diluted  flame.  For  example,  the  reduction  in  the  fuel 
concentration  term  comparing  an  ethene  flame  diluted  with  50%  nitrogen  to  its  undiluted  counterpart 
is  0.81  (=(0.5/ 1.0)®-^®)  while  the  corresponding  reduction  in  the  exponential  temperature  term  for 
temperatures  of  2310  K  and  2369  K  is  0.60.  The  ratio  of  the  reduction  in  soot  concentration  due  to 
temperature  effects  to  the  reduction  due  to  concentration  effects  (l-0.60)/(l-0.81),  indicates  a 
temperature  effect  which  is  2. 1  times  more  important.  It  should  be  noted  that  this  estimate  is  believed 
to  be  a  lower  limit  on  the  temperature  dependence  since  other  factors  which  would  likely  increase  the 
temperature  dependence  were  not  included  in  the  analysis.  Specifically,  thermocouple 
measurements  in  flames  which  have  equal  calculated  flame  temperature  show  that  in  the  region  where 
soot  is  being  formed,  the  argon  diluted  flames  actually  have  lower  temperatures.  Also,  neglecting 
changes  in  residence  time,  which  appears  to  be  more  strongly  influenced  by  temperature  than 
concentration,  would  likely  lead  to  an  underestimation  of  the  temperature  effect. 
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A  second  assessment  of  the  dependence  on  temperature  can  be  made  by  direct  comparison  of 
the  soot  formation  rate  in  pairs  of  flames  with  equal  concentration.  In  this  comparison,  the  parameter 
a  need  not  be  evaluated.  The  activation  energy,  Ea.  was  found  to  be  235,  105,  and  97  kcal/mole  for 
flame  pairs  1  and  la,  3  and  3a,  and  6  and  7  respectively.  The  larger  temperature  dependence 
displayed  for  the  2.75  cm^/s  flames  (1  and  la)  may  indicate  a  more  complex  relation  between  fuel 
flow  rate  and  the  dependence  of  soot  formation  rate  on  temperature  than  indicated  by  Equation  4. 
One  should  be  aware  that  since  these  values  are  based  on  comparison  between  only  two  flames  they 
are  subject  to  considerable  uncertainty.  However,  in  each  case  these  values  indicate  a  dominant 
temperature  effect.  These  results  can.  therefore,  be  viewed  as  indicating  an  upper  limit  on  the 
temperature  sensitivity  of  the  soot  formation  process  in  fuel  dilution  studies. 

The  relative  effects  of  concentration  and  temperature  on  the  soot  formation  rate  are  shown 
graphically  in  Figure  10.  In  this  representation,  the  measured  soot  formation  rate  for  diluted  flames 
normalized  by  th^ir  undiluted  counterparts  is  shown  with  the  corresponding  ratio  of  the  exponential 
temperature  term  in  Equation  4  for  the  diluted  flames  and  their  undiluted  counterparts.  The 
exponential  temperature  dependence  is  shown  for  apparent  activation  energies.  Ea,  of  40.  95.  and  150 
kcal/mole  for  the  calculated  adiabatic  flame  temperature  of  each  of  the  ethene  flames  in  this  study. 
The  dashed  line  represents  a  slope  of  unity,  or  the  total  reduction  in  the  soot  formation  rate  by  both 
temperature  and  concentration.  For  a  given  apparent  activation  energy,  the  reduction  in  soot 
formation  rate  due  to  temperature  reduction  is  the  difference  between  the  exponential  temperature 
dependence  (the  solid  line)  from  unity.  The  contribution  of  concentration  to  the  reduction  in  soot 
formation  rate  with  dilution  is  then  given  by  the  difference  between  the  solid  line  and  the  dashed  line. 
For  example,  for  an  activation  energy  of  40  kcal/mole,  concentration  is  nearly  twice  as  important  as 
temperature  in  reducing  soot  formation  in  diluted  flames.  Figure  10  also  shows  that  an  activation 
energy  of  150  kcal/mole  represents  an  upper  limit,  since  the  calculated  reduction  in  soot  formation 
rate  due  to  temperature  fully  accounts  for  the  measured  reduction  in  the  soot  formation  rate.  The 
present  results  argue  for  an  apparent  activation  energy  closer  to  95  kcal/mole  to  represent  the 
temperature  sensitivity  of  the  soot  formation  process  in  these  diluted  ethene  flames.  However,  this 
single  value  can  not  be  rigorously  justified  for  all  flames  conditions  as  the  above  reported  variation  in 
Ea  attests.  The  present  results,  however,  do  indicate  that  under  conditions  where  diluents  are  added  to 
the  fuel  flow,  the  temperature  effect  is  more  impx)rtant  than  the  concentration  effect. 

Any  attempt  to  further  quantitatively  narrow  the  estimate  of  the  temperature  sensitivity  of  the 
soot  formation  process  in  terms  of  an  activation  energy  is  complicated  by  the  nature  of  the  dilution 
experiments.  Dilution  of  the  fuel  results  in  only  small  changes  in  temperature  and,  as  shown  in  this 
work,  the  local  fuel  concentration.  Thus,  highly  accurate  measurements  of  the  changes  in  soot 
particle  concentration,  temperature,  and  gas  species  concentration  are  required.  This  observation 
highlights  the  difficulty  in  studying  concentration  and  temperature  effects  through  a  fuel  dilution 
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approach  where,  based  on  the  current  results,  coupled  interactions  involving  thermal  and  mass 
transport  are  seen  to  be  important.  The  complexity  of  this  flame  environment  may  contribute  to  the 
variety  of  interpretations  regarding  the  effects  of  inert  dilution  on  soot  reduction. 

Conclusions 

Measurements  of  the  decrease  in  soot  volume  fraction  at  various  dilutions  have  shown  that 
soot  is  reduced  considerably  as  both  the  adiabatic  flame  temperature  and  initial  fuel  concentration  are 
simultaneously  decreased.  Both  temperature  reduction  and  fuel  concentration  reduction  were 
considered  as  likely  explanations  for  the  reduction  in  the  soot  formation  rate.  An  examination  of  the 
decrease  in  the  soot  formation  rate  in  flames  with  equal  fuel  concentration,  but  different  degrees  of 
initial  dilution,  show  that  the  more  heavily  diluted  argon  flames  produce  less  soot  than  the  nitrogen 
diluted  flames.  Rapid  insertion  thermocouple  measurements  made  in  these  flames  near  the  region 
where  soot  first  forms  show  that  the  argon  diluted  flames  display  lower  temperatures  than  their 
nitrogen  diluted  counterparts.  This  point  attests  to  the  need  for  local  measurements  when  trying  to 
quantitatively  deduce  temperature  and  concentration  effects  in  diffusion  flames.  The  use  of  adiabatic 
temperatures  to  quantify  small  changes  in  temperature  due  to  dilution  appears  to  underestimate  the 
effect  of  temperature  when  comparing  soot  formation  rates  in  argon  and  nitrogen  diluted  flames  with 
the  same  calculated  adiabatic  flame  temperatures. 

The  results  from  this  study  suggest  that  the  role  of  temperature  on  the  soot  reduction  process 
exceeds  the  reduction  due  to  a  decrease  in  fuel  concentration.  These  results  are  consistent  with  the 
measured  and  calculated  effect  of  dilution  on  the  local  fuel  concentration;  that  is,  in  the  coannular 
geometry,  the  local  fuel  concentration  quickly  approaches  a  similar  value  regardless  of  the  initial 
dilution  conditions.  In  fact,  in  the  region  where  soot  is  first  formed,  diffusion  greatly  reduces  the  fuel 
concentration  ratio  when  comparing  diluted  flames  to  undiluted  flames  as  compared  to  the  initial  fuel 
tube  exit  conditions.  A  comparison  of  argon  and  nitrogen  diluted  ethene  flames  at  equal  dilution 
and  different  calculated  adiabatic  flame  temperatures  shows  that  temperature  significantly  alters  the 
soot  formation  rates.  One  difficulty  in  accessing  the  relative  contribution  of  temperature  is  that  very 
small  changes  in  temperature  need  to  be  measured  accurately.  In  addition,  dilution  experiments  offer 
only  a  narrow  range  over  which  temperature  can  be  modified  and  the  resulting  changes  in 
temperature  and  concentration  are  observed  to  be  small.  Although  the  present  analysis  does  not  lead 
to  a  singular  assessment  of  the  relative  effect  of  concentration  and  temperature  on  soot  formation 
rates  in  diluted  flames  (one  may  not  expect  this  as  these  effects  may  be  functions  of  fuel  flow  rate, 
fuel  etc.  for  these  diffusion  flames),  these  results  support  a  temperature  effect  which  is  more 
important  than  the  fuel  concentration  effect  on  the  soot  reduction  associated  with  inert  dilution  to  the 
fuel  flow  in  coannular  laminar  diffusion  flames. 
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Figure  Captions 


Figure.  1.  Uncorrected  thermocouple  measurements  of  the  temperature  profiles  as  a  function  of 
radial  position  for  several  axial  heights  in  an  undiluted  (flame  8)  and  diluted  (flame  12)  ethene 
diffusion  flame.  The  ethene  flow  rate  is  4.90  cm^/s  while  the  argon  diluent  flow  rate  is  4.50  cmVs 

Figure  2.  Comparison  of  uncorrected  thermocouple  measurements  of  the  temperature  profile  at 
an  axial  height  of  5  mm  as  a  function  of  radial  position.  Results  are  presented  for  an  undiluted 
ethene  flame  (flame  12)  and  two  pairs  of  diluted  flames  (flame  5  and  6  and  flames  7  and  8) 
which  have  identical  calculated  adiabatic  temperatures  (2346  and  2333  K  respectively). 

Figure  3.  The  total  integrated  soot  volume  fraction.  Fy,  as  a  function  of  height  in  the  flame  for  a) 
flames  12,  5.  and  6;  b)  flames  12,  7,  and  8;  c)  flames  12,  3.  and  4.  The  fuel  mole  fraction  is 
listed  in  parentheses  next  to  the  symbol  identification.  The  lines  represent  a  spline  fit  to  the  data. 

Figure  4.  The  total  integrated  soot  volume  fraction.  Fy  as  a  function  of  height  in  the  flames  for 
pairs  of  flames  (3  and  3a)  and  (6  and  7)  at  equal  dilutions  of  50%  and  64%  respectively.  The 
lines  represent  spline  fits  to  the  data. 

Figure  5.  Mass  Spectrometric  mole  fraction  measurements  of  N2,  Ot.  H2O,  and  CH4  as  a 
function  of  radial  position  in  an  undiluted  methane  diffusion  flame  burning  in  air.  The  axial 
measurement  location  is  12  mm  above  the  burner  exit  The  methane  fuel  flow  rate  was  5.7  cm^/s. 

Figure  6.  Mass  Spectrometric  mole  fraction  measurements  of  CO,  CO2,  H2,  and  C2H2  as  a 
function  of  radial  position  in  an  undiluted  methane  diffusion  flame  burning  in  air.  The  axial 
measurement  location  is  12  mm  above  the  burner  exit  The  methane  fuel  flow  rate  was  5.7  cm^/s. 
Note  negative  CO  mole  fractions  are  a  result  of  a  differencing  procedure  in  which  the  mole 
fraction  of  CO  and  N2  are  deduced  from  measurements  at  two  different  electron  voltages. 

Figure  7.  The  maximum  volumetric  scattering  cross  section,  Qw,  as  a  function  of  height  for 
diluted  and  undiluted  ethene  diffusion  flames.  Flow  rate  conditions  in  cm^/s  are  shown  in 
parentheses  next  to  the  symbol  identification. 

Figure  8.  Comparison  of  the  calculated  fuel  mole  fraction  as  a  function  of  height  along  the 
center  line  of  a  diffusion  flame  for  two  ethene  fuel  flow  rates  (2.75  cmVs  and  6.58  cmVs). 
Calculations  compare  an  undiluted  and  diluted  (50%  N^)  flame  case.  Flow  rates  in  cm^/s  for  C2H4 
and  N2  are  listed  in  parentheses  next  to  the  identifications  of  the  line  symbols. 

Figure  9.  Comparison  of  the  calculated  fuel  mole  fraction  as  a  function  of  radial  position  for  a 
diluted  (50%  Nj)  and  an  undiluted  ethene  diffusion  flame  at  heights  of  a)  0.3  mm,  b)  2.7  mm,  and 
c)  14.7  mm  above  the  fuel  tube.  Flow  rates  in  cm  Vs  for  C2H4  and  N2  are  listed  in  parentheses 
next  to  the  identifications  of  the  line  symbols. 

Figure  10.  Comparison  of  temperature  and  concentration  effects  on  soot  reduction  with  dilution 
for  the  ethene  flames.  Plotted  is  the  ratio  of  the  exponential  temperature  dependence  of  the 
diluted  and  undiluted  flames  versus  the  ratio  of  the  soot  formation  rate  for  the  corresponding 
diluted  and  undiluted  conditions.  The  dashed  line  represents  the  total  reduction  in  soot  formation 
rate.  The  contributions  of  temperature  and  fuel  concentration  are  shown  for  Eg  of  40,  95,  and 
150  kcal/mole.  The  contributions  of  the  temperature  and  concentration  dependence  for  Ea=40 
kcal/mole  are  indicated  by  the  arrows.  The  present  results  support  a  value  of  Ea  close  to  95 
kcal/mole. 
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Table  1.  Flow  conditions  for  the  diluted  and  undiluted  laminar  flames  specifying  fuel, 
fuel  flow  rate  (Q),  diluent,  initial  fuel  mole  fraction  (Xp)  and  calculated  adiabatic  flame 
temperature  (Tad)-  The  total  integrated  soot  volume  fraction,  F„  and  the  ratio  expressed 
in  percentage  of  the  measured  F,  to  the  F,  for  a  pure  undiluted  fuel  case  is  also  given. 


Exp# 

Fuel  (Q) 
(cm^/s) 

Diluent 

Xf 

Tad 

(K) 

Fv 

(cm“) 

%  of  Pure 

1 

C2H4(2.75) 

N2 

0.50 

2310 

1.66E-06 

46.7 

la 

C'>H4(2.75) 

Ar 

2330 

2.29E-06 

64.7 

2 

C2H4(2.75) 

At 

0.37 

1.05E-06 

29.6 

3 

C2H4(4.90) 

N2 

0.50 

2.91E-06 

51.8 

3a 

C2H4(4.90) 

At 

mSM 

3.70E-06 

65.8 

4 

C2H4(4.90) 

Ar 

0.37 

2310 

2.12E-06 

37.7 

5 

C2H4(4.90) 

N2 

0.74 

2346 

4.90E-06 

87.2 

6 

C2H4(4.90) 

At 

0.64 

2346 

456E-06 

81.2 

7 

C2H4(4.90) 

N2 

0.64 

4.08E-06 

72.7 

8 

C2H4(4.90) 

At 

0.52 

3.77E-06 

67.2 

9 

C2H4(6.58) 

N2 

0.64 

2333 

5.30E-06 

74.5 

10 

C2H4(6.58) 

At 

0.52 

4.48E-06 

63.0 

11 

C2H4(175) 

— 

2369 

334E-06 

100.0 

12 

C2H4{4.90) 

- 

2369 

5.62E-06 

13 

C2H4(6.58) 

•• 

1.00 

2369 

7.10E-06 

100.0 
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Table  2.  Comparison  of  thermocouple  temperature  measurements  at  an  axial  position  of 
5  mm  above  the  fuel  tube  for  ethene  flames  diluted  with  argon  and  nitrogen  at  the  same 
calculated  adiabatic  flame  temperature.  Maximum  temperatures  (Tmax)  as  well  as 
temperature  differences  (ATr.5„„)  between  the  nitrogen  and  argon  diluted  flames  (Tn2* 
Ta,)  are  shown. 


Exp. 

# 

Fuel  (Q) 

(cm^/s) 

Diluent 

Xf 

Tad 

(K) 

Tmax 

(K) 

^TR=5niin 

(K) 

1 

C2H4(2.75) 

N2 

0.50 

2310 

1934 

56 

2 

C2H4(2.75) 

At 

0.37 

2310 

1924 

5 

C2H:4(4.90) 

N2 

0.74 

2346 

1994 

33 

6 

C2H4(4.90) 

At 

0.64 

2346 

7 

C2H4(4.90) 

N2 

0.64 

2333 

1988 

56 

8 

C2H4(4.90) 

At 

0A2 

2333 

1988 
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Table  3.  Comparison  of  fuel  concentration  effects  for  flames  at  similar  calculated 
adiabatic  flame  conditions. 


Exp. 

# 

Fuel  (Q) 
(cm^/s) 

Tad 

(K) 

AXp* 

% 

AC-  ** 

AFv 

% 

.  *** 

D 

1  &2 

C2H4(2.75) 

2310 

26.0 

36 

lAO 

3  &4 

C2H4(4.90) 

2310 

26.0 

27 

1.06 

5&6 

C2H4(4.90) 

2346 

13.5 

6.9 

0.36 

7&8 

C2H4(4.90) 

2333 

18.8 

7.6 

0.38 

9&10 

C2H4(6.58) 

2333 

18.8 

15 

0.81 

14&15 

C3H8(2i6) 

2240 

18.0 

7.7 

0.40 

*  AXp  is  the  difference  in  fuel  mole  fraction  and  is  based  on  a  comparison  of  Ar  and  Nj 
diluted  flam  at  the  same  calculated  temperature. 

**  AFv  is  the  difference  in  total  integrated  soot  volume  fraction  and  is  based  on  a 
comparison  of  Ar  and  Nj  diluted  flames  at  the  same  calculated  temperature. 

***  b  is  based  on  assuming  AF,,  is  proportional  to  AXp. 
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Table  4.  Mitchell  flame  model  comparisons  for  fuel  mole  fraction  (Xp)  in  50%  diluted 
and  undiluted  ethene  diffusion  flames  at  selected  radial  (R)  and  axial  (Z)  coordinates. 


Fuel  Flow 
Rate 
(cm^/s) 

R 

(mm) 

Z 

(mml 

XplPure) 

Xf(50%) 

Xp(50%) 

Xp(Pure) 

6.58 

4.63 

0.3 

0.482 

mm 

o.m 

6.58 

4.63 

2.7 

0.256 

0.849 

6.58 

0.00 

14.5 

0.504 

0.401 

0.195 

4.90 

4.63 

0.3 

0.410 

0.304 

0.740 

4.90 

4.63 

2.7 

0.207 

0.170 

0.823 

4.90 

0.00 

14.5 

0.383 

0.330 

0.860 

2.75 

4.63 

0.3 

0.298 

0.220 

0.738 

2.75 

4.63 

2.7 

0.132 

0.101 

0.764 

2.75 

0.00 

14.5 

0.215 

0.189 

0.879 
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Fuel  Mole  Fraction 
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OPTICAL  CROSS  SECTIONS  FOR  AGGREGATES 

H.  A  DOBBINS 
Division  of  Eufiincerina 
Brown  I'niversitij 
Providencr.  Bl  (>2f)I2 


H  J.  SANTOBO 

Department  of  Mechanical  P.niiincentiu 
The  Pennsylvania  State  I'niversitii 
Vniversity  Park.  P.\  I(i.S()2 

\\l) 

H.  C.  SEMERJIAN 
Center  for  Chemical  Technology 
Sational  Institute  of  Stanilanis  ami  Tcchnoloiiy 
Caithcrslnins.  MD  2t)SW 


Soot  tbrmfcl  in  tluint-s  ii.sually  consists  ol  aujjrciiatcs  oliislcrs  or  auiilomcratcs)  ol  a  \arial)lc 
niiinl)cr  ot  nearly  spherical,  inonoilisperse  primary  particles  unonoiners  or  sphenilesi.  In  this 
work,  the  optical  |)roperties  of  polydis|wrs«*  anKretjatcs  are  used  to  analyze  liiiht  scatterini: 
data  troTii  a  coannular  ethene  diffusion  flame.  In  previonsU  reported  studies,  tlata  have  been 
obtained  on  the  local  extinction  and  volumetric  scattcrinii  cross  sections  from  laser  scatteriim 
experiments,  on  the  flame  velocity  field  from  la.ser  velocimeirs ,  and  on  the  primary  particle 
sizes  determined  by  electron  microscopy.  The  present  anaKsis  \ields  the  aseraue  number  of 
primary  particles  per  aggregate,  the  mean-square  radius  of  gyration,  the  soot  volume  fraction 
and  the  aggregation  rate,  ft  is  found  that  sustaiiu-d  collisional  growth  of  the  aggregates  occurs 
while  their  primary  particles  grow  through  helerogeneons  reactions  low  in  the  llame.  and 
contract  through  surface  oxidation  in  the  upper  half  of  the  llame.  .\  recent  \alne  of  the 
refractive  index  gives  internally  consistent  iioment  ratios  of  the  aggregate  size  distribution 
lunction.  This  method  ot  analysis  pro'  ides  .'  more  detailed  .md  lomplete  description  ot  the 
tirrmation.  growth  and  oxidation  of  soot  aggreuates  in  a  diffusion  llame. 


Introduction 

Tlic  use  ol  laser  liglit  scattering  experiments  to 
probe  the  development  of  the  soot  field  in  the 
coannular  laminar  ethene  difiusion  flame  has  been 
previously  reported. '  This  work  was  based  on  the 
progress  made  In  various  investigators  to  measure 
the  optical  properties  of  the  soot  particles  field  with 
high  spatial  resolutioir  '  and  the  application  of  a 
tomographic  Inversion  fechnitjue’  to  obtain  local 
values  ot  the  extinction  coefficient  for  axisymetric 
diffusion  flames.  The  measured  extinction  and  vol- 
uinetric  scattering  cross  sections  were  introduced 
into  expressions  derived  from  the  Rayleigh  descrip¬ 
tion  for  the  scattering  of  light  by  small  spheres. 
Polydispersifv  was  included  by  assuming  the  size 
distribution  to  be  self  preserving.  However,  Rav- 


leigh  theorv  im|)hes  tiiat  the  vcrticallv  polarized 
comiXHient  of  scattering  is  independent  of  the  scat¬ 
tering  angle  ()  measured  from  the  forward  direc¬ 
tion.  This  condition  is  not  fulfilled  in  the  annular 
region  ol  the  llame  where  high  soot  concentrations 
were  observed  and  the  scattering  in  the  forw.ird  di¬ 
rection  was  more  pronounced.'’ 

.Subsequentiv  .i  detailed  mapping  of  the  ll.iine 
velocity  and  temperature  fields'  provided  iulorm.i- 
tion  on  the  soot  transport  and  growth  |trocesses  In¬ 
formation  on  particle  morphologv  was  ,ilso 
developeiP  In  means  of  a  theruiophorctic  s.un|ilmg 
teehni()ue.  Tliese  ri'sults  showed  the  soot  partules 
in  the  ethene  diffusion  llame  to  consist  ot  aggre¬ 
gates  of  primarv  particles  whose  shape  vvas  .ipitrov- 
imatelv  spherical  and  whose  diameters  r/,,  in  .mv 
one  flame  location  was  quite  uniform.  The  optie.il 
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properties  of  anjiregates  was  further  studied  throiitdi 
eomputational  simulation'  of  both  the  eluster-elus- 
ter  growth  and  the  light  scattering  processes.  Tliesr- 
results  together  with  the  observations  of  the  soot 
field  by  electron  microscopy  (TEM)"’  '"  were  used 
to  formulate*'^  the  optical  cross  sections  of  i-wjlytlis- 
perse  aggregates.  Polydispersits'  results  when  the 
number  n  of  primary  particles  per  aggregate  is  vari¬ 
able  within  a  given  aggregate  population,  and  is  the 
inevitable  result  of  the  cluster-cluster  aggregation 
process. 

These  developments  now  make  possible  the  anal- 
\'sis  of  the  measured  optical  properties  of  soot  using 
the  cross  sections  for  randomly-oriented  aggregates. 
The  presence  of  angular  dissv  mmctry  of  the  scat¬ 
tered  light  provides  the  opportunity  to  perform  a 
more  detailed  study  than  would  otherwise  be  pos¬ 
sible.  .\dditionally,  the  accumulated  data  on  this 
flame — scattering  and  extinction  mexsurements.  ve¬ 
locity  data  from  laser  Doppler  anemometry.  and 
primary'  particle  sizes  as  determined  by  TEN! — in 
the  region  of  maximum  soot  concentration  allow  for 
a  detailed  examination  of  the  dynamics  of  the  pro¬ 
cesses  which  influence  the  soot  field.  In  particidar, 
with  the  development  of  an  aggregate  mcKlel  for  the 
optical  cross  sections,  the  eflects  of  the  aggregation 
process  on  particle  morphology,  surface  area  and 
number  concentration  can  be  investigated.  The 
combination  of  detailed  optical  measurements  along 
w  ith  an  aggregate-based  analysis  provides  an  im¬ 
proved  understanding  into  the  evolution  of  s(M)t 
particle  fields  in  flames.  Given  the  recent  interest 
in  the  modeling  of  the  formation  and  growth  of  soot 
particles  and  their  relationship  to  flame  radiation, 
improved  approaches  for  the  interpretation  of  tip- 
tical  measurements  of  soot  aggregates  is  a  topic  of 
significant  importance. 

Theoretical  Background 

.Aggregate  Statistics: 

•A  single  aggregate  is  characterized  by  the  pri- 
marv'  particle  diameter  r/,„  the  number  of  (irimary 
particles  within  the  aggregate  n,  and  its  radius  of 
gyration  that  is  given  by 


where  r,  is  the  distance  of  the  ith  primary  particle 
to  the  center  of  mass  of  the  aggregate,  fjrge  ag¬ 
gregates  are  mass  fractals'^"’  which  display  scale 
similarity  and  obey  the  relationship 


where  kf  is  the  prefactor  anil  1)(  is  the  fractal  or 
Ilausdorf  dimension.  When  aggregates  grow  In 
cluster-cluster  collision  processes,  the  fractal  di¬ 
mension  is  in  the  range  of  1.7  to  1.9.  .Although  Kip 
i2)  applies  to  a  ixipulation  of  particles  on  a  st.itis- 
tical  basis,  even  relatively  small  aggregates  in  =  5' 
that  do  not  display  scale  similaritv  have  been  fouiui 
to  follow  this  relationship.'" 

In  order  to  describe  a  population  of  aggregates, 
it  is  necessarv  to  introduce  the  probability  distri¬ 
bution  function  iPDE'  pin)  that  represents  the  rel¬ 
ative  lret|uencv  of  occurrence  of  aggregates  lon- 
taining  n  primarv  particles,  fhe  i/' '  moment  of  the 
PDF  is  given  by 

ii'^  =  V  n'  ;)(n)  lA' 

Furthermore,  lor  each  value  of  n  there  is  an  as¬ 
sociated  volume-e(|uivalent  diameter  D  and.  lor  the 
total  aggregate  population,  the  corresponding  vol¬ 
ume-mean  diameter  Di„  is  given'"  by 

1)!,,  -  i/;v 

where  is  the  average  number  of  primary'  parti¬ 
cles  per  aggregate.  Since  the  diameters  of  the  pri¬ 
mary  particles  in  a  given  region  of  a  laminar  flame 
are  narrowly  distributed  (see  Ref.  12  and  the  ref¬ 
erences  cited  therein),  in  this  work  we  treat  the 
primary  particles  as  monodisperse  at  each  flame  lo¬ 
cation.  With  this  definition,  the  volume  fraction  /\ 
of  an  aerosol  of  N„  ansiraiatrs  per  unit  volume  can 
be  expresseil  as 

A  =  -  ''A' 

() 

The  volume  fraction  is  a  ijuantitv  tor  wliich  optical 
techniques  provide  direct  information  and  is  a  pri¬ 
mary  quantitv  of  interest  in  the  present  analvsis. 

Finally  a  series  of  relationships  exist  relating  pri¬ 
mary  and  aggregate  particle  properties  to  e.ich  other 
and  to  important  observables.  .A  moment  ratio  that 
is  important  in  the  reduction  of  light  scattering  data 
is which  is  a  measure  of  the  width  of  the  PDF. 
and  is  defined  by 


F’or  monodisperse  aggregates  /„  is  nnitv  and  it  at¬ 
tains  a  value  near  2.0  lor  the  contmuons  self  pre¬ 
serving  distribution  function,  rhis  (|uantitv  explic¬ 
itly  enters  into  the  analysis  of  light  scattering 
measurements  for  determining  Dy,  which  is  de¬ 
scribed  in  the  next  section. 
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The  number  ot  priiuar\  particles  per  unit  volume 
can  be  related  to  the  number  ot  auareijates  bv 

N,,  =  u'  N„  .7) 

and  the  suriace  area  per  unit  volume  S,  can  be  cal¬ 
culated  it  it  is  assumed  that  the  primary  particles 
are  in  point  contact  with  one  another,  \  ieldinu 

S,  =  ■ndj.n'  N„.  iS) 


tion  ot  polydisperse  auiireuates.  ‘  ’  The  determina¬ 
tion  of  p^.,,  which  is  otten  small  but  is  significant  at 
inidflame  heights  in  the  ethene  dittusion  tlame,  is 
given  in  .\ppendix  11. 

Because  the  pnxfuct  ot  i?.  and  N„  appears 
in  E(|.  (10).  the  particle  soluine  fraction  is  directly 
obtainable  from  a  measurement  ot  extinction  with¬ 
out  any  requirement  of  knowledge  of  the  size  dis¬ 
tribution  tunction.  By  combining  Kqs.  i.5)  and  (lOK 
we  find  the  particle  xolume  fraction. 


.Since  surface  growth  leads  to  the  tormation  ot  men¬ 
isci  that  reduce  the  surface  area,  the  \alue  of  S, 
vielded  by  Eq.  (8)  is  an  upper  bound  estimate. 
IMipulation  of  polydisperse  aggregates  has  a  mean- 
sipiare  radius  of  gyration  Rj  that  is  defined  by 
weighting  R"(n)  by  n~.  see  .Appendix  1.  The  final 
important  relationship  relates  the  mean-square  ra¬ 
dius  of  gy  ration  R;  with  d,,  and  D  ),,, 

.  f^w 
K  R: 

wlu're  H  is  a  high  luoment  ratio  of  tlu-  I’DE  (!<•- 
fined  in  .Appendix  I. 

The  above  ecjuations  contain  eight  properties  rel¬ 
evant  to  the  description  ot  a  suspension  of  aggre¬ 
gates — </,„  Djo,  f\  .  in.  Rj.  and  S,.  .As  five 

algebraic  relationsh^s  exist,  (Eqs.  (4),  i.5i,  (7).  (8' 
and  (9)]  and,  it  /„,  H.  and  Dy  are  known,  there  are 
three  degrees  of  freedom.  Measurements  of  three 
optical  properties  which  are  independent  functions 
ot  three  or  more  ot  the  above  eight  properties  are 
required  to  completely  specify  the  aggregate  jMip- 
ulation  when  the  refractive  index  of  the  aggregate 
material  is  known.  In  the  procedure  employed  be¬ 
low.  we  will  consider/,,  and  K.  which  are  not  iii- 
ilependent.  as  unknown  and  will  use.  as  substitute 
information,  a  determination  of  (/,,  bv  TEM. 

Optual  Propertivs  of  Aiiiirruatcs: 

The  extinction  cix-tficiiMit  of  aggregates  which  Ixith 
absorb  and  scatter  light'*  is  given  bv 


p,,i 


lirN,,!!'  .r,',Eimi  _ 

- - ,  1  *  ,  10) 


Eim,  =  -Ii 


with  .tp  —  TTf/p/X.  k  -  2  it/\.  Ill  —  iiir  —  iirii  is 
the  refractive  index,  and  is  the  ratio  of  the  scat¬ 
tering  to  absorption  cross  sections  for  the  popiila- 


■JkEimid  +  p^j,  5  k  Elm) 

The  volumetric  scattering  cross  section  Q\x  (O)  for 
vertically  polarized  scattered  light,  when  the  inci¬ 
dent  beam  is  also  vertically  polarized  with  respect 
to  the  scattering  plane,  is  the  product  of  N'a  and 
the  mean  cross  section  of  the  aggregates.  It  is  given 
bv 


y .V  ,Ol  - 


hi:  R;) 


Fim)  = 


nr  -  1  i 


and  q,  is  the  modulus  oi  the  scattering  vector  de¬ 
fined  by 

llT  (f, 

(/,  =  —  sin  —  1 15) 

K  2 

The  tunction  /u/r  R|;)  results  Ironi  the  e.xamination 
of  the  angular  scattering  ot  light  by  aggregates"*  and 
its  lunctional  form  depends  on  q:  R-  and  f),  as 
summanzcHl  below  For  small  and  intermediate  sizes. 
<i:  R;;  <  I  ..5  1),.  the  function  fu/r  R''  is  given  by 


f\q:  R;:‘  =  -■\p| 


while  tor  larger  sizes,  ip  R;: 
expression  applies. 


I  5  I),.  the  following 


q,  R. 


where  (),,  =  ,1.5  D,/ei"' '  =  1)  99.19  for  1)^  =  l.SO 
as  is  representative  tor  cluster-iiuster  aggregation 
Eipiations  ,  16)  and  <  17)  result  Irom  the  detailed 
investigation'*  of  the  iiitlucnce  of  |x)lydispersity  on 
the  optical  properties  of  randomly  oriented  aggre¬ 
gates.  and  thev  are  related  to  expressions  which  ap- 


0  99.W  for  I),  =  l  .SO 
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|K'ar  ill  the  literature  ot  fractals.**"’  The  inclusion 
ol  the  ettects  ot  polyclispersity  of  the  a(»greijates  is 
an  imperative  consideration  in  iinderstandint;  the 
soot  formation  processes  in  tlames. 

The  mean-square  radius  of  gyration,  R^,  is  di¬ 
rectly  yielded  by  solving  Eqs.  (13),  (16)  and  (17)  as 
appropriate  for  the  ratio  (0.)/Qvv  iB,'. 

which  depends  onl\  on  k~  Rn  for  specified  values 
ol  ^  and  0,  since  all  other  quantities  cancel.  When 

‘l7  K;  <  1.5  iRj,  <  hii),  Eq.  (16)  is  applicable 
.ind  we  find 

R:  =  ptnR„.  ,181 

Somewhat  larger  aggregates  will  result  in  qf  rJ,  > 

1.5  D,(R„  >  ('ll)  and  Rn  is  then  found  by  using  both 
Eqs.  1 16)  and  (17)  to  formulate  Ry.  In  this  event, 
it  becomes  necessars  to  solve  the  transcendental 
('(piation 


'()'  R;))  e\|) 


i 

fMJ 

1 

\  1.5  IV  ' 

\  dj 

2/D, 


,lh. 


The  \alues  of  the  quantities  a,,,  by,  Cy,  and  d,i  are 
defined  in  general  terms  in  Table  1  where  they  are 
also  tabulated  fior  0,  =  4.5°  and  0y  equal  to  either 
)K)°  or  1.35°  tor  a  fractal  dimension  Df  =  l.H.  Two 
determinations  of  R;;  \ielded  by  the  use  of  rivo  sets 
ot  diss\mniefr\'  ratios  provide  replication  of  this  da¬ 
tum.  but  they  do  not  afford  independent  informa¬ 
tion  about  the  aggregate  population.  In  the  data 
anab'sis  givenjielow .  we  report  the  average  of  the 
two  values  ot  R;:  found  from  the  dissymmetrv  ratios 
Ri2  and  R|v 

We  find  the  solume  mean  diameter  by  dividing 
E(j.  (10)  by  (1.3)  aiid  using  Eq.  (6)  to  eliminate  ir 
in  favor  of/„  and  n*.  The  result  is 


\  j  fTrEdii)  (,)v\  10,1  ' 

TT  I  f’,  I'i/ii’/iq,"  R:,)lv,'„,' 


I'rom  K(|s.  ,5.  anil  20)  tor  f\  .mil  Dn)  we  liiid  \ 
w  hich  can  be  e.spressetl  as 


f-n-.r;,. 


where  .Vi,, 

W  hen  i/rRr,  1.  Eq.  il3'  reduces  to  the  i  siires- 
sion  applicable  to  Rayleigh-si/.e  spheres  even  though 
the  particles  are  actually  aggregates  of  small  spheres. 
In  this  particular  case.  Rj  and  i/,,  do  not  appear  in 
the  equations,  and  their  values  cannot  be  deter¬ 
mined  from  the  above  measurements.  .Although  the 
values  of  f\.  Dyi,  and  ire  yielded  by  the  e\- 
tinction-seattering  measurements,  no  direct  iiifor- 
niation  on  whether  the  observed  particles  are  spheres 
or  aggregates  is  provided.  Independent  ohserv.i- 
tions  .ire  then  required  in  order  to  determine  the 
.letiial  particle  morphology.  Without  this  .idded  m- 
lorniation  the  determination  of  .S,  is  impeded. 


Data  .Analysis  and  Discussion 

'Vbe  rtame  Imrner  which  is  described  elsewhere 
consisted  of  two  coannular  brass  tubes  ot  1 1 . 1  .mil 
101.6  mm  i.d.;  ethene  flows  throiigli  the  central  tube 
■It  the  rate  of  .3.85  cm ’/s  while  the  air  How  is  713 
em’/s  corresponding  to  the  nonsooting  (NS)  flame 
previously  described.’  In  this  e.xperiment  a  laser  light 
scattering  apparatus  is  ii.sed  to  determine  the  local 
extinction  K,.,,  anil  the  vohunetric  scattering  cross 
sections  v  (0|)  at  0,  =  45°.  ))()°.  and  13.5°.  These 
quantities  were  nieasured  .is  a  funetion  ot  height 
.iloiig  a  soot  p.irtii'le  path  in  the  .imiiilar  region  ot 


T,il,le  I 

Value  ol  C. (instants  lor  I  se  in  the  Determination  ol  tlie  Me. in  .Square  Iv.ulins  ol  ttv  ration  Ironi  l.ip 
IS)  .uid  (19)  for  O,  =  4.5°.  anil  D  ‘Mr  or  1  >.5''  when  D,  =  I  SI) 


r„ 

/>,. 

1 

(/., 

(j  =  12 

1)  5412 

2.121 

1  S)M) 

)  020 

i  039 

1/  =  13 

0.4142 

1.061 

2.  lOS 

t  SS(i 

»  910 

rhe  ahovt*  (juantitu’s  art*  (iffined  bv 


=  sin  (0,/2)/sin  (B,/2)  a.  ^  O  ,  |  mu 

h.,  =  exp[().5D,(l  -  r^l]  r..  ^  r,.'  " 

<L  =  rJC„  C,.  =  (15  I),/.')'”- 
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limh  MM)t  coiiccntratioii  in  tlie  toanmilar  latninar 
tlamf.  In  the  lollnwinu  analysis  we  take  advantaue 
nl  tliree  htnlies  nt  experimental  data — the  seattcr- 
ine-extinction  measurements,'  the  laser  veloeity 
data, '  and  obserx  ations  ot  primary  particle  size  ob¬ 
tained  by  thermophoretic  sampling  and  TEM  ex¬ 
amination.''  "  The  availability  ol  this  information  for 
the  ethene  diffusion  flame  allows  the  soot  field  to 
lie  examined  alontj  a  soot  particle  path  which,  ow- 
ine  to  the  action  of  thermophoresis,  does  not  co¬ 
incide  with  a  Has  streamline  in  the  lower  (lortion  of 
the  flame  where  velocities  are  lower  and  thermal 
Gradients  are  larne.  In  addition,  the  velocity  data 
provides  a  conversion  of  heiHht  above  the  burner, 
;,  to  time  and  allows  differentiation  of  certain  cpian- 
tities  to  obtain  rate  data. 

The  ({uantity  ft  contains  a  ratio  of  hiHher  mo¬ 
ments  of  the  angreuate  PDF  about  which  only  lim¬ 
ited  information  is  available.  The  values  of/„  and 
ft  are  not  independent  as  both  are  influenced  by 
the  width  ot  the  angreitate  PDF.  We  have  there¬ 
fore  analyzed  the  <lata  In  treatine  ^  and  f„  as 
i|uantities  to  be  adjusted  to  nive  the  best  auree- 
ment  between  </,,  based  on  optical  measurements 
nsiiiH  E<|.  lb'  and  from  TEM  oliservatioiis.  The  val¬ 
ues  of  Dx)  and  \„  are  independent  of  'M  and  are 
weakly  influenced  In  /,,,  while  d,,  is  sensitive  to  /„ 
and  to  ft  that  has  a  broader  ranije  of  variation.  On 
the  other  hand,  the  (juantities  fy  and  Rj  are  cal¬ 
culated  without  reference  to  k. 

The  values  of  d,,  obtained  from  Eq.  (9)  are  also 
influenced  bv  the  refractive  index  m  that  is  used  in 
the  analysis  of  this  data.  For  a  Hiven  value  of  m, 
we  find  a  unicpie  value  of  the  product  ( f„k)  that 
vields  the  minimiim  rms  of  (d,,  iTE.M)  -  d,,  (opt.'] 
=  Cr.  In  Table  11.  we  list  the  optimum  values  of 
f„  and  M  that  are  mutuallv  consistent  with  a  simjle 


value  ot  the  Heometrie  mean  vt.uulard  deviation  <r . 
lor  a  discrete  Ion  normal  distribution  fliree  values 
of  refractive-  index*'  .ire  used  in  this  anaivsis.  In 
Table  II  we  tind  ineonipatible  \, lines  of  f„  and  H 
result  vvlu-ii  III  -  I..j7  -  i0..5(i.  On  the-  other  hand, 
when  the  leeeiit  value'’  ot  in  =  2.10  —  /0..t5  is 
used.  /„  —  1..J  and  ft  =  .5.1  are  loiiiid  to  be  mii- 
tiially  eonsisti-nt  with  a  narrow  distribution,  ir .  - 
0.1.5.  This  combination  of  parameters  was  sc-leeted 
to  obtain  the  data  used  in  Fins.  1  to  4  for  wliieli 
C.  =  IS-T. 

Recentiv  it  has  been  reported'""'  that  the  v.ihie 
of  m  depends  on  the  hvclroHen  content  ol  the  soot. 
We  have  explored  this  possibility  bv  treatiiin  ot  the 
refractive  index  to  be  heiiiht  dependent  to  simulate 
the  annealinii  action  of  the  hot  flame  uases.  The 
retractive  index  in  is  allowed  to  varv  wi-akiv  vvith 
heiuht  r  and  achieves  the  value  ol  the  in  =  2. 10 
-  iff. 55  near  the  tip  of  the  flame  as  described  in 
the  IcKitnote  in  Table  11.  In  addition.  H  is  varii^ 
to  be  eonsistent  with  the  vertical  variation  of  ii' 
shown  in  Fiu.  5  below  tiir  a  Ion  normal  distribution 
with  (T^  -  0  :50.  rhese  interestiiiH  results  .ire  shown 
in  the  last  line  ol  Table  II  and.  .ilthoiiuh  tentative, 
they  siiHHest  .l  liitiire  direction  tor  research. 

Fitiure  I  shows  the  values  of  N„  and  D^j  ob¬ 
tained  from  Havieiuh  theory  ,ind  from  aniirenate 
theory  plotted  anainst  heiiiht.  .Vt  fi  =  2.5  ms  the 
()iiaiititv  /:(/■  R:'  is  equal  to  iinitv  .md  the  two 
methods  uive  the  same  values  tor  these  two  prop¬ 
erties.  The  diHerences  between  the  values  ot  Dm 
and  \„  that  are  predicted  bv  the  two  methods  is 
significant  at  all  other  times. 

The  quantities  /V.  Dxi.  and  N,,  are  displayed  in 
Fjn.  2,  and  the  (jiiantities  R^..  D«),  (/,,iTEM'.  and 
u'  are  niven  m  Fin.  5.  These  quantities  follow  .i 
consistent  p.ittern.  Hie  soot  voliniu-  traction  f\  in- 


Table  II 

\  allies  ol  H  and  K.,  recovered  from  minimization  ol  lrac  fioii.il  H\I,S  nl  li/,.  I'EM  il  ,  opt  |  lor 
selected  values  ol  refractive  index,  in  .il  \  =  "514  •')  mil. 


tn 

Ret 

J: 

it 

iT 

1  ' 

K 

2' 

1  57  -  lO  .50 

17 

I  t) 

10  1 

4'JO  ■  1  "'.  III  /s 

1  'JO  -  i0.,>5 

IS 

I  t 

50 

005 

ISO 

2  10  -  (0  .55 

I'f 

1  .5 

5  1 

O  1.5 

fMI  CJ 

— 

1  44 

2.‘J7  to  1  05 

0  50 

T'J 

note  5 

mott'  4' 

1  \  allies  ol  ueonietric  mean  standard  dev  iation  <r  .ire  estiin.ileci  Iroin  disc  rete  loi;  iioniial  clistribiitioiis 
which  vielil  the-  iiiclic .itecl  values  ot  f„  and  ft 

2  K,  is  the-  .miireiiation  r.ite  as  defined  bv  Kc).  li'}). 

5  Refractive  index  varies  accordine  to  in<z>  =  in.  I  -  i  nr  .v  z/li'""'  ■  I  with  m.  2  10.  m 

0  .5.5.  and  li  =  70  mm  corresponding  to  the  highest  obsc-rvation.  _ 

4.  The  values  ot  and  ft  are  cleterinined  to  be  consistent  with  n  i~i  shown  in  Fiu  5  bv  .issminnii  .i 
discrete  lou  normal  distribution  with  rr.  =  0..'40  and  K,  =  5  S 
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Fir.,  1.  C.'omparison  of  Data  IttHliutioii  nl 
l«Mi;h  Theory  for  Spheres  and  .\i'Kreuat<-  I'lieors  lor 
the  (,'oanniilar  Laminar  Diihision  Flaim*  (NS  ll.ime, 
(^)  =  .1.8.5  cin’/s).  (N,  and  D-^,  are  plotted  sersiis 
residence  time  t.  Time  cipials  zero  .it  the  Inniiei- 
inonth.  I 


•  *  •  :  0  -  . . 

20  *c  4. 

r:  -^s 

Fir.  2.  The  (.Quantities  D„,.  /, .  .ind  N  \s  t  .is 
Di'ternnneii  I  sinii  ( iross  .Sections  l.ir  Vmireiiati's. 


Fir  1.  SnrI.ice  Area  per  t  int  \  oinine  S,  from 
Kaylemii  Tlieors  (or  Spiieres  .ind  (roin  .Auureiiate 
rlieory  Based  on  Foint  (aint.ict  ot  I’riniarv  (’articles 
\s.  t 

creases  durinii  tlie  region  where  surface  growth 
dominates  t25  <  t  5fi  ms)  and  tlien  decreases  as 
particle  oxiilalion  oieiirs.  Both  Du,  and  r/,,  parallel 
this  trend.  The  aucrenate  nnniher  concentration  N„ 
■^I'lieralK  deere.ises  .is  uonid  occur  il  auiireuation 
continues  in  the  presence  of  either  siirtace  urowth 
or  oxidation  of  the  primary  particles.  The  averace 
miinher  ot  primary  particles  per  ai:jrrei;atc  tT  and 
die  niean-srpiare  radius  of  I'xralion  K|  increase  as 
.1  result  of  cluster-eliister  aiiureuation.  Primary  par¬ 
ticle  concentration  \,,  is  not  shown  hut  is  found  to 
remain  close  to  1  x  10 '‘^m  '  for  ail  t. 

_  riie  particle  surface  area.  S,.  calculated  nsinil 
n'.  N„.  .md  </,,  I’fF.M  is  shown  m  Fimire  t.  .Also 
(lisplated  is  the  \ahie  ol  .S,  were  it  to  he  calculatetl 
on  the  basis  that  the  particles  ari'  spheres  rather 
than  auureuates.  It  is  .ipparent  that  the  predicted 
surface  areas  ditter  lis  ,i  laetor  ot  two  to  lour.  ’I'liis 
ditierence  is  ipiite  sicmlieant  in  the  exahiation  ol 
the  surface  urowtli  r.itcs  .md  points  to  the  inipor- 
fanei-  ot  .iccoimtinii  tor  the  .lunreuate  stnietiire  ol 
the  particles 

From  the  .iimri’U.Uc  iimnher  concentration  that  is 
shown  III  Fin  I  we  caicni.ite  llic  .inureuation  rate 
K„  which  IS  defined  tlironnii  the  expression'' 

d\  .  K„  . 

— ^  - N:  22' 

dt  2 


I  lie  fitter  c(]nalion  i  .in  he  ic.irr.nined  .is 


»f 

Fl(.  ,3  The  (Qnaiitilies  D„,.  ,i‘  .md  iB'.)'  '  ss.  t  .is 
Determined  I  sinn  (.ross  Sections  tor  .Alturenates, 
\  aloes  ot  (/,,  iTE.M)  are  (rom  electron  niicroscops  ■ 


K., 


2  il 


2:3' 


.ind  therelore  the  slope  ot  the  liMst-sipiare  strainht 
line  tit  to  the  uraph  ot  the  ratio  /N„iL 's  time 
t  yields  an  .ixerane  v.ilne  ol  K„.  \. lines  ot  K„  that 
result  (rom  use  ot  the  several  \. lines  of  retractive 
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iiulex  are  mven  in  I'ahlc  II.  value  ol  x 

10~*"  i  tn’/s  at  29S  K  is  tvpieal  of  spherical  parti¬ 
cles. “The  lowest  value  ot  K„  =  79  x  U)  cm ‘/s 
quoted  in  Table  11>  appears  to  he  reasonable,  hut 
replication  ot  these  results  is  desired. 


Summary  and  Conclusions 

A  data  reduction  method  that  ’.s  based  on  the  use 
ol Optical  cross  sections  lor  augreijates  has  been  ap¬ 
plied  to  measurements  ot  local  e.xtinction  and  dif¬ 
ferential  scattering  at  three  ;ingles  by  soot  aggre¬ 
gates  in  a  laminar  ethene  diffusion  tlame.  As 
compared  to  the  Ravleigh  data  reduction  for  spher¬ 
ical  particles,  this  method  yields  higher  values  of 
volume-mean  diameter,  a  larger  surface  area  per  unit 
volume,  and  lower  values  of  the  aggregate  number 
concentration.  The  analysis  also  provides  the  mean- 
s()uare  radius  of  gyration  and  the  average  number 
of  primary  particles  per  aggregate.  Both  these 
ipiantities  increase  monotonically  with  time  as  is 
consistent  with  the  (Kcurrence  of  cluster-cluster  ag¬ 
gregation.  Frimarv  particle  diameter  is  also  ob¬ 
tained  from  the  optical  observations  and  it  has  been 
independently  measured  by  thermophoretic  sam¬ 
pling  fiillowed  by  TRM  analysis.  By  requiring  the¬ 
rms  difference  between  the  two  observations  of  pri- 
mar\  particle  size  to  be  a  minimum,  we  obtain  in¬ 
formation  on  the  allowable  values  of  refractive  in¬ 
dex.  Higher  values  of  the  real  portion  of  the 
refractive  index  give  more  plausible  moment  ratios 
of  the  size  distribution  of  the  aggregates.  The  pres¬ 
ent  results  indicate  that  a  self-consistent  interpre¬ 
tation  of  the  light  scattering  properties  of  the  soot 
aerosol  is  afforded  by  recognizing  its  aggregate 
structure.  In  particular,  significant  differences  in  the 
surface  area  and  number  concentration  are  ob- 
si-rved  for  soot  aerosols  consisting  of  aggregates  of 
primary  particles,  The.se  differences  have  important 
ramifications  in  deriving  chemical  and  physical  rate 
data  from  the  laser  scattering  extinction  experi¬ 
ment.  Furthermore,  when  combined  with  TEM 
.inalvsis  of  the  primary  particles  which  constitute 
the  aggregates,  the  present  analysis  provides  the 
basis  for  ev  aluating  the  importance  of  other  particle 
properties  such  as  refractive  index.  Thus,  this 
inethiKl  of  data  analvsis  provides  a  more  detailed 
and  comprehensive  description  of  the  development 
and  character  ot  the  sinit  aggregate  population  in 
laminar  diffusion  flames. 
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.Appendix  I  _ 

Relationship  Between  d,,  and  R"  For 
Polydisper.se  .Aggregates 

For  a  polvilispcrsion  ol  .iggregates.  the  mcan- 
sipiarc  radius  of  gvration  is  defined  bv 


^  R311)  11"  p  n> 


V  piir 


.ind  from  Kqs.  i2',  1  H  .md  i.Al',  it  follows  that 


V  pi/l>  V  ri 

\  "I 


When  D,  =  2.0.  F.q.  i.A.il  reduces  to 

H  =  cA4t 

The  value  of  the  prcfactor  k,  vv.is  found  through 
computer  simulations'  to  be  .A.S  lor  narrow  si/e 
distributions  H  —  k/.  while  for  discrete  log  normal 
distributions  ol  moderate  width  H  dcireases  to  less 
than  two  when  I),  =  I  S.  An  cvammation  ol  two 
aggregate  samples,  which  were  taken  in  the  ethene 
diffusion  llame'-  .it  a  height  :  =  1)  mm  .ibove 
burner,  vields  /„  =  l.S  and  H  ~  I  .5  to  I...  in 
general  agreement  with  calculations  based  on  .1  dis¬ 
crete  log  normal  PDF. 


.Appendix  11 

Correction  of  Observed  Flxfinction  to  Obtain  tlie 
Volumetric  Rayleigh  .Absorption  Cross  Section 

For  aggregates  ol  miKlerate  size.  (/'Hr;  .ibout  e<|ual 
to  iinitv  or  larger,  the  measured  extinction  con¬ 
tains.  in  addition  to  absorption  which  dominates  in 
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the  tast'  ot  ihf  Miuilli'i  .liZUri'Uatcs,  a  Mumluant  loii- 
fnlHitidii  ()\^  illl'  to  s(  attonmt  Tlu-  I'xttiii  tioii  can  lie 
expressed 

is.  o  d  ^  p...  V5' 

where  p^.,  is  the  ratio  ot  the  seatteriTiu  i  ross  M-elion 
to  the  absorption  eross  settion.  Itie  latter  i'  i;i\en 
ln“ 

p7,  -  f,i  n'  oj,,  i:-k'  Vfi' 

where  oj,„  the  allu-do  ot  the  iMiinarx  partnie.  is 
liixen  In 

2  ,  Kim  I 


and  the  tnnetion  >iik~  H5;  t;i'en  In 

I),  j 

-ik-  W:'  =  I  -  —  il-  H: )  \.S' 

'  )D,  ' 

Kroin  K<is,  .  1’  and  i.\2t  the  tirst  moment  ot  the  a^- 
^reijate  si/e  distribution  can  be  expressed  In 


The  ratio  ot  scattcrini!  to  extinction  p7.  can  be  cal¬ 
culated  iisimt  I'EM  oliserxations  ot  </,,  and  the  val¬ 
ues  of  Rn  tiiiind  iisimj  the  observeil  disv  nimetrv  ra¬ 
tios,  R|j  .ind  Rn  ith  s|)ecifieil  values  ot  It  and 
the  ipiantities  iT  io,,.  2  and  p,„  can  be  calcu¬ 
lated.  In  the  non-sinokinii  ethene  tlaine.  tin-  value 
ot  p,„  is  as  hiizh  as  3()‘7f  in  the  nndllame  leuion 
where  the  primarx'  particles  are  33  nm  m  iliameter 
Thus,  the  peak  values  ot  siKit  volume  traction  f\  are 
reduct*d  bv-  about  dO'T  below  the  value  which  would 
be  predicted  by  the  conventional  Ravieiuh  analvsis 
lor  spheres  in  which  a  contribution  to  extinction  bv 
scatterinc  is  not  included. 
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Abstract — Buoyancy  is  an  important  factor  in  the  dynamic  behavior  of  jet  diffusion  flames.  In  order  to 
determine  the  exact  role  that  buoyancy  plays,  a  simple  procedure  is  described  for  varying  in  isolation  the 
relative  buoyancy  force  in  stationary  laboratory  jet  diffusion  flame  experiments.  This  procedure,  which  is 
denved  from  a  theoretical  model  of  these  flames,  merely  requires  that  background  pressure  be  varied  while 
maintaining  constant  mass  flows  of  fuel  and  oxidizer  into  the  burner.  It  is  shown  that  the  sole  result  of  these 
pressure  vanations  in  the  theoretical  model  is  that  the  effective  gravitational  acceleration  acting  upon  the 
flame  varies  as  the  square  of  the  pressure.  Comparisons  are  made  between  the  structure  of  a  low  speed 
laboratory  methane/air  flame  at  various  pressures  and  the  results  of  a  direct  numerical  simulation  of  the 
same  flame  with  various  gravitational  accelerations.  Similar  evolutions  in  flame  structure  are  observed  in 
both  cases. 

INTRODUCTION 

Investigations  into  the  dynamic  structure  of  jet  diffusion  flames  have  become  more 
common  in  the  last  decade  [Yule  el  al.  (1981),  Eickhoff  and  Winandy  (1985),  Roque- 
more  et  al.  (1989),  Takahashi  el  al.  (1988),  Chen  el  al.  (1988)].  This  is  because  these 
geometrically-simple  flames  exhibit  many  features  which  are  generic  to  more  complex 
combusting  flows.  The  interactions  between  the  flame  surface  and  outer  surrounding 
vortex  structures  are  of  particular  concern,  especially  in  low  speed  and  transitional 
(buoyancy-dominated)  jet  diffusion  flames.  It  is  these  interactions  which  are  apparent¬ 
ly  responsible  for  various  forms  of  flame  flicker  [Chen  et  al.  (1988),  Roquemore  et  al. 
(1989),  Buckmaster  and  Peters  (1986)].  Since  the  outer  vortex  structures  form  in  a 
buoyancy-driven  shear  layer  outside  the  flame  surface,  it  becomes  highly  desirable  to 
vary  the  buoyancy  force  (relative  to  inertial  and  viscous  forces)  in  isolation  from  all 
other  parameters  in  order  to  observe  the  resulting  effects,  especially  on  flame  flicker. 
However,  this  has  been  difficult  and  expensive  to  do,  requiring  accelerating  test  rigs 
[e.g..  Altenkirch  et  al.  (1976)].  It  is  the  purpose  of  the  present  paper  to  outline  a  more 
practical  procedure  for  isolating  buoyancy  effects  in  experiments  involving  jet  dif¬ 
fusion  flames. 

The  procedure  to  be  described  here  is  derived  from  a  theoretical  model  for  jet 
diffusion  flames  based  on  the  flame  sheet  and  conserved  variable  approximations.  It 
will  be  shown  that  buoyancy  effects  in  this  model  can  be  isolated  by  simply  maintain¬ 
ing  constant  mass  flows  of  fuel  and  oxidizer  into  the  burner  while  varying  the 
background  pressure.  In  order  to  test  the  validity  of  this  procedure,  results  obtained 
from  numerical  solutions  of  the  equations  of  motion  are  compared  with  experimental 
results  obtained  in  a  pressure  chamber.  A  similar  evolution  in  large-scale  flickering 
behavior  is  observed  both  numerically  and  experimentally  as  pressure  increases,  albeit 
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at  different  pressure  levels  due  presumably  to  stability  differences  between  the  com¬ 
puted  and  actual  flames.  Previous  experimental  investigations  of  jet  diffusion  flames 
at  varying  pressures  have  focused  on  the  structural  effects  of  Reynolds  number 
variations  in  the  excited  case  [Strawa  and  Cantwell]  or  on  soot  formation  [Flower 
(1988)].  The  present  study  is  believed  to  be  the  first  investigation  of  these  flames  to 
focus  solely  on  isolating  buoyancy  effects.  It  is  not.  however,  the  purpose  of  this  paper 
to  thoroughly  document  these  effects  but  instead  to  open  the  door  to  future  more 
exhaustive  studies. 


THEORETICAL  BACKGROUND 

The  nondimensional  continuity,  momentum  and  state  equations  for  the  unsteady 
axisymmetric  flow  of  a  multicomponent  fuel  Jet  surrounded  by  a  coflowing  airstream 
are 


g  +  V-{pq)  =  0.  (1) 

p  ^  -I-  p(q  •  V)q  =  -Vp  -  pR/f.  +  ^  { -  V.x(V.vpq)  -I-  4/3V(pV  •  q)  (2) 
-I-  V(q  •  Vp)  -  qV^p  -f-  Vp.x(V.xq)  -  (V  •  q)Vp}, 

pTl^R^Y,  ==  po  =  K.-  (3) 


Here  p  is  density;  p  is  viscosity;  T  is  temperature;  q  =  (v,  v),  where  v  and  u  are 
velocity  components  in  an  axisymmetric  reference  frame  (r,  r);  Yj  is  the  mass  fraction 
of  species  i  with  gas  constant  R,  ;  and  e.  is  the  unit  normal  in  the  axial  direction.  The 
pressure  consists  of  a  constant  background  pressure,  po.  pltJS  ^  perturbation  pressure. 
p(r.  r,  r),  where  Po  >  >  p.  As  discussed  by  Rehm  and  Baum  (1978),  the  omission  of 
p  from  Eq.  (3)  is  entirely  appropriate  here,  resulting  in  the  suppression  of  acoustic 
waves  while  admitting  the  low  frequency  motions  due  to  buoyant  effects.  All  quan¬ 
tities  in  Eqs.  (1-3)  have  been  nondimensionalized  with  respect  to  conditions  in  the 


U  L 

coaxial  airstream  at  the  burner  inlet.  Thus  Re  —  Reynolds  number  =  — ^  and 


Ri 


=  Richardson  number  = 


ur 


where  and  are  the  entering  airstream  velocity 


and  kinematic  viscosity;  L  is  the  fuel  jet  radius,  and  g  is  gravitational  acceleration. 

In  order  to  reduce  the  complexity  of  the  problem  while  still  providing  a  good 
approximation  to  the  physical  situation,  a  flame  sheet  model  (infinite  rate  kinetics) 
will  be  employed  here.  In  addition,  the  Lewis  number  (ratio  of  thermal  to  species 
diffusivities)  is  assumed  to  be  unity.  Therefore,  conserved  variables.  P,.  for  species 
mass  fractions  and  enthalpy  can  be  utilized.  These  obey  the  simple  nondimensional 
convection-diffusion  equation 


^  +  V-(pq/?,)  =  Lv-ipDVP,). 


(4) 
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where  D  is  binary  diffusion  coefficient  and  Pe  =  Peclet  number  =  UnLIDf.  with  Df 
taken  as  the  diffusion  coefficient  of  fuel  into  nitrogen  at  the  burner  inlet.  The 
conserved  variables  are  defined  as 


/?>  = 


= 


aK.,  +  y'f  -  Y, 

oiYi  +  t;, 

^Q  +  h‘- 
2 


Here  h  is  enthalpy;  Q  is  heat  release  per  unit  mass  of  oxygen  consumed;  2 


where  v  is  stoichiometric  coefficient  and  M  molecular  weight;  subscripts / and  ox  refer 
to  fuel  and  oxidizer  (with  nitrogen  accounted  for  ia  tnc  enthalpy);  and  superscript  / 
refers  to  burner  inlet  conditions.  The  analytical  model  for  the  jet  diffusion  flame  under 
discussion  here  consists  of  Eqs.  (l-4i  along  with  appropriate  boundary  conditions.  A 
flame  sheet  exists  at  locations  where  )S|  ==  ///(aT',  +  T/).  The  solution  method  will 
be  discussed  in  the  next  secf  on.  The  remainder  of  the  discussion  here  will  focus  on 
the  overall  characteristics  of  the  model. 

I*,  is  apparent  from  Eqs.  (1-4)  that,  given  a  fixed  set  of  boundary  conditions,  the 
solution  field  is  dependent  upon  only  three  parameters:  Re,  Ri,  and  Pe.  Of  particular 
concern  here  is  the  Richardson  number,  Ri,  which  represents  the  ratio  of  buoyancy 
forces  to  inertial  forces.  By  varying  this  parameter  independently,  it  is  possible  to 
isolate  the  effects  of  a  varying  relative  buoyancy  force  (or  gravitational  acceleration) 
on  the  theoretical  jet  diffusion  flame.  If  the  model  here  is  assumed  to  be  a  good 
approximation  to  reality,  then  this  isolation  of  buoyancy  effects  could  also  be  per¬ 
formed  experimentally  if  Ri  could  be  independently  varied  in  the  laboratory.  A 
procedure  for  accomplishing  this  will  now  be  discussed. 

A  simple  method  for  varying  Ri  while  keeping  Re,  Pe  and  the  boundary  conditions 
fixed  is  merely  to  maintain  constant  mass  flows  of  both  fuel  and  air  into  the  burner 
while  varying  the  background  pressure.  Since  po  Uq  is  fixed  and  po  is  pressure  indepen¬ 
dent,  Re  will  remain  constant  under  this  procedure.  Both  Uq  and  Df  vary  inversely 
with  pressure,  thus  maintaining  a  fixed  Pe.  However,  Ri  will  vary  as  the  square  of  the 
pressure.  Since  the  other  thermophysical  properties  inherent  in  the  model  (specific 
heats,  Lewis  number)  can  be  assumed  independent  of  pressure,  this  variation  in  Ri  is 
the  only  effect  of  the  background  pressure  changes.  Thus,  the  effective  gravitational 
acceleration  will  vary  as  pressure  squared,  e.g.,  a  background  pressure  of  2  atmos¬ 
pheres  corresponds  to  4  g's.  This  relationship  is  the  basis  of  this  paper  and  will  be 
investigated  both  computationally  and  experimentally.  It  is  reasonable  to  believe  that, 
at  least  over  a  moderate  range  of  pressure  changes,  this  procedure  can  provide  a 
simple  means  of  varying  the  relative  buoyancy  force  in  these  flames  without  resorting 
to  expensive  moving  test  rigs. 


NUMERICAL  MODEL 

The  analytical  model  consisting  of  Eqs.  (1-4)  and  appropriate  boundary  conditions 
is  solved  numerically  by  means  of  a  modified  version  of  the  finite  difference  scheme 
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employed  previously  by  Davis  and  Moore  (1985)  to  simulate  cold  mixing  layer 
dynamics.  This  scheme  employs  a  variably-spaced  staggered  mesh  in  which  pressures 
are  defined  at  cell  centers  and  normal  velocities  at  cell  faces.  Quadratic  upwind 
differencing  is  used  for  convection,  and  an  explicit  Leith-type  of  temporal  differencing 
is  employed.  At  each  time  step,  a  Poisson  equation  for  perturbation  pressure  is  solved 
by  a  direct  method.  Variable  thermophysical  properties  are  calculated  as  described  by 
Mitchell  (1980),  and  no  turbulence  or  radiation  models  are  employed.  Species  binary 
diffusion  is  taken  as  fuel  into  nitrogen  inside  the  flame  and  oxygen  into  nitrogen 
outside.  All  relevant  variables  are  specified  at  the  burner  inlet  (z  =  0);  axisymmetry 
is  enforced  along  the  burner  centerline;  and  zero  or  constant  gradients  are  employed 
at  very  large  axial  and  radial  distances  (Le.,  no  shroud).  Mesh  dimensions  are  5000 
fuel  jet  radii  axially  (145  mesh  cells)  and  50  radially  (84  cells),  both  attained  by 
employing  rapidly  expanding  mesh  cell  sizes  away  from  the  flame.  Run  times  on  the 
NIST  CYBER  205  are  typically  several  hours,  with  initial  solution  fields  obtained  by 
means  of  a  steady-state  flame  simulation  [Davis  and  Moore  (1987)].  After  awhile, 
flame  flickering  commences  in  the  unsteady  computation  without  the  triggering 
required  in  the  cold  Jet  case  [Davis  and  Moore  (1985)].  The  numerical  model  Just 
described  has  also  been  used  to  study  the  dynamic  structure  of  propane/air  Jet 
diffusion  flames,  the  results  of  which  have  compared  well  with  experimental  visual¬ 
izations  [Davis  et  ai]. 


EXPERIMENTAL  APPARATUS 

The  experimental  portion  of  this  study  was  conducted  in  a  high  pressure  diffusion 
flame  burner  (Figure  1 )  designed  for  operation  at  pressures  up  to  twenty  atmospheres. 
The  coannular  burner  employed  here  consists  of  a  fuel  tube  of  1 . 10  cm  inner  diameter 
surrounded  by  a  10. 16  cm  inner  diameter  air  passage.  The  air  passage  is  partially  filled 
with  glass  beads  followed  by  a  series  of  fine  mesh  screens  to  provide  for  flow 
conditioning  of  the  incoming  air.  A  honeycomb  section  2.54cm  in  thickness  with 
0. 1 5  cm  square  cells  is  located  at  the  burner  exit  to  provide  a  uniform  exit  air  flow  field. 
The  fuel  tube,  which  extends  0.25  cm  above  the  ceramic  honeycomb,  is  also  partially 
filled  with  glass  beads  for  flow  conditioning. 

The  burner  is  housed  in  a  high  pressure  chamber  which  consists  of  two  cylindrical 
sections.  The  lower  section  is  constructed  of  20  cm  nominal  diameter  pipe  (8  inch 
diameter,  schedule  40  carbon  steel  pipe)  which  is  61  cm  in  length.  This  section  houses 
the  motorized  translational  stage  which  provides  for  vertical  movement  of  the  burner. 
Electrical  connections  as  well  as  gas  lines  for  the  air  and  fuel  are  introduced  through 
a  flange  located  on  the  bottom  of  this  section.  The  bottom  flange  is  attached  to  a  ball 
bearing  stage  which  supports  the  weight  of  the  entire  chamber.  Horizontal  movement 
of  the  chamber  can  thus  be  achieved  using  a  small  motorized  translation  stage.  Both 
the  horizontal  and  vertical  movements  of  the  burner  are  controlled  by  stepper  motors 
which  are  interfaced  to  a  laboratory  personal  computer.  For  the  present  study  only 
the  vertical  motion  capability  was  required. 

The  upper  section  of  the  chamber  contains  four  6.35  cm  diameter  windows  which 
are  1 .27  cm  in  thickness  and  are  at  the  same  axial  height  90°  apart.  These  windows 
provide  access  for  high  speed  video  recording  as  well  as  laser  diagnostic  measure¬ 
ments.  This  section  is  constructed  of  a  76.2  cm  long,  15  cm  nominal  diameter  carbon 
steel  pipe  with  a  wall  thickness  of  2.19  cm.  This  thicker  wall  pipe  is  required  in  order 
to  maintain  the  proper  material  strength  after  machining  the  holes  for  the  windows. 
The  operating  pressure  for  the  burner  is  controlled  using  manual  valves  located  on  a 
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where  D  is  binary  diffusion  coefficient  and  Pe  =  Peclel  number  =  UnLjDf.  with  Df 
taken  as  the  diffusion  coefficient  of  fuel  into  nitrogen  at  the  burner  inlet.  The 
conserved  variables  are  defined  as 


ar„,  +  Y'r  -  Yr 

+  y'r^ 


P2  = 


Here  h  is  enthalpy;  Q  is  heat  release  per  unit  mass  of  oxygen  consumed;  x  = 

where  v  is  stoichiometric  coefficient  and  M  molecular  weight;  subscripts / and  at  refer 
to  fuel  and  oxidizer  (with  nitrogen  accounted  for  in  the  enthalpy);  and  superscript  7 
refers  to  burner  inlet  conditions.  The  analytical  model  for  the  jet  diffusion  flame  under 
discussion  here  consists  of  Eqs.  ( 1  -4)  along  with  appropriate  boundary  conditions.  A 
flame  sheet  exists  at  locations  where  ^1  =  T//(aT'j  +  T/).  The  solution  method  will 
be  discussed  in  the  next  section.  The  remainder  of  the  discussion  here  will  focus  on 
the  overall  characteristics  of  the  model. 

It  is  apparent  from  Eqs.  (1-4)  that,  given  a  fixed  set  of  boundary  conditions,  the 
solution  field  is  dependent  upon  only  three  parameters:  Re,  Ri.  and  Pe.  Of  particular 
concern  here  is  the  Richardson  number,  Ri,  which  represents  the  ratio  of  buoyancy 
forces  to  inertial  forces.  By  varying  this  parameter  independently,  it  is  possible  to 
isolate  the  effects  of  a  varying  relative  buoyancy  force  (or  gravitational  acceleration) 
on  the  theoretical  jet  diffusion  flame.  If  the  model  here  is  assumed  to  be  a  good 
approximation  to  reality,  then  this  isolation  of  buoyancy  effects  could  also  be  per¬ 
formed  experimentally  if  Ri  could  be  independently  varied  in  the  laboratory.  A 
procedure  for  accomplishing  this  will  now  be  discussed. 

A  simple  method  for  varying  Ri  while  keeping  Re,  Pe  and  the  boundary  conditions 
fixed  is  merely  to  maintain  constant  mass  flows  of  both  fuel  and  air  into  the  burner 
while  varying  the  background  pressure.  Since  Po  Uq  is  fixed  and  Pq  is  pressure  indepen¬ 
dent,  Re  will  remain  constant  under  this  procedure.  Both  and  Df  vary  inversely 
with  pressure,  thus  maintaining  a  fixed  Pe.  However.  Ri  will  vary  as  the  square  of  the 
pressure.  Since  the  other  thermophysical  properties  inherent  in  the  model  (specific 
heats.  Lewis  number)  can  be  assumed  independent  of  pressure,  this  variation  in  Ri  is 
the  only  effect  of  the  background  pressure  changes.  Thus,  the  effective  gravitational 
acceleration  will  vary  as  pressure  squared,  e.g.,  a  background  pressure  of  2  atmos¬ 
pheres  corresponds  to  4  g  s.  This  relationship  is  the  basis  of  this  paper  and  will  be 
investigated  both  computationally  and  experimentally.  It  is  reasonable  to  believe  that, 
at  least  over  a  moderate  range  of  pressure  changes,  this  procedure  can  provide  a 
simple  means  of  varying  the  relative  buoyancy  force  in  these  flames  without  resorting 
to  expensive  moving  test  rigs. 


NUMERICAL  MODEL 

The  analytical  model  consisting  of  Eqs.  (1-4)  and  appropriate  boundary  conditions 
is  solved  numerically  by  means  of  a  modified  version  of  the  finite  difference  scheme 
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FIGURE  1  ‘  -'limatic  representation  of  high  pressure  diffusion  flame  burner. 

flange  mounted  to  the  top  of  the  pressure  chamber.  Steady  operating  pressures  over 
a  range  of  one  to  twenty  atmospheres  have  been  maintained  in  this  system. 

Gas  flows  are  measured  using  a  mass  flow  metering  system.  For  the  fuel  flow  rates, 
a  mass  controller  is  utilized  which  maintains  the  mass  fuel  flow  rate  constant  regard¬ 
less  of  the  chamber  pressure.  The  air  flow  rate  is  monitored  using  a  mass  flow  meter 
and,  thus,  the  air  flow  has  to  be  manually  adjusted  to  maintain  a  constant  mass  flow 
rate  for  the  air  as  the  chamber  pressure  is  varied. 

In  the  present  study  methane  was  burned  with  air  supplied  from  an  in-house 
compressor  capable  of  providing  source  pressures  of  300  psig  and  up  to  300  SCFM 
of  air  to  the  burner.  These  capabilities  exceed  the  pressure  and  flow  rate  requirements 
for  this  study.  Prior  to  introducing  the  air  into  the  burner  it  is  passed  through  filters 
to  remove  particles  and  moisture.  Methane  is  supplied  from  high  pressure  cylinders 
with  a  purity  of  99%. 

High  speed  video  recording  of  the  flame  structure  was  obtained  using  a  Spin  Physics 
high  speed  video  camera  operated  at  a  500  Hz  framing  rate.  Observations  of  the  flame 
were  made  through  one  of  the  four  windows  located  in  the  chamber  and  provided  a 
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TABLE  I 

Experimental/Numerical  Comparison  of  Tip-Cutting  Heights  and  Frequencies 


Pressure.  p„  (atm) 

Tip-Cutting  Height.  H,  tcm) 

Frequency.  /  ( Hz) 

Experiment 

Computation 

Experiment 

Computation 

Experiment 

Computation 

1.0 

0.50  {Ri  =  0.34) 

H  =  »' 

H  =  10' 

Stead  V 

Steady 

1.5 

0.75  {Ri  =  0.76) 

8 

9 

18.5 

15 

:.o 

1.20  {Ri  =  1.94) 

6 

7 

15 

15 

3.0 

1.41  (/?/  =  2.70) 

4 

6 

16 

15 

'H  denotes  full  height  of  steady  flame 


5.1  cm  open  aperture  viewing  region.  The  flame  structure  images  were  obtained  using 
the  soot  luminosity  of  the  flame.  Since  the  flame  dimensions  typically  exceeded  the 
5.1  cm  viewing  area  provided  by  the  window,  the  burner  was  raised  or  lowered  to  an 
appropriate  region  of  interest  for  study.  In  some  cases  a  series  of  video  sequences  were 
required  to  characterize  the  entire  flame.  The  Spin  Physics  camera  also  provides  the 
capability  to  play  back  the  recorded  images  at  slow  speeds  or  a  single  frame  at  a  time. 
This  capability  was  used  in  the  analysis  of  the  recorded  images  which  follows. 


RESULTS  AND  DISCUSSION 

As  a  test  of  the  aformentioned  procedure  for  varying  the  relative  buoyancy  force  in 
jet  diffusion  flames,  a  coordinated  series  of  computations  and  experiments  has  been 
carried  out.  The  selected  configuration  was  a  methane/air  diffusion  flame  with  fuel 
and  air  burner  inlet  velocities  at  one  atmosphere  of  10  and  20  cm/sec,  respectively,  and 
a  fuel  jet  radius  of  0.55  cm.  This  flame  appears  virtually  steady  at  atmospheric 
pressure  while  entering  into  a  simple  tip-cutting  mode  at  higher  pressures.  This  mode 
is  quite  regular  in  frequency  and  extremely  reproducible  in  terms  of  the  flame 
structures  which  result.  The  flame  is  observed  to  periodically  pinch  off  at  a  particular 
axial  height,  with  the  cut-off  tip  being  convected  upward  while  the  lower  region 
recedes.  The  lower  part  of  the  flame  subsequently  increases  in  length  and  the  cycle 
repeats.  When  viewed  in  real  time,  this  periodic  tip-cutting  appears  as  a  flickering 
flame.  The  large-scale  features  of  this  flickering  phenomenon  are  highly  dependent  on 
pressure  level,  thus  allowing  use  of  simple  flame  visualizations  to  make  comparisons 
between  experiment  and  computation.  Due  to  the  low  flow  velocities  involved, 
buoyancy  is  clearly  an  important  factor  in  this  flame's  behavior,  while  small-scale 
turbulence  is  almost  certainly  not.  The  low  level  of  soot  production  minimizes 
radiative  losses,  and  the  Lewis  number  can  reasonably  be  assumed  as  unity.  Thus,  for 
all  these  reasons,  this  flame  was  deemed  to  be  a  suitable  candidate  for  the  present 
study. 

Experiments  involving  this  methane/air  flame  have  been  carried  out  at  pressures 
ranging  from  1  to  8  atmospheres.  In  these  experiments  the  mass  flow  rates  of  the  fuel 
and  air  were  maintained  constant  as  the  chamber  operating  pressure  was  increased. 
Thus,  the  volumetric  flow  rates  and  linear  velocities  at  the  burner  exit  were 
proportionately  reduced.  Corresponding  computations  (with  slug  flow  burner  inlet 
velocity  profiles)  have  been  performed  by  simply  varying  the  gravitational  accelera¬ 
tion  in  the  model  as  the  square  of  this  background  pressure;  i.e..  Ri  x.  pi,  where 
Ri  =  1 .35  for  =  1  atm.  The  results  are  summarized  in  Table  I.  with  the  evolution 
in  flame  structure  as  a  function  of  (or  Ri)  shown  in  Figure  2.  This  figure  compares 
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recorded  images  of  the  experimental  flame  at  four  pressures  with  computed  flame 
sheets  (solid  squares)  and  their  surrounding  vortex  structures  visualized  via  non- 
dimensional  isovorticity  contours.  Note  that  the  (nondimensional)  scales  in  the 
experimental  and  computational  visualizations  shown  here  have  not  been  matched 
and  that  the  computational  visualizations  extend  lower  than  the  experimental  ones. 
It  can  be  seen  from  Figure  2  and  Table  I  that  the  experimental  and  computational 
evolutions  of  flame  structure  with  increasing  pressure  are  strikingly  similar.  The 
computational  evolution,  however,  occurs  at  lower  pressures,  a  subject  which  will  be 
discussed  subsequently.  The  particular  computational  pressures  chosen  for  these 
comparisons  were  selected  so  as  to  match  the  experimental  flame  visualizations  as 
closely  as  possible,  with  phase  angles  arbitrarily  picked.  The  flame  structure  seen  at 
low  pressure  (1  atm  experimentally;  0.5  atm  computationally)  in  Figure  2a  is  essen¬ 
tially  steady,  although  there  is  slight  tip  wavering  in  both  cases.  Some  weak  outer 
vortices  can  be  seen  near  the  flame  tip  in  the  computational  visualization.  At  the 
higher  pressures  (1.5  atm  experimentally;  0.75  atm  computationally)  of  Figure  2b. 
tip-cutting  has  begun  in  both  cases.  The  tip-cutting  height.  //,,  and  frequency,/,  are 
indicated  in  Table  I  for  both  the  experimental  and  computational  cases.  The  external 
vortex  structures  have  strengthened  and  moved  downward  as  compared  with  Figure 
2a.  This  is  due  to  the  increased  buoyancy-induced  velocity  gradients  at  this  pressure, 
with  the  effective  gravitational  acceleration  here  being  greater  than  that  for  Figure  2a 
by  a  factor  of  2.25.  At  the  next  higher  pressure  (2  atm  experimentally;  1.2  atm 
computationally),  the  cut-off  tip  has  assumed  a  basically  triangular  shape  with  a 
cusped  region  flaring  out  from  its  base.  Both  the  experimental  and  computational 
tip-cutting  heights  have  decreased  from  Figure  2b,  with  the  outer  vortices  also  moving 
down  commensurately.  The  experimental  frequency,  but  not  the  computational  one, 
has  dropped  slightly.  The  situation  at  even  higher  pressure  (3  atm  experimentally; 
1.41  atm  computationally)  is  illustrated  in  Figure  2d.  The  base  of  the  cut-off  tip  has 
straightened  out  to  produce  a  highly  triangular  shape.  The  tip-cutting  heights  have 
dropped  again,  the  frequencies  are  about  the  same,  and  the  outer  vortex  structures  are 
somewhat  lower.  It  is  clear  from  both  Table  I  and  Figure  2  that  the  experimental 
evolution  of  flame  structure  with  increasing  pressure  compares  well  with  the  evolution 
computed  via  Richardson  number  increases,  thus  tending  to  confirm  the  basic  validity 
of  the  buoyancy-isolation  procedure  at  issue  here. 

Although  comparisons  between  numerical  and  experimental  results  have  been 
made  for  a  pressure  range  of  one  to  three  atmospheres,  high  speed  video  recordings 
have  been  obtained  for  pressures  as  high  as  eight  atmospheres.  Comparisons  are 
limited  to  the  lower  range  of  pressures  because  three-dimensional  effects  become 
significant  at  pressures  near  and  above  four  atmospheres,  thus  rendering  comparisons 
with  the  axisymmetric  model  of  little  value.  This  loss  of  axisymmetry  in  the  experi¬ 
mental  flame  is  most  evident  in  the  evolution  of  the  upper  section  of  the  flame  after 
a  "tip-cutting"  event  has  occurred.  Under  all  pressure  conditions  where  flickering  is 
observed,  the  luminous  cut-off  region  is  observed  to  disappear  with  time.  This  is  due 
to  the  oxidation  of  the  soot  particles  which  are  responsible  for  the  observed  lumino¬ 
sity.  For  the  pressure  range  of  one  to  three  atmospheres,  this  process  is  observed  to 
proceed  in  a  regular,  symmetric  manner  with  the  lowest  regions  of  the  cut-off  tip 
oxidizing  first.  As  the  pressure  is  raised  to  four  atmospheres  and  above,  the  rapid 
oxidation  process  is  observed  to  favor  one  side  of  the  cut-off  tip  and  to  proceed 
diametrically  upward  and  across  this  upper  flame  structure.  At  pressures  between  five 
and  six  atmospheres,  the  upper  flame  structure  is  actually  split  into  two  or  four 
sections  depending  on  the  chamber  pressure.  In  each  case,  these  observations  are 
highly  reproducible  from  cycle  to  cycle  and  clearly  demonstrate  the  strengthening  of 
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FIGURE  2  Comparisons  between  experimentally  recorded  images  and  computational  flames  (solid 
squares  with  surrounding  isovorticity  contours)  at  various  values  of  background  pressure.  Length  scales 
are  in  terms  of  fuel  jet  radii. 
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azimuthal  instabilities  with  rising  pressure.  This  increasing  instability,  however,  is 
precisely  the  type  of  behaviour  that  would  be  expected  as  Ri  increases,  as  found  by 
Subbarao  (1989)  in  experimental  studies  of  buoyant  helium  jets.  It  is  also  noted  that 
the  addition  of  full  three-dimensionality  to  the  model  employed  here  would  in  no  way 
affect  the  conclusions  previously  reached  concerning  the  relationship  between  back¬ 
ground  pressure  variations  and  the  relative  buoyancy  force. 

Finally,  the  lower  pressures  (or  effective  gravitational  accelerations)  found  in  the 
computational  evolution  in  Figure  2  require  some  explanation.  This  is  probably  due 
to  the  use  of  the  flame  sheet  approximation  and  the  lack  of  any  radiation  in  the  model, 
thus  resulting  in  higher  temperatures  computationally  than  experimentally.  This,  in 
turn,  should  lead  to  higher  velocities  and  sharper  gradients  in  the  computed  flame 
than  in  the  experimental  one.  thus  rendering  the  former  more  unstable.  Thus  it  should 
require  less  effective  gravitational  acceleration  to  trigger  a  particular  stage  of  the 
structural  evolution  process  computationally  than  experimentally.  Therefore,  com¬ 
parisons  between  computational  and  experimental  results  at  the  same  pressure  would 
not  be  particularly  meaningful.  This  can  be  seen  by  comparing  the  computational 
visualization  in  Figure  2d  with  the  experimental  visualization  in  Figure  2b.  which  are 
both  at  approximately  the  same  pressure.  There  is  clearly  little  similarity  between 
these  two  visualizations,  yet  a  great  deal  of  similarity  is  observed  when  comparisons 
are  made  at  appropriate  differential  pressures. 

CONCLUSIONS 

A  simple  procedure  has  been  described  for  isolating  buoyancy  effects  in  jet  diffusion 
flame  experiments.  All  that  is  necessary  is  that  background  pressure  be  varied  while 
maintaining  constant  mass  flows  of  fuel  and  oxidizer  into  the  burner.  A  theoretical 
model  for  these  flames  indicates  that  this  procedure  isolates  Richardson  number  (or 
relative  buoyancy  force)  as  the  only  variable  parameter,  which  is  equivalent  to  varying 
the  gravitational  acceleration.  A  joint  series  of  experiments  and  computations  involv¬ 
ing  a  pressurized  low  speed  mcthane/air  diffusion  flame  has  provided  strong  evidence 
to  support  the  theory. 

It  is  clear  that  much  work  needs  to  be  done  in  order  to  clarify  the  important  role 
that  buoyancy  plays  in  jet  diffusion  flame  dynamics.  The  results  presented  here  only 
open  to  door  to  future  studies,  which  can  now  be  performed  with  reasonable  effort 
and  cost-effectiveness.  Future  research  is  also  needed  in  order  to  determine  the  limits 
of  usefulness  of  the  buoyancy-isolation  procedure  described  here.  Clearly  use  of  this 
procedure  to  attain  very  small  effective  gravitational  accelerations  appears  unreason¬ 
able.  as  chemical  kinetics  could  certainly  not  be  ignored  at  the  low  pressures  required. 
Additionally,  enhanced  soot  production  at  low  g  [Bahadori  el  al.  (1990)]  would  not 
be  simulated  by  this  pressure  variation  technique.  However,  relaxation  of  the  require¬ 
ment  that  Lewis  number  be  unity  appears  feasible  since  the  Prandtl  and  Schmidt 
numbers  can  both  be  considered  as  independent  of  pressure.  Finally,  since  the  basic 
model  described  here  is  not  dependent  on  a  particular  geometry,  other  types  of 
diffusion  flames  in  different  burner  configurations  could  also  be  considered  as  can¬ 
didates  for  buoyancy-isolation  experiments.  Thus,  the  possibilities  for  future  research 
efforts  in  the  area  of  buoyancy  effects  in  diffusion  flames  appear  rather  substantial. 
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